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SEVEN-YEAR PERIOD 
M., G, Chilikin and I, I. Petrov 


The far-reaching plan for the development of the 
national economy of the USSR in the period 1959-1965, 
put forth at the Twenty-first Congress of the Communist 
Party of the Soviet Union, envisions an all-around in- 
crease in the productivity of national labor on the basis 
of a strengthening of its electrical sinews and the automa- 
tion of production, The electrical motive power avail- 
able determines the degree to which production is mech- 
anized, since the basic and principal form of the prime 
mover for operating machinery is the electric drive. The 
productivity of labor depends directly on its electrifica - 
tion (i.e., on the electric motive power available to it): 
over the past ten years, the productivity of industry in the 
USSR increased by about the same percentage as its elec - 
trification. 

By an electric drive we understand a complex device 
which transforms electrical energy to mechanical energy, 
and which provides electrical control of the transformed 
mechanical energy. In this connection, the future devel- 
opment of electric drives must be directed to the imple- 
mentation of all forms of motion without the use of any 
mechanical means in achieving these purposes (rotary and 
reciprocating motion; prolonged discontinuous variable 
and tracking motion; with fixed speed, and with speeds 
which change sharply, which vary continuously in accord- 
ance with a given law, or which change arbitrarily). The 
development of electric drives must provide for the broad- 
ening of their functions in the automatic control and regu- 
lation of the operating machines of technological aggre - 
gates (control of the starting, reversing, and braking of 
operating machines; maintaining constancy of speed, 
torque and power; program control of operating machines, 
etc.), 

The science of electric drives derives basically from 
the work of Soviet scientists and, today, this branch of 
technical knowledge has taken definitive shape in the 
form of an independent section of applied electrical 
engineering. In the Soviet Union there has been created, 
and introduced into industry, a system of electric drives 
of a high technological level. However, the total volume 
of research work carried out by research and planning in- 
stitutes, by higher educational institutions, by the design 
offices of enterprises and other organizations is still in- 
sufficient, There remain unsolved many problems con- 
cerning electric drives which are of great significance for 
the country’s national economy. 

Among the unsolved problems, the first that should 
be mentioned is the creation of controlled economical ac 


SCIENTIFIC AND ENGINEERING PROBLEMS OF AUTOMATED 
ELECTRIC DRIVES DURING THE CURRENT 


electric drives. Of particular significance here are the 
questions of developing simple and reliable systems of 
frequency control for ac electric motors. and, in particu- 
lar, of creating ac supplies with _ controlled frequency. 
Here, special attention must be given to the development 
and investigation of static~frequency transformers using 
ionic, electronic, and semiconductor technologies. 

The situation is somewhat better with the develop- 
ment of controlled ac electric drives with choke control. 
In the last few years there have been created, and intro- 
duced into industry, a number of such electric drives, and 
the scientific and engineering methods for designing them 
have been adequately developed, However, the domain 
of possible applicability of electric drives with choke con- 
trol has been somewhat overestimated, which sometimes 
leads to incorrect choices in selecting the type of electric 
drive for productive machines, Thus, for example, on in- 
sufficient grounds, a choked electric drive was used for 
excavators of medium productivity, which led to their un- 
economical and unreliable performance. Special resolu- 
tions of the State Scientific-Engineering Committee of 
the Soviet Ministers of the USSR recommended that the 
production of excavators with such electric drives be halted, 
and that they be replaced by controlled dc electric drives, 
The attempt to use choked electric drives for materials- 
handling vehicles was also adjudged unsuccessful, Seem - 
ingly, it is necessary to find reasonable limits of applica - 
bility of choke-controlled electric drives, taking into ac~- 
count their relatively low energy efficiency for tight and 
prolonged control of production-machine speed, and the 
peculiarities of these machine*s load characteristics, 

Also requiring further investigation are the questions 
of the dynamics of choked electric drives, with the non- 
linearities of their individual elements taken into account, 
and also the question of electric-drive heating as the load 
varies cyclically, 

In connection with the intensification of productive 
processes, the necessity arises for developing electric 
drives with asynchronous short-circuited electric motors 
designed for production mechanisms with partial starts 
and reverses, However, this question remains unsolved, 

In the near future there must be created, based on normal- 
ly manufactured asynchronous short-circuited electric 
motors, a modification of the electric motor with a heat- 
insulating base and with an extended motor, designed for 
a high frequency of switching in, and there should alsobe 
developed a methodology for designing electric drives 
with partial starts and reverses, With this, particular at- 
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tention should be given to the transient responses in elec - 
tric drives, and to the question of their heating. 

There has been an improvement recently in the status 
of the investigations into the operating properties of ac 
electric drives supplied from autonomous supplies of com- 
mensurable power, Even here, however, there remain 
many problems in connection with increasing the power 
of asynchronous electric motors supplied from individual 
generators. The questions as to the dynamics of the pro- 
cesses transpiring in such electric drives deserve particular 
attention, 

It is also going to be necessary to investigate the ques- 
tions related to the use of synchronous electric motors for 
uncontrolled (operating at constant speed) and controlled 
(operating with variable speed) industrial electric drives. 
The transition to circuits with direct starting of synchro- 
nous motors and starting with permanently connected ex- 
citers led to a certain broadening of the domain of ap- 
plicability of synchronous motors, However, despite their 
specific properties and their high energy efficiency, syn- 
chronous motors have still not been as widely used as they 
should be, Still, a wide dissemination of synchronous elec - 
tric motors for drives of operational machines can signifi - 
cantly increase the power factor of electrical stands and 
provide a high rigidity of electric-drive mechanical char- 
acteristics within a wide range of speeds, The use of 
motors with rotary stators permits the problem of synphas- 
ing a synchronous shaft to be solved, etc. There might, 
therefore, be great interest in works devoted to elucidat- 
ing new domains of applicability of synchronous electric 
motors, and to the investigation of new automatic con- 
trol systems for synchronous electric drives, operating both 
with constant and, especially with variable speeds of 
rotation, In particular, interest inheres in the questions of 
the economic speed control of synchronous electric motors 
with frequency control, and in questions of automatic ex- 
citation control of electric motors for drives with impact 
loads, Also necessary are further investigation of the tran- 
sient responses in synchronous electric drives and the de- 
velopment of a scientific method of calculating these pro- 
cesses, 


Many problems in the domain of controlled dc elec - 
tric drives also require solution, In the widely used “mo- 
tor-generator™ (M-G) systems, only de motors with par- 
allel excitation are used, these being less liable to torque 
overheating than motors with series and with hybrid ex- 
citation, The use of power semiconducting rectifiers, con- 
nected in the series-excitation winding circuit, permits 
one to create new types of M-G systems with improved 
dynamic characteristics, Analogous systems can also be 
created by supplying the motor from controlled mercury - 
arc rectifiers (CMR). 

Significant success has been achieved in the devel- 
opment and dissemination of electric drives with ionic 
transformers (for controlled nonreversible drives), An es- 
sential improvement of the dynamic characteristics of 
an M-G drive system has been attained by the use of con- 
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trolled mercury-arc rectifiers for the excitation of the 
system’s machine. However, the development and dis- 
semination of such systems have been carried out extreme- 
ly slowly. 

Great difficulties ensue when a CMR-M system is 
used for reversible electric drives, In this connection, 
great interest inheres in the question of further develop- 
ment and investigation of an automatic control system 
for an ionic electric drive, where the ionic transformer 
must be considered as an element of the over-all contro] 
system. These systems must possess a high energy factor, 
permit operation in the rectifying and inverting modes, 
be fast-acting, and have little control power. 

Development and investigation are required for {onic- 
drive systems, increasing the power factor for tight speed 
control, and also for systems for the automatic controlof 
the power factor of ionic aggregates for the purpose of 
providing the very best conditions for supplying electricity 
to enterprises, 

Recently, magnetic amplifiers have been finding in- 
creasing application in electric drives, They are used in 
"magnetic amplifier — motor” electric drive systems, in 
M-G systems for generator excitation control, for linecon- 
trol by mercury -arc rectifiers, and also as industrial am- 
plifying links, 

The employment of semiconductor amplifiers has 
also been broadened, these being used advantageously in 
the form of preamplifier stages in conjuction with rotary 
and magnetic amplifiers. However, the domain of ap- 
plication of these new types of amplifiers has still not 
been defined sufficiently clearly. It is necessary to carry 
out a resolute replacement of rotary amplifiers in electric- 
drive systems by the more reliable and long-lasting mag- 
netic and semiconductor amplifiers, 

Simultaneously, it is necessary to work on the crea- 
tion of power semiconductor rectifiers for *semiconduc- 
tor amplifier — motor® electric-drive systems, and also for 
the control of generator excitation in M-G systems, The 
new means of automation (rotary and magnetic ampli- 
fiers, semiconductor and electronic devices) allow the 
creation of highly efficient electric-drive control systems, 
providing the maximum usage of the-machines, the ac- 
curate maintaining of a given trajectory of motion, etc, 
Of particular significance here are the so-called optimal 
systems in which, given limitations on trajectory, speed, 
acceleration, or heating, there is attained a minimum 
time of execution of the drive mechanism's operating 
function, 

The investigation of the transient responses in such 
essentially nonlinear systems presents significant difficulty. 
Today we are completely capable of solving these com 
plicated tasks by the method of mathematic modeling 
using analog and digital computers, However, for the 
creation of automated electric-drive systems, these ca- 
pabilities are used very weakly. This is explained, on the 
one hand, by the insufficient quantities of such machines 
at the disposal of research and planning organizations and, 
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on the other hand, by the shortage of scientific and engi - 
neering cadres familiar with the methods of mathematical 
modeling. The preparation of such cadres is unsatisfac - 
tory in the higher institutions of learning which turn out 
engineers with the “electric drive” specialty. Until now, 
all the textbooks on electric drives have been based on 
the classical methods of investigation and do not expound, 
in adequate measure, the more modern methods, includ- 
ing the methods of mathematical and physical modeling, 
The laboratory courses in the educational institutions do 
not provide work with computers, In connection with 

this, the important problems requiring solution are: equip- 
ping research and educational institutes, planning -design- 
ing, and other organizations with analog and digital com- 
puter-models; improving the knowhow of the scientific 
and engineering cadres in the domains of mathematical 
and physical modeling; further developing of research 
work on the study of transient responses in complex sys- 
tems of automated electric drives, constructed with the 
use of the new means of automation, and on the synthesis 
of the modern, essentially nonlinear, systems which pro- 
vide optimal electric-drive transient responses; modern- 
ization of the university courses in *electric drives* and 
the related disciplines in the direction of developing in 
them the more modern methods of theoretical and experi - 
mental investigation of electric drives; the creating, in 
the near future, of manuals and textbooks on electric 
drives which are designed to make the transition from 

the older methods of studying electric-drive dynamics to 
the newer ones, including the methods of mathematical 
and physical modeling, 

Certain problems in the domain of multimotor elec - 
tric drives are also not solved, Thus, further investigation 
is required of electric shaft systems with variable speed 
ratios of the individual elements, for mass-production- 
line electric drives, It is necessary to develop a theory 
of the nonsteady-state processes in electric drives with 
friction coupling between the motors, and to investigate 
questions of oscillation damping in electric shaft systems 
with frequency transformers and dual-source machines. 
There has yet to be created a generalized theory for multi - 
motor electric drives with electrical connections, provid- 
ing a constant tension on the material being processed, for 
aggregates with continuous technological processes, 

Insufficient attention has been given to the produc- 
tion, and to the use in automated electric drives, of the 
various electromagnetic clutches, including those with 
ferromagnetic packing, Nonetheless, their use could lead 
to the creation of a rational electric-drive system, partic - 
ularly in cases where it is necessary to have smooth start- 
ing of the drive, acceleration of large inertial masses, 
control of torque, speed, etc, They can also play the 
roles, in automated electric drive systems, of brake dy- 
namometer, position indexer, torque limiter, etc. Today, 
electromagnetic clutches are comparatively widely used 
only in machine-stand construction, It is necessary to 
translate the experience in machine -stand construction to 


other branches of the national economy, and to carry out 

a number of investigative tasks, directed toward the ex- 
plication of the possibilities of constructing rationalelectric- 
drive systems with electromagnetic clutches for diverse 
productive mechanisms, 

During many years, various branches of industry, and 
particularly the machine -building industry, have put for- 
ward the problem of creating electric drives with rotary - 
reciprocating and vibratory motion, Unfortunately, this 
problem is still unsolved for practical purposes, Research 
organizations must necessarily increase the exploratory 
work in the domains of creating reciprocating electric 
motors, and of constructing autooscillatory systems with 
simple control systems for the amplitude and frequency of 
the vibrations, 

Of great interest in both the theory and the practice 
of electric drives might be an investigation of systems 
containing elastic links, Today, all system elements are 
considered to be ideally solid bodies, and therefore the 
phenomena and processes in which elastic deformations 
play a part are not taken into consideration. The deforma - 
tions of elastic links lead, in many cases, to significant 
errors, It would thus be of great advantage to develop 
general methods of analyzing automatic-control systems 
containing elastic links, and to solve the problem of syn- 
thesizing electric drives and their control systems for 
mechanisms with elastic links, 

Despite the significant development of this theory in 
recent years, certain questions in the theory of electric- 
drive heating are still unsolved, Thus, further investiga - 
tion is needed of the heating processes of ac motors in 
normal condtions of use, including the asymmetric mode 
of supply of three-phase asynchronous electric motors, 
Also requiring further development are the questions of 
electric-motor heating in transient responses and for cycli- 
cal loads, with account being taken of insulation aging. 
Work should also be continued on the direct use of heat~ 
conduction equations for the thermal computations for 
electric-drive motors, 

There are also many questions in the domain of 
microdrives to be resolved, This domain is little studied 
but, at the same time, this energetically evolving. In the 
next few years, millions of electric motors will be re- 
quired for microdrives, In connection with this, there is 
necessary a further development of the questions of cre- 
ating economical systems of actual microdrives, of in- 
creasing their reliability and noiselessness in use, andalso 
the accuracy and speed of response of microdrives as used 
in the devices of automation. 

The modern electric drive plays a very large role in 
systems of automated productive processes, The partic- 
ularly essential question is that of program control of the 
electric drives of individual productive mechanisms and 
aggregates, Today, a certain measure of success has been 
achieved in the introduction of program-control systems 
for a number of metal-working stands, in the screw-down 
devices of reversible rolling mills, in the electric drives 
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of charging systems for blast furnaces, etc. However, the 
general state of work in the areas of creating and dissemi - 
nating program control systems is unsatisfactory, In these 
systems, the new means of automation - semiconducting 
and magnetic elements - are used with inadequate com- 
pleteness, and one finds little reflection of the recent 
achievements in the area of creating self-adjusting and 
self-learning systems, Asa rule, the programmed systems 
in use today have a rigidly given program, which does 
not always answer to the required technological-econom - 
ical indicators of the process under conditions of noncon- 
stant external disturbances, The use of devices providing 
automatic program correction in correspondence with 
changes in an aggregate’s external operating conditions 
would permit the creation of an automatic-control system 


which would establish the most efficient mode of operation 


A significant obstacle to the widespread introduction 
of automatic control by mechanisms supplied with elec - 
tric drives is the inadequacy of existing position trans - 
ducers of the machines’ individual organs and of the ma- 
terial to be processed, and also the inadequacy of manu- 
factured relay-contact apparatus which allows only a lim- 





ited number of switchings in, Thus, the number of con- 
tact operations on am automatic machine -stand line reaches 
tens of thousands, and sometimes even hundreds of thou- 
sands, per hour, which sharply reduces the reliability of 
electric control systems, In this connection, great im- 
portance inheres in the development of new types of con- 
tactless automatic control devices which use highly reli- 
able semiconducting and magnetic elements. 

Insufficient use is made, in control schemes with 
drives in production lines, of electric apparatus used for 
communciations devices, It is necessary to practice a 
wider use of this apparatus which has small dimensions 
and increased reliability. 

For the purpose of accelerating the large-scale use 
of modern automatic-control systems by electric drives, 
it becomes necessary to accelerate the utilization, in the 
electrical industry's standard assemblies, of new control 
systems, including digital and programmed devices, and 
also complex control systems for standard electric drives, 

The enumeration of problems in the domain of elec- 
tric drives which has been given here is far from exhaus- 
tive, and contains only the first-priority problems, 
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CHARACTERISTICS 
A. A. Fel'dbaum 


(Moscow) 


STATISTICAL THEORY OF GRADIENT SYSTEMS OF AUTOMATIC 
OPTIMIZATION FOR OBJECTS WITH QUADRATIC 


(Translated from: Avtomatika i Telemekhanika, Vol, 21, No. 2, pp. 167-179, February, 1960) 


Original article submitted June 12, 1959 


The paper formulates the problem of investigating gradient systems of automatic optimization in the 
presence of random factors, Methods are given for determining the steady-state error and estimating the 
duration of the transient response for gradient systems with several input variables, The results are illus- 
trated by examples of systems with one and with two variables, 


INTRODUCTION 


The near future will be marked by a significant de- 
velopment of self-adjusting (adaptive) systems, Great 
value,therefore,inheres in the development of a theory of 
automatic serach, which is the basic process in such sys- 
terns, and which differs essentially from the process of 
error liquidation in automatic-control systems. One of 
the classes of adaptive systems which promises to be of 
importance practically is that of automatic optimization 
systems, We present below the theory of one of the forms 
of such systems, 

The problem of an automatic optimization system 
consists of automatically finding and maintaining a min- 
imum of some quantity y with, generally speaking, addi - 
tional conditions imposed, The quantity y is a function 
of the variables u;: 


Y = y (Uy, Ug... Um). (1) 

The process of automatic search for a minimum is 
slowed when there is random noise present, and the min- 
imum itself will be determined, and maintained, with 
some error under these circumstances, The duration of 
the transient response of the search process, and the steady- 
state error, are the fundamental symptomatic indicators 
of the system's process, 

The most general problem consists of synthesizing 
such a controlling portion which will implement the best 
(in some sense) process of automatic search, We shall 
consider below the more restricted problem of determin - 
ing the indicators of the process for a given controlling 
portion, and also of choosing the optimal values of this 
postion’s parameters, 

In the automatic serach for a minimum, tentative 
movements are made — “experiments® with the object, 
as it were — and thereafter the operational movement is 
made on the basis of an analysis of the results of the ten- 
tative trial motions, Various combinations of operationa 
and trial movements are possible, In the sequel, we con- 
sider systems of the discrete type in which, during onecy- 
cle, a series of trial steps is taken, and then followed by 





an operational step; in the following cycle, there isagain 
taken a series of trial steps followed by an operational 
step, etc, 

Figure 1 shows an example of an automatic optimiza- 
tion system. Here, C is the controlling portion; at the 
output of part Fy of object 0 there is the quantity y which 
must be minimized, It is joined, in part F of the object, 
by an error z, At the object's input, the controlling quan- 
tities uj are joined by random-noise quantities, z,.;. The 
sums thus obtained 

Xi = Uji + Zoi (2) 
are applied to the input of part F, of the object, In carry- 
ing out trials on the ith input variable, one adds to the 
quantity xj the quantity +o, in the first trial interval 
and -&%,; in the second trial interval, Thus, a trial con- 
sists of two intervals, If the total number of input vari- 
ables is m, the entire cycle is composed of 2m intervals, 
We shall assume that the variations of the object's input 
quantities require the expenditure of no time, but that the 
carrying out of an "experiment? asa result of which the 
quantity y appears at the object's output, requires some 
time interval of constant duration, 

The presence of the noise z,, at the input gives rise 
to a random displacement of the minimum point y in the 
SPACe Uy, Ug,... » Up, Which imitates a slow “drift” of 
the characteristics which occurs in actual objects. It is 
possible to set Zoj = yyn, where n is discrete time or the 
number of cycles (n = 0, 1, 2,. . .) and the y; are random 
variables, The noise z at the object’s output (for exam- 
ple, in a measuring block) is a random noise of compara - 
tively high frequency. We shall suppose that the incre- 
ments of these quantities are independent random varia - 
bles with known distribution law ¥(v), zero mean value 
and finite dispersion @*, Such. a form for the noise was 
introduced for these systems in [1]. 

The Markov Process in a Gradient System 





of Automatic Optimization 





Let the dependence of the output y* of portion Fy, of 
the object on the input quantities xf be described by the 
formula ° 
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> Dik Li Le-+ 2 Cik 2, +d’. (3) 
i, k= k=1 
By applying a linear transformation to each variable 
separately, we can take this formula to the form 


m 
y =f (%1,.--52m) = > biz Li Lr, (4) 
i, k=1 
where by, = by. We shall assume that quadratic form (4) 


is positive definite, Then y has a unique minimum at 
the origin of coordinates, 

During the ith pair of intervals, whose _ num - 
bers within the cycle are j and j + 1, where j = 2i-1, 
measurements are made of the quantity y, We He 
the values of the qmaniigs x; in the jth interval by x), 
and the value of y by y' In the first interval of the ith 
pair, there is the quantity xf) + O,; at the ith input, 
while the second interval finds the quantity _ ~ Oo; 
at this imput. 

After these two intervals, we determine the differ- 
ence 


(3) x\?)” (3) 


Aiyn = f (Zi n,-->> + Goi, --. » Im, n}) — 
(j+1 j (j+1) 
1 (at ,, eres fit — Aoi,.- . oe ). 


(5) 
The subscripts n denote that the measurements were 
made in the nth cycle, The corresponding quantity at 
the output of controlling portion C (Fig. 1) is defined by 
the expression 


Aw, = Ain + Az, = Aiyn + Y. (6) 
The gradient method consists of this: that after a 
series of trials, as a result of which one determines the 
direction of the gradient of the quantity y — the vector 
with components @ y/0 x; & A,y;, the operational stepis 
taken in this direction. Obviously, the increments of the 


quantities See he where n is discrete time, must be 


proportional to the quantities AW found in the nth cycle; 


Yi, npr — Ui,n = — % (Ajwn) = — a;(Aiyn + 2). (1) 
By taking into account (2) and (5), and the formula 
= 7,0, we find that 


(2i—1) {tit 
Sd ae pees + Ai, 
(8) 
(2i—1 2i (2 2 
+92 Zm.n) — jm”. . tf — ag, . + wb r vi] 
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Here, we have denoted x) simply by x, ,. Since 
® ’ 














(j) (1) j—i i— 4 
x = 2 "—ieie Te a he 
; 9; 
0 a ee 
Li eee ee a, (9) 
then 
( fa 3 i—{ 
Ti, n+1 =%n- Yi — &; Kath — + 22 + Tt —, 
i—1i, {= 
t+ ti + Api, - »<m Ym —})— 
(10) 
“2i—1 2i—1 
—f(m+%1 2m pores Dit Ti Tl Sob. 


2i — 1° 
-»>Im-> Ym =) + vs| . 
For brevity, we denote Xin = &; and x; = X43 we 
introduce the vector x in the m-dimensional space ox 
the vector € with coordinates €; inthe space M = 
and the vector y with coordinates y; and the vecnint 
again in space Q,. Then (10)can be rewritten in the form 
i— 1 
% = Et yi— 04 (7 (F + agi Si +7 —") = 
(11) 
— (Fai Goi re — ~) + v;). 


We now take the form of function f(x) [cf., (4)] into 
account, By substituting (4) in (10) or (11), we obtainaf- 
ter some transformation, 


m 
m= Y— arity) Ayes + Bi, (12) 
j=1 
where 
Ay = = bi; bea xs bx ke + aida (= — 4a, |, 
k=1 
k+i 
Ye wae te 
i— 
B= =, % oo \ Ose 5K + 
j. k=1 
ju koi 


Pig (= —_— Zou) bii j7 - bi Xi (= — 2a)! ° (13) 


m 
» } 

Here, 54; = 1 for i = j and 6;; = 0 fori # j; ie, 
85 is a Kronecker delta. 

The sequence Xn (n= 0,1,2,...), oF the vector 

, is a random Markov process (cf, [2,3]) since the prob- 

ability distribution for x, , , depends on the value of x,. 

Initially, let the y,; be fixed quantities, In view of 
the independence of the y; and the ¥; , we can compute 
the distribution W(y) for y at any subsequent stage, Fur- 
ther, let the v; now be all zero. We first investigate the 
convergence of the automatic search process when the 
random noises are lacking. Equation (12) can be rewrit- 
ten as: 


x = AE +B, (14) 





3 we 


es 


$E ois 
> form 


(11) 


in af- 


(12) 


>ctor 
f Xp 
ww of 
pute 


»writ- 


(14) 





where A is the matrix with elements A;;. and B is the vec- 
tor with coordinates yj; + B;. If the system converges to a 
steady state x then this latter state may be found from 

the equation 


ste’ 


Xero = AXste + B. (15) 


Let E be the unit matrix, and let the inverse (E-A)~! 
of the matrix E-A exist. Then 


Xe = (EZ — A)"B. (16) 
We consider the differences 
Xnti — Xste = Ong, Xn—Xete = An. (17) 
If follows, from (14), (16),and (17) that 
Ana = Adg. (18) 


The necessary and sutticient conditions for the con- 
vergence of the search process, i.e., for the holding of the 
condition A, ~ 0 asn~ o, are the inequalities 


lre] < 4, (19) 


where the Ay (k = 1,.. .,m)are the roots of the charac- 
teristic equation | A-EA| = 0 of the matrix A. As is 
well known, these inequalities reduce to condtions of the 
Routh-Hurwitz type. We shall assume that these condi- 
tions are satisfied,* With this, the speeds yj; can not be 
arbitrary; the vector y must lie within a certain region 
ay of the m-dimensional space, 

We turn now tothe general equation, in which ran- 
dom noise is present, Let the transition probability from 
the point € ={X,n, Xyn,--.,. Xm,nJ} in the interval 
4S Xi.n+1.= xit dy (i =1,...,m) equal 9 (E, x). 
dxy. « - Xm = QE, x) dQy, where dQy = dx, . 

Then, Pp +4(x), the probability density that x, 4, 
appear in the interval x; =< Xi n+1= i+ dx; (i=1,...,m), 
is defined by the formula 


» ~ Oe 


. ( ? (§, x) Pn (é) dé,. o- din= 


—So 


? (§, X) Pn (§) dQ:, (20) 


where dQe = d&4,...d&_.. If this linear integral trans- 
formation has a fixed point to which the sequence Pn(x) 
tends, then by denoting 


P(x) = lim pp (x), 


we may find this function from the linear integral equa - 
tion 


(21) 

P (x) = \ (§, x) P (§) dQ. (22) 
Qe 

We find ¥(&, x), assuming now that the y; are fixed, 


from (12) we get 


m 
y= = [lv H+ >} Aj &; + Bi). (23) 


j=1 


The conditional probability that xj, n+ 4 is found in 
the interval xj S xj ny ys %j + dx if xj, = & (j= 1, 
«++,» m) equals ¥(v;) dx lo, where the v, are deter- 
mined from (23), Since, for different x;, the vj are inde: 
pendent, we have that 


m 


2x) =[] 2p I] tofgln—ai+ 





i=1 t i=l 


™m 
+ YAvh+Bi ft. (24) 
j=1 
Determination of the Indicators of the 


Automatic Search Process 








We now find the mathematical expectation of the 
random variable y,, by which we understand the arith- 
metic mean over the nth cycle; 


eo i—4 
Miyal=\ 3} Ue +ayby + 1st) 4 

x i=1 : 
, . 2i—1\ 
+4 (x— ay fo + ¥ =) | Pa (R) dQ. 

We consider separately the expression in the square 
brackets, By substituting (2) for f(x) we obtain, after some 
transformations, 


>) (7 (& + ai Soi +7 thf (XK — a4 Soi +7 “+ \ = 


i=1 


(25) 








$y: > bin 252 -- (1— 5) > bj, Aajx+C, (26) 





j. k=1 j, k=1 
where . 
=~ (m+i1)(4m—7)+5 ™ 
fn bse tite + 
C= a 2, ik Tsk 


1 m { m 
+ = >} bit %*?— mt > bij Mot Ti 
= 


(27) 
= i, j= 
Therefore, 
co co m 
M [Yn) = eee bi», Zs = Se. 0 
J J {3 be a0 
oe 
+i) >) bie 25.0% + Ch + 
j. k= 1 


Pn(2y, ni Le, ni --- Lm, nd, ndXe, n- . -d2m,n{28) 


or ‘“ e 


M |ynj = > us,4) Lj XL Pn (x) dQ, + 


j, k=1 x 


y 4 ee a 
+ (1-25, Y bet) typ dQe + 6, 


j. kK=1 


(29) 


* The stability conditions for continuous gradient systems, 
derived in [6], have well-known similarities to those given 
above for y; = 0. 
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since 


\ Pn (x) dQ, = 


x 


We introduce the notation 


Tix, a= \ ZLjLe Pn (x) dQ.., 


a, 
Jin = 25 pn (x) dQ. (30) 
a, 
Then, 
M [yn] = > Oj Tjx, n + 
j. k= 
4 m 
+(1— =) p bin Xk Jj, n+ C. (31) 


j, k=1 
Equations (30) for the integrals can be found from 
(20), We multiply both members by xj;x;,dQ, and integrate 
over Q,. Then, Ij, 143 appears in the left member of 
the equation, and the expression takes the form; 


Tin, np = \ 2X \ > (E, x) Pn (&) dQ; |dQ = 


" 
=§ 252 l\ Il ¥%) >) dQ, ]dQ,. (82) 
x — i=) “s 

We now reverse the order of integration, and change 
from the variables x;, Xg,... , Xp, to the new integrand 
variables v;, Vg,...,V,,- Since dx; = —ajdv;, and the 
change in sign of the differential compensates for the 
changes in the signs of the limits of integration, we then 
find, by letting dQy = dv}... . dv,p, that 


Isr, a] pr Ot) 36000 (y— 252; + 


v i=] 


+ > Aju. + Bi) (te — om OK + 
v=1 





m 
+ 3 Anke + Ba) Os) dO — | pn (8) Lx dQz.(33) 
me 
We now consider the interior integral L; ik included 
within the square brackets. By expanding the parentheses 
in it, we obtain 


+20 +00 om 
Ly = \ vee \ (I ¢ (r)| [Oey 5 Y_ V5 + V5 (— 25 Te — 
—oo —o i=] 
m 
— a; 5} Ay &y — &; By) + 
voi 
+, (—Oxnxj— a>) Ay & — Ar B)) + (34) 
v=] 
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+ (45 Te + ye By + 15 Be + Be Bi) + 





+(15 + TH) >) (Amy + Aj) & + (Bi + Be) > (Aw + 


v=] v=] 
m 
+ Ajy ) Sv + 7 Axy A ye - a dv), . . dm. 
B, v= 
Now, we take into account that 
+00 
| didn = 1, 
mo (35) 


+00 
\ vid (vi) dv; = 0, 
—oo 


and denote the dispersion of the noise at the output by 
_ 
o? = \ v8 } (x) dyj. 


—0o 


(36) 


By taking (35) and (36) into account, we may rewrite 
(34) in the form 
Ly, = 55, af o? + (15 %q + T7By + 1B; + B,B;) + 
m m 
+ Dl; + BA t+ y+ BYAle+ DY - 


v=] »,v=1 


Ai Ai, baby 
(37) 





By substituting (37) in (33), we obtain 
Ti, npr = Oju Me F + (15%y + 158, + 1B; + BjBy) + 
m m 
+ Dery + BD) Any + te + BAN. nt DY Ape Ain tons 
v=1 p, v= 
(38) 
It is possible to write m(m + 1)/ 2 such equations for 
determining tik n+1 fomIyy y- The m quantities 
Jy UY =leeees m) also appear in these equations, 
These latter oneitian can be found from (20) in an anal- 
ogous way. Indeed, by multiplying this expression by 
xj,n+14Q, and integrating over Q,, we arrive at the for- 
mula 


Inti =) 2 Ms (&, ¥) Pn (@) dQ, |dQ, = 
Q, 
¥(%) 
=| 


“| = a (§) dQ, dQ. = a 


x i=] 


(39) 





. eS of dy 
=} pat){\ 25 Jf eae] aQ; = \ pn (Ll 
Q- Qy i=1 2: 
With this, 
: m $(v;) a 
bg = ) 2 ll -. dQ. = mm I] ¢@ (45 — 417 
x i=1 . —co —0o i=! 
(40) 
+ ») Aye + Bj) dv,...dvum=1; +S} Ajviy + Bi. 


v=] v=] 








35) 





(39) 
(&) L; dd 


— a; Uj T 
(40) 


+ Bj. 





By substituting (40) in (39), we find that 


Jinsr =) pa) (45+ Aj’ + 


QE v=l1 


™m 
+ Bi) dQ, =(7; + Bi) + DS Awd, n- 
v=] 

Using (41), we can determine J aes from. the jy n 
and then substitute the result in (383 The presence of the 
noises vj is completely unreflected in (41), For vj = 0, 
the distribution densities reduce to unit pulse functions, 
but (41) remain unchanged by this, Therefore, the nec- 
essary and sufficient conditions for convergence for the 
sequence Jy yn asn->@ remain the same inequalities (19), 
The sequence of second moments lik, n of the distribution 
also converges for v; # 0. 

In fact, the presence of the noise w in (38) is reflected 
by the additive terms 6.) or o*, For 0? = 0 (absence of 
noises v,), conditions (19) are necessary and sufficient for 


(41) 


convergence of the sequence ly, ,. It is obvious that 


the presence of the additional terms $4) af a? only changes 
the steady~state value, but has no effect on the fact of 
convergence, since the transformation matrix of the lik, n 
[cf. (38) ] depends only on the Ajxe Therefore, inequalities 
(19), when random noises Vv; are present, also guarantee 
convergence of the first and second moments and, conse - 
quently, convergence of the sequence M[yy] [cf. (31) ]. 
We set 
LT jx = lim Tx, ni J; = lim J;, ns 
n-—-co n--co 


y = lim M [yp). (42) 


It follows, from (31), (38), and (41), that the quan- 
tities Tip, J i’ and 7 are defined by (43), (44), and (45): 


Tig = Sig ORS + (1 5%y + GBR + eB; + BB) + 
7 (rj + By) Ag, + (ty + By) Ay) J, + 
v=1 


m 


4. be A,, A; ‘ » G&. &=1,...5m). (43) 
p, v=1 
™m 
J,=(4 +B) + SY Ap dv(i=1,..-.m), (44) 
v=1 
m m 
] , . 
r= Sete +d) Sen yee. a 


}, k=1 j. k=] 


Thus, the problem of determining the mathematical 
expectation 7 of the quantity y in the steady-state mode 
is solved, A consequence of (43), (44), and (45) is that 
the error is defined only by the dispersion 9? of the noise, 
and does not depend on the form of the function ¥ (v), This 


property is a consequence of the quadratic characteristic 
of the object. The mathematical expectation of the de- 
viations from the minimum (in the given case, Ymin = 9 
in the steady state may be considered as a certain meas- 
ure of the error in this mode, As for the duration of the 
transient response, its mathematical expectation may also 
serve as its measure, This criterion was introduced in[4}, 

However, another measure turns out to be simpler — 
not the mathematical expectation of the duration of the 
transient response, but the duration of the transient re- 
sponse for the mathematical expectation, which is not the 
same thing, Consider the dependence of M[y,] on n [cf. 
(31)]. This is a transient response in some equivalent dis- 
crete system which may be studied by the ordinary meth- 
ods, for example by frequency methods, integral methods, 
etc, In the simple cases (cf,, below), one can obtain an 
expression for the duration of the transient response in such 
equivalent systems, 

If the vector y is not fixed, but has some probability 


distribution W(y), we may then, by denoting da, = dy; 
- +. dy,,, average the values of n over a, in accord- 
ance with the formula 

Nav = \ 9W (7) dQ,. (46) 


Q, 


It may be more advantageous for computational pur- 
poses to use (45), where the maximum expected valuesof 
the y; are substituted, 


Example of a Simplest System (m = 1) 





We illustrate the use of the formulas obtained by 
the example of a simple system with one input variable. 
We set 

A=, Or =O, =z, FL =F, t= 
= 7, bi; = b, B, = B, Ay =A. 

In this case, (19) for search stability takes the form 
| A] < 1 or, in accordance with (13), 


| A | = el — Saad + aybye< 1. (47) 
It is clear from this that the speed y must not go be- 
yond certain limiting values, i.e., y,;= y= y 2, where 
(Ys ¥2) ® 0. 
Equation (45) takes the form 


by ~2 
m= bly + Int b(a ay Se 2 TY (48) 
The integrals I,, and J, are determined from (43) and 


(44), which are rewritten in the form 








Thy = 262 + (y+ By)? +92A (y+ B) + A®%Jy, (49) 

Here, in accordance with the second of Eqs, (13), 
{ ~ “ 

B=— aby 5 (2a — <) = al; (+ —a) : (51) 

From these expressions, we find that 

1+B a%c? + (7 + BY (y+ BY 
A={oa ue “joa 4g a 

(52) 
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The transient response is determined by (31), (38), 
and (41), which take the following form: 


1% 7), & 
My = OL, + (a2 — ' ths T+ 


Jing = O70? + (7 + BYP +24 (7 + B) Fin t Aly» ; 


Jintp=Tt+Bt+ad, x (58) 
We set 
b+ bo 
emily. tas) —. k = adob, e=—, 
(54) 


Then, after substitution of the values of I,, and J; 
from (49) and (50), (48) assumes the form 


e= (4, ~ +) + () (ay + 
ky 
; (4) 1{—(1—4k + By 
We rewrite stability condition (47) in the form 





(55) 


i) 
1— 4k >” 


5 (56) 
Ecan j =1—4k+k7, TO K(4—z) = 


0 


=i—j=g. 


) 
If j=1—4k+k — pel-iza 


In a stable = Sl “1<j<i ar 0<q< 2, The 
boundary of stability is reached for 


zz: ——_ 
4-SA, (57) 
The other boundary is obtained from the condition 
that k > 0, i.e., 


i) 


4—G->0. (58) 


Thus, the region of stability has the form shown in 
Fig. 2. 

As an example, let us set 5 = 0,2, For this value, 
and by giving A», we find from(55) the minimal value 
€min 43 4 function of k, Let it be attained for k = Kopt 
(Fig. 3). If we trace the curve €pin = Emin (Ao (Fig. 
4a) and the curve Kopt = KoprlAo) (Fig. 4b), ~ can find 
the minimum possible value of the steady-state error 
€min 2nd the values of the parameters A, and k which 

min 


correspond to it. It is clear from Fig. 4 that €min ~ 
m 


0.3, Ay = 0.4, and k = 0,15, 

The optimal parameter values depend on the fact 
that for very large steps the error increases due to the 
significant deviations from the minimum; however, for 
very small steps the system can not compensate the effect 
of the "drift" y of the characteristics, thanks to which the 
error again increases, One easily convinces oneself that 
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for y = 0, the optimal values of and a, are infinites- 
imally small and, in the limit, the error ¢€ tends to zero, 
We now find the duration of the transient response 


where, for simplicity, we limit ourselves to the case when 
y = 0. Then from (51) and (53), we find that B = 0 and 








Ny = OT, an + b031 yng = O20? + AU), QA = 1— 4k. 
(59) 
The steady-state values turn out to be 
a2c2 ba2c? 
In=g Go)’ ™=BRAG— mit ote (60) 


Since 


Ly n= ate? + AY 


iS i. n= = O20? + A® [a%o? + AP 


++ = a2o? [1 + AP +---4 A? MD) 4 4] 


11, nal =" 


16°" 





—. 
= ato? - —5 +A™, = In (1—A™) + A™ yy 
(61) 
the difference 
6, = —21=5(1,) g— Ju) = bA™ (1, gp — 111). (62) 


It may be that at the initial moment I,; 9 < I, if 
the initial value of x is, for example, at the origin of co- 
ordinates, However, if we assume that the system is sig- 
nificantly distant from the steady-state position at the 
initial moment of time, we can assume that I;;,9-Iy 7 9 
Then, the quantity 6, tends to zero with increasing n, 
while remaining positive. We now find the discrete time 
Ng after whose expiration the following relationship is vali¢: 


6 


n 


zy; S1, (63) 


where the constant . << 1 is some admissable threashold 
It follows from (60), (62), and (63), that 


Ink (Ju. + af) —InU/,, 
InA 


The quantity ny serves as a measure of the duration 
of the transient response in the system. 


o—fu) 


(64) 
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Fig. 3. 
Example of a System with Two Variables 





(m = 2) 


Let the object’s equation have the form: 


y= x? + 2Ba x2 + x. (65) 
For |8| < 1, the unique minimum of the function y 
is found at the origin of coordinates, This case corre- 
sponds to m = 2, by; = bys = 1, and by = B in the general 
formulas, We further set y; = 0, &y = Gg = &, gy = Ag= 
= Xp, HyOXo; = Ayte, = k. Then, from (13), we find that 
Ay = Ag =1i—4k, Ajy=Ayn=—4Bk. (66) 
Equation (45) takes the form: 


7 = Ti + 2Blia + Ia2 + 8, (67) 


and ly, I, , and }, are found from (43), which can be 
written in the following forms: 


11,8k (4 _— 2k) + 8BA (4 _ 4k) Ih. _ 1687k?/ 99 = a%c?, 


I1)4Bk (4 — 4k) + 8k (1 — 2k — 28%) Jip + Ieg4Bk (1 — 4k) = 0 


— 168%k2J,) + 88k (1 — 4k) Jag + 1998k (1 — 2k) = a220?. 
(68) 
In these equations, the quantities I,, and Igg enter 
in identical fashion, as was to be expected, We set 


8Bk(i1—4k)=D, 8k(1—2k—2p%*)=£. (69) 
Then, solving (68), we find that 
aot a*c*D 
In =/2=jFa—pi la = fie (70) 
and, as follows from (67), 


2a*s? (E — 8D 
oF +03 (71) 
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We now introduce the relative error from the random 
noise (for this, we assume that aa = k) 


Ss 204, 


%—@5 2k*(E — BD) 
(ew. wil 
a) 





= 0(k, B). (72) 


By substituting the values of E and D from (69), we find 
that 


k (1—2k 
6 (k, p) = Ey 


1— B? 


47 — Dk) — tp (1 — 2k) + 4p? — Bp? (1 — 4h)? « (73) 





We find the stability conditions from the character- 
istic equation | A-EA]| = 0, or 


Au—h Ay aide 
An An—d|~ ve 


By substituting the values of Aij here, we arrive at 
the equation 


(1 — 4k — A)*? — 168%? = 0), (75) 
whose roots are 
Avg = (1 — 4k) + 48k. (76) 


The search is stable if | 4,,| < 1, It follows from 
this that k must be found within the limits 


1 
°<*<204TBD mm) 


For a given k, the quantity 8 may vary within the 
limits 


IBi<ap—1. (78) 
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Fig. 5. 


Figure 5 shows the function @( 8,k). It is clear 
from the graphs that the error from the random noise, 
generally speaking, increases with increasing |B]. Indeed, 
for B = 0 the lines y = const on the (x, xg) plane are con- 
centric circles, With increasing | 6 |, they become el- 
lipses whose axes are turned by 45° with respect to the 
coordinate axes, The greater| 8| is, the more 
“elongated® these ellipses become, which worsens the 
conditions for minimum search and increases the steady- 
state error from the random noise. 






For «->0 and a,->0, the total steady-state errorcan 
be made arbitrarily small. Certainly, if y, =k 0, there 
exist optimal values of a and a» which are greater than 
zero, The essential dependence of the error on &, anda 
(or k) (Fig. 5) shows that automatic tuning of a» and a, 
providing a minimum of the error, is advantageous in 
automatic optimizers [5). 

The author wishes to express his gratitude to B. M, 
Levitan for his counsel in carrying out this work, 
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The paper deals with a computer for the computation of estimates of certain statistical characteristics 
of stationary random functions and of linear objects’ characteristics, 


In the experimental investigation of the dynamic char- 
acteristics of controlled industrial objects by the method of 
determining an object's reaction to standard stimuli, the 
presence of “noise” limits, to a significant extent, the ac- 
curacy of the estimates of the high-frequency characteris- 
tics. 

The investigated object itself is frequently a source of 
“noise.” In a number of cases (for example, in the chem- 
ical industry), the “noise” level is so high, and the require- 
ments on the technological regime so rigid, that the inves- 
tigation of the dynamics of such an object by means of 
sufficiently large standard stimuli becomes impossible. 

Today, one of the methods of investigating the dy- 
namics of industrial objects under their actual working con- 
ditions is a statistical method, based on the correlation 
analysis of realizations of the input and output quantities, 
with the condition that these quantities, as well as the in- 
vestigated object's dynamic characteristics, are stationary. 
The basic difficulty in analyzing industrial systems by sta- 
tistical methods is that these methods require a large ex- 
penditure of labor in carrying out the computational work, 
in particular, in computing correlation functions or spectral 
densities. 

With the aim of decreasing the volume of computation 
al work, special instruments for statistical analysis were de- 
veloped, both in the USSR and abroad. The instruments de- 
veloped can be divided into two groups on the basis of their 
method of recording results: 

a) recording on magnetic tape, 

b) recording on paper. 

There are a number of special features in the techni- 
cal characteristics of instruments of the first and the second 
type. Instruments that record on magnetic tape are superi- 
or to instruments that record on paper as regards speed of 
processing statistical data; but under commutation condi- 
tions, or in the presence of magnetic fields, they are sur- 
pasced by the latter in reliability and in simplicity of use. 
Moreover, instruments that record on magnetic tape are 
designed, as a rule, for the processing of comparatively 
high-frequency random processes. From this point of view, 
it might sometimes be more convenient, in investigating 
the dynamics of industrial objects, to use records obtained 
on ruled (graph) paper. 


Computers of such type, developed in the USSR and 
the Amsler company ,were constructed on mechanical prin- 
ciples. Multiplication and integration in these instruments 


. are implemented mechanically with elements which com- 


plicate preparation and increase the instrument's cost. 
Moreover, the existing instruments which record on paper 
are designed for narrowly circumscribed goals, i.e., for the 
computing of correlation functions, while the problem of 
investigating the dynamics of industrial objects requires a 
more complete statistical analysis. 

Below, we describe an instrument for the statistical in- 
vestigation of the dynamic characteristics of objects in the 
presence of noise, and for the analysis of infralow-fre- 

iency random processes. Considered also are the possibi- 
lities of using such an instrument for estimating statistical 
characteristics, weight-function parameters and for the de- 
termination of objects’ frequency characteristics by statis- 
tical methods. 


Instrument Capabilities and Their Uses 





The instrument permits an operation of the following 
form to be applied to the functions x(t) and y(t): 


Ryz (47) = x\ zy t—*+)d. () 


ows 


1. Computing Estimates of Correlation 





Functions and Dispersions 





Expression (1) is an estimate of the cross-correlation 
function of the random processes X(t) and Y(t). Thus, the 
instrument allows one to obtain estimates of cross-correla- 
tion functions. This is an unbiased estimate which con- 
verges in the mean to the true correlation function as Too. 

If one sets X(t) = Y(t), then (1) gives an estimate of 
the autocorrelation function of the process X(t) and, in 
particular, for r= 0 we obtain an estimate of the dispersion. 


2. Computing Coefficients of Fourier 





Series 


If it is necessary to compute the coefficients of the 
Fourier series for an arbitrary function x(t), then one should 
choose, in(1), y(t) = cos wt andr = 0. For this, a reali- 
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zation of x(t) in the given interval (0,T) is written on the 


gtaph paper. The realization is reproduced in the form of 


a voltage on potentiometer brushes, supplied from a sine- 
wave generator which produced voltages of the form 
sin 2 t k/T)(k=1,2,...). 


3. Computing Estimates of Spectral—Den- 





sity Functions 





If x(t) is a realization of the stationary random pro- 


cess X(t), then the value of the spectral density at the point 


«, is computed on the basis of (1) 


T 
2 44° (2) 
M {ah) = 7\ Rex (2) cos wx + de-+ 0(4) 
0 
9 
a, = 7 \ x(t) cos w, t dt. (3) 
0 
If the interval (0,T) is sufficiently large, then 
(4) 


4 
M [2,7] = aw Gas (wx) + 0 (F:) ’ 
where G,, (w,) is the spectral density of the random pro- 
cess X(t). 


The value of the dispersion Mfaf J is computed approx - 
imately, starting from a finite number of realizations of the 


process X(t). 


4. Determination of Frequency Characteris- 





tics by Statistical Methods 





Even with a significant noise level, statistical methods 
permit one to find an estimate of the frequency characteris- 


tics of an object at whose input a harmonic oscillation of 


known frequency and amplitude is applied. Such a method 


is based on the following relationships. 
Let there be three signals: 


/,(t) = Dcos ot, 

/2(t) = Bsin (wt + 9) + n(t) 
and 

/3(t) = Dsin ot, 
where f (t) is the signal at the object's input, f,(t) is the 
signal at the object's output which includes the reaction 


Bsin(wt+ @ )to stimulus f(t) and the “noise” n(t), with the 
condition that M[mt)] = 0. 


It is then easily verified that 


T 
Ray 2(0) = lim 7-\ fa (0) fa (dt = 7 BD cos 0,(8) 
0 


T 
Ry,» (0) = lim 7-( fa(t) fo (t)dt = + BDsin6.(6) 
0 


On the basis of (5) and (6) we get 





1—2,V Ri, (0) + Ry g8(0),0~are tan 2) 
? ~~ Ja 12 (0) + ts, 2 (0),@=are a Ry 2 (0)’ (7) 
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where A = B/D is the value of the amplitude characteris. 
tic at frequency w, and @ is the value of the phase charac. 
teristic at the same frequency. 

It is easily seen that Ry 9(0) and R, (0) are estimates 
of the correlation functions forr = 0. These estimates: 
are increasingly accurate, the larger the time of inte- 
gration T. 

It is essential to mention that the computation of the 
integrals must be carried out for the function f,(t) for 
t = t, where t, is the time after which the natural motion 
of the systern will, for all practical purposes, have ceased, 

The instrument described in this paper permits one 
to execute the enumerated mathematical operations on 
the estimates of frequency characteristics. A scheme 
for taking off frequency characteristics by the statistical 
method is shown in Fig. 1. 
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Fig. 1. 1 is a sine and cosine generator; 2 is 


the object; 3 is the computer. 


The signal foom the servosystem, which continuously 
reads the realization f,(t) of the random process under 
investigation, is applied to potentiometers which are 
supplied from the sine-cosine generator. The quantities 
thus obtained, f,(t) sinwt and f,(t)cos wt, are applied 
to integrating motors. The indications of the counters 
on the motors’ shafts correspond to the quantities Rg x0) 
and R,,(0). 


5. Estimating the Parameters of Linear 





Objects’ Weight Functions 





Sometimes, to simplify a problem, one takes a rough 
general form of an object's transfer function; it is then 
necessary to determine its coefficients. As previously, 
we assume that there is “noise” in the object,and, con- 
sequently, we can only determine estimates of the de- 
sired parameters on the basis of recordings of the signals 
at the input and output of the object being investigated. 

We shall consider here a method of estimating the 
moments of the object's weight function by means of 
the instrument described, for the case of objects with 
monotonic transient responses* (the reactions of such 
objects to unit disturbances are monotonic functions of 
time). 

Industrial objects with monotonic transient responses 
occur quite frequently. In order to describe, in the fimt 


*The method considered here can also be extended to 
the case of nonmonotonic transient responses. 
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approximation, the weight function of such an object, 
one need estimate only several parameters. 

Let kt) be the object's weight function. By the 
zeroth-order moment of the function kt) we mean the 


number 
+co 


a = \ k(t) dt. (8) 


—co 


Since o» = © (jw) |,.9, where (jw) is the object's 
frequency characteristic, a is then the static gain. 
By the nth-order moment we mean the number 
+00 
an = \ i*k(t)dt (n=0,1,...). (9) 
—co 
For objects with monotonic transient responses, the 
weight function is completely determined, under suffi- 
ciently general conditions, by giving its moments. We 
note that for such objects the weight function obviously 
possesses the properties of a probability-density function, 
and the frequency characteristic possesses the properties 
of a characteristic function. For a complete analogy, 
we should introduce the normalized weight function 


k (ty = = k(t) (10) 


and, correspondingly, 
“— 
Cn = \ t"k(t)dt (n=0,4,2,...). (1) 
—oo 


Then, the analogs of the mathematical expectation 
and the dispersion will be the quantities 
— as! dy 


a, = om (12) 


in B-(Q-e—ao, 


The quantities ag, G&, and Gp turn out to be very 
useful in the investigation of objects, and these quanti- 
ties are frequently sufficient for the description, in a 
first approximation, of an object's functioning. These 
quantities have the following physical meanings: cy is 
the object's static gain; a, is the average lag of the trans- 
ient response; Ggg Characteristizes the distortion of the 
input signal. 

If the object possesses a pure lag, and does not dis- 
tort the form of its input signals, its weight function then 
obviously has the form: 


k(t) =8(t—T), (14) 


where T is the lag and 6(t) is a Dirac delta function. 
Then, on the basis of (8), (11), and (13), we obtain 
%=1, Oy = a, = T, and Cig = 0. 


For an object with the weight function 


—(t—T) a2 
k(t) = nm T, (15) 





we have 


Oy = 1,41 = T+Tyandaieg = T?4+277 4273, (16) 


This weight function corresponds to an object con- 
sisting of a series connection of a pure delay I'nk and an 
aperiodic link. We note that, on the basis of (16), we 
can easily compute the time constant T, and the true 
lag T of the object as functions of the average lag a 
and Of Gigs . 

Now, let Kt) be the weight function of an object 
consisting of two series-connected links with weight 
functions k¢t) and kg t). 

Then, 

+00 
k(t) = | ky(t— +) ha (s)de, (17) 


—co 


and for the transfer functions 


® (s) = D, (s) D, (s), (18) 
where 
+00 
®; (s) = \ ek, (t)dt (i =1,2), 
we have 


® (s) = { [tase + Fp... ]k (at =a,— 


— as +3... +. (19) 


On the basis of (18) and (19) we have, by equating 
coefficients of like powers of s, 
% = Yo Bo, 
1 = YoP1 + YiBo, (20) 


Oe = Yoo + 24:81 + Toho, 
from which it follows that 


a= 7 _ Bi. (21) 
20 a 120 + Boos 


where Yn and By are the moments of the functions 


kt) and k(t), andy, and By are the moments of the 
normalized functions kt) and left) [Cf.(10) }. 


We now consider the expression 
T 


Ryx (t) = \ k (t) Ryy(t—*) de. (22) 
0 


We assume here that T in (22) is chosen sufficiently 


large. By analogy with (20), we can write 
Ay = Abo, 
A, = box, + by%, (23) 


As = bot, + 26,0, + bore, 













































+00 
a, = \ t” Ruyx (t) dt (24) 
and 
+00 co 
bn=\ t” Ruy (t) dt=2 \ t” Ryy(t) dt(n = 0, 2, 4,...). 
ae7 ? (25) 
Since Ryy (rT) is an even function, we have 
A, = boy, 
(26) 


a, = bots te boty. 


Thus, the problem of estimating the a, reduces to es- 
timating the moments apy and b,, of the corresponding cor- 
relation functions, and to the solution of the system of al- 
gebraic equations (26). 

The instrument considered below allows the moments 
a, and b,, to be estimated. We consider the estimate of 
the an: 

On the basis of (24), we can write, starting from the 
record of a realization on the interval (t-Ts, t), 

a t 
\ ™ dt \ z(u)y(u—) du, (27) 
—) 


t—T, 


| 


The inaccuracy of the estimate in (27) is related to 
two circumstances: first, the infinite limits in (24) are re- 
placed by finite ones in (27) and, secondly, instead of the 
cross-correlation function T), we used its estimate, 
obtained by an averaging operation on the finite interval 
(t-Tg,t). 

It is necessary to bear in mind that the moment of 
time t in (27), at which one can begin averaging to ob- 
tain the estimate of the correlation function, must be 
chosen subject to the condition that | = T, where T is 


the object's “memory” (duration of the transient response). 

















With this, t = 0 is taken as the moment when the system 
begins to operate. 
By reversing the order of integration in (27), we get 





t a 
On’ = 7 \ z(u) du | (t)"y(u — t)dt = 
t—T, —b 
t a 
= \ tz@)) @ry@—ade+ 
t—T, 0 
b 
+ y (u) (— 1)" x (u —*) dt] du. (28) 


An estimate in the form of (28) can be computed in 
the following way. The realization y(t), written on one 
strip of paper, is read off in the form of a voltage and 
applied to a model, where one has chosen a system with 
finite memory a (with weight function (rT )"}. The voltage 
at the model's output is applied to a potentiometer whose 
moving contact is moved in accordance with the law x(t), 
The voltage at the potentiometer's output is applied to 
an integrating motor. The number of rotation of the 
motor gives the first term in (28); the second term is com- 
puted analogously. 


The functional block diagram for the instrument 
which computes these estimates is shown in Fig. 2. For 
estimating the moments of the correlation function in , 
(28), it is necessary to set 


x(t) = y(0), 
and then 
t a 
b,"= ra \ y (u) du \ ty (u — t) dt(n = 0, 2, 4,...). 
T, 


bu Se 
The conditions for choosing a, b, and Ty such that 
the estimates of the moments a, and b, will be sufficient: 
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ly accurate are easily obtained if the functions Ry yh t) and 
Ryx(T ) are known [1]. A more detailed discussion of the 
questions involved in the choice of a, b, and Ty may be 
found in [5]. 


Description of the Instrument 








The instrument's functional-block schematic is shown 
in Fig. 3, and its outside view is shown in Fig. 4. The in- 
strument is constructed physically of two blocks. In the 
first block are mounted the servosystem (reading), the 
tape-handling mechanism, and the multiplication poten- 
tiometer. 























Fig. 3. 1 4 a reading device (servo- 
system); 2 is a multiplier; 3 is a pow- 
er amplifier; 4 is an integrating motor, 
5 is a rotation counter; 6 is a stabilized 
power supply and 7 is a control panel. 


The operation of the servosystem is based on the prin- 
ciple of tracing by a photohead along the edge of a 
curve drawn on the graph paper. (By the edge of a curve 
we mean the boundary between the black and white re- 
gions.) The photohead consists of a photoimpedance, 


which is the sensitive element, and a small light bulb. 

The photoimpedance is placed on the front side, 
and the light bulb on the internal side, of the graph-paper 
strip. The photoimpedance is connected in one of the 
arms of the servosystem's bridge circuit. The bridge cir- 
cuit is supplied from a transformer. 

Depending on whether the illuminated area of the 
photoimpedance is increased or decreased, its impedance 
will decrease or increase, causing a signal of the proper 
phase to be applied to the input of the phase-sensitive 
amplifier, made up of twin triodes in a bridge circuit. 

The servosystem is executed on the basis of program- 
ing device RU5-01 (02). The block schematic of the 
servosystem'’s mock-up was changed, permitting the in- 
strument's dynamic characteristics to be increased, i.e., 
the servosystem's processing speed was increased from 12 
to 100 mm/sec. 


Depending on the phase of the input signal to a re- 
versible two-phase induction motor, whose control wind- 
ing is connected in the amplifier’s cathode circuit and 
which is turned in the proper direction, the carriage on 
which the photohead is mounted will move so as to fe- 
crease the imbalance in the bridge circuit. 


At the moment when the program's edge covers that 
portion of the light-sensitive layer of the photoimpedance 
for which the bridge circuit is balanced, the error signal 
will equal zero. 


Synchronous motors serve as the tape-handling 
mechanism of the instrument. A realization is carried 
separately on two reels of graph paper. Reading is im- 
plemented separately by the two reading servosystems. 
Such a kinematic scheme allows the lag Tr to be estab- 
lished comparatively simply and within broad limits. The 



























speeds of the tape-handling drives are synchronized. Mul- 
tiplication is affected by potentiometers, of which there 
are two pairs. The moving contacts of the paired poten- 
tiometers are rigidly connected to each other. 

The resistance ratio of the multiplication potentiom- 
eters was designed to be 1:20, for the purpose of avoiding 
the distortions of their characteristics caused by the load- 
ing of one potentiometer by the other. The potentiom- 
eters are made up of bridge circuits for the compensation 
of the voltages corresponding to the mathematical expec- 
tations of the random functions. 

In this block, the mechanism for reversing the tape- 
handling drives is also mounted, The second block contains: 


1) two ferroresonant voltage stabilizers with two 
half-period rectifiers and filters; 

2) two semiconductor stabilized—voltage sources 
based on Zener pedestal diodes, used for supplying the 
multiplication potentiometers; 

3) two semiconductor phase-sensitive power ampli- 
fiers for controlling the integrating motors, whose inputs 
are the multiplication potentiometers’ output circuits; 

4) two integrating motors with rotation counters; 

5) the control panel. 

The control panel is used for starting, reversing, 
and stopping the instrument. Instrument stopping is also 
effected by limit switches on the tape-handling mechan- 
isms. When the limit switches operate, the tape-handling 
mechanism and integrating motors, whose windings are 
shorted, are disconnected. The limit switches operate 
when an index from the pencil on the graph paper hits 
the transducer of the limit switch. Execution of the vari- 


ous computational task: on the instrument is implemented 
by various combinations of its elements. 


Technical Characteristics 





1 Speed of travel of the graph paper is 15 mm/sec. 
Maximum tracking speed is 10 mm/sec. 
Static-tracking accuracy is 0.5%. 

Supply voltage is 127/220 v. 

Supply-voltage frequency is 50 cps. 

Power required is 30 w. 

. The reading device's scale is graduated from 0 to 
100 %. 

8. The maximum length of graph paper is 16 m. 

9. The basic parameters of the ferroresonant voltage 
stabilizer are: input voltage 12 v, and output voltage 
22/25 v. 

10. The basic parameters of the semiconductor 
stabilized-voltage supply are: input voltage 22 v, and 
output voltage 8.7 v. With variations of the input vol- 
tage by +20%, the output voltage remains virtually un- 
changed. 


I OaPow 


11. The basic parameters of the phase-sensitive power . 











Experimental Characteristics 








Figure 5 shows the curve for the number n of rotation, 
of the integrating motor during 30 seconds as a function 
of the voltages u, and u, at the inputs of the amplifier 
aad motor respectively. 
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Fig. 5. The instrument's integral 
characteristic. 
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Fig. 6 

Figure 6 shows the results of applying operator(1) to 
the function sinw t + 6), it being known that the form of 
this curve is invariant (as T > o) with respect to trans- 
formation (1). Consequently, sinusoidal (or cosinusoidal) 
curves are conveniently used for verifying the instrument: 
accuracy. The small circles denote the ordinates of the 
curve computed by means of the instrument. As is clear, 
the deviation of the experimental function from the com: 
puted one does not exceed 4%. 

Figure 7 shows a correlation function computed with 
the help of the instrument. The realization of the corre 
ponding random process is shown in Fig. 8. 


Conclusions 





This work considered methods of analyzing low-fre- 
quency random processes, estimating the parameters of 
objects in the presence of noise, and also described an 
instrument for solving such problems. 

The instrument's computing device permits the fol- 
lowing quantities and functions to be computed auto- 
matically: 

1) estimates of mathematical expectations, auto- 


amplifier are: input voltage + 6.7 v, output voltage + 16 v, correlation functions, cross-correlation functions and dis 


current gain of 10, amplifier characteristic linear to with- 
in 2%. If the amplifier’s voltage supply varies from -10 
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to 20%, the voltage on the load remains virtually unaltered. 





persions of stationary random functions; 
2) Fourier coefficients of individual realizations: 
3) estimates of spectral-density functions; 
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4) estimates of systems’ frequency characteristics; 

5) estimates of weight-function moments. 

The instrument developed is very simple to use, and 
is very retlable, These two qualities make it particular- 
ly convenient for use in the analysis of industrial objects, 
and also nicely distinguish it from electronic computing 
devices with complex circuitry and low reliability. 

Among the superior features of this instrument is 
its capability of using realizations written on paper, 
since widely used industrial apparatus permits recording 
on paper tape. 

Among the disadvantages of the instrument is the 
large amount of time necessary for an analysis. 
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METHOD OF INVESTIGATING THE ACCURACY OF ESSENTIALLY 
NONLINEAR AUTOMATIC-CONTROL SYSTEMS BY MEANS OF 


EQUIVALENT-TRANSFER FUNCTIONS 
K. A. Pupkov 


(Moscow) 
(Translated from: Avtomatika i Telemekhanika, 
Original article submitted April 29, 1959 


Vol. 21, No. 2, pp. 191-200, February, 1960) 


A method is considered for approximating essentially nonlinear functions by equivalent frequency charac - 


teristics on the basis of a comparison of the spectra of the ra 
nonlinear element, 


One of the complicated problems in automatic-control 
theory is that of analyzing the dynamic accuracy and the 
forced motions of closed nonlinear automatic systems act- 
ed on by random disturbances. Today, the most effective 
method for investigating such systems is that of statistical 
linearization, developed by I. E. Kazakov [1], which starts 
from the representation of inertialess nonlinear elements by 
constant equivalent coefficients. This coefficient is most 
simply defined by one or another of two assumptions and, 
with this, two values are obtained for the coefficient 

Definition of the coefficient by the condition that the 
dispersions at the output of the true nonlinear element and 
of the equivalent linear one be equal leads to the expression 


(1) 


S 

Mo = 3, 
where o is the mean-square deviation of the random 
function at the input of tne nonlinear element and Oy is 
the mean-square deviation at its output. 

The definition of the coefficient starting from the re- 
quirement that the mean-square difference between the 
true and the approximating random functions be a minimum 
leads to an expression of the form 


pia ae neelee? (2) 
1 32 
x 
where Ry,(r ) is the cross-correlation function of the random 
processes at the nonlinear element's input and output. 
In both cases, the nonlinear element's transfer factor 
by average value is defined as the ratio of the mathemat- 


ical expectations at the output and input of this element: 
ky = (3) 


where m, and my are the mathematical expectations of 
the random processes at the nonlinear element's input and 
output respectively. 

Approximation by the coefficients in(1) and (2) gives 
good results with the condition that a narrow-band filter 
stands beyond the nonlinear element. These coefficients 
are conveniently used only in the analysis of the operation- 
al accuracy of simple automatic-control systems which do 
not contain complicated filters in their structure since, in 
this case, the representation of the true spectral content of 


m, 


mt 


? 
c 
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ndom process at the input and at the output of the 


the random process in the nonlinear system is lost. This 
latter circumstance also makes the synthesis of optimal 
nonline1r systems difficult. 

We turn now to the consideration of a method of ap- 
proximating nonlinear functions by equivalent frequency 
coefficients, or by equivalent “transfer” functions, which 
more completely mirror the true outlines of the random 
process's passage through the nonlinearity. 

This method allows one to compute, with as great an 
accuracy as desired, the dispersion of any coordinate of a 
nonlinear system, and opens up a path for the synthesis of 
optimal filters. 

An equivalent frequency-transfer factor for a nonlinear 
element is defined as the square root of the ratio of the 
spectral densities of the random process at the input and 
output: 


J % ur (@) 


Sin (@) ’ 
where Sj,(w) is the spectral density of the random process 
at the nonlinear element's input, and S,,,{) is that at 
its Output. 

Naturally, K(w, 0,, Mm ) is also a function of o, and 
m,, and depends on the type of nonlinear characteristic. 

For the determination of the frequency characteristic, 
it is necessary to know the spectral density of the process 
at the output. We assume that we know the two~dimension 
akdistribution law and the correlation function of the ran- 
dom process at the input of the nonlinear element, where 
the distribution law is a normal one. 

For a stationary random process, the correlation func- 
tion (the autocorrelation function) is defined by the ex- 


K (w, ex, mz) = )/ (4) 


pression 
R,()=\ \ Ly — M Lp — Mz)W,, (1, Xo, t) dy, Azo, 


where W(X}, X_,7) is the random process's two-dimension- 
al probability density. 
We now consider the nonlinear functional dependency 
y =f (2). (6) 


The autocorrelation function RIT) at the nonlinear 
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element's output will be the average value of the product 
xdy(t+ 7), which depends on the form of nonlinear func- 
tion (6). 


If the values of the input signal equal x(t) and x(t+ rT), 


then the values of the output signal will be f[x(t) ] and 
f{x(t+ T)]. 

Consequently, the probability of the product f[x{t)] 
f(x(t+ T)] will determine the two -dimensional prob- 
ability density Wg (xX, X, T). 

The mean value 


M {(f[z())flz(t +9) m 


is the autocorrelation function of the random process at 
the nonlinear element's output. 


Expression (7) can be written in the alternative form 
Ry() =) \ 4) f (es) We (21, 22, 9) dey day, © 


where 
r= 2(t), = 2z(t+7). 


In (8), the random process is taken to have zero mean 
value. 

For the computation of the integral in (8), we approx - 
imate the nonlinear function of (6) by rectilinear segments 
and,for Wg (X4, %g, T), we use Cramer's expansion by 
terms of orthogonal Hermite polynomials [2]. 

The Cramer expansion has the form: 


2 2 
- exp | — = > Pay 2m |= 
1— p* 2 (1 — p*) 








- > H, (2) Hy (2) _ 4 


v! 





(9) 


v=0 


R(r) 


where ps RO) is the normalized autocorrelation func- 


tion, and H,,(x) (v =1,2,3,...) are orthogonal Hermite poly- 
nomials. 


For the two-dimensional normal probability-density 
function, we get 


ee ex peat + pre — 2pa re Seat zt + zs 
®Vi-p 2 (1 —p*) err 


2 27] © 
= exp|— 14/5 H, (%;) H, (22) ” (10) 
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x + 4) © Hi, (x1) H, (22) 


2 > 2nv! 


= il — 
L ong 


p”. (11) 


We now normalize the random variable x with respect 
to the mean square value o, i.e., we setE 1 =X, /O, and 
Es = Xq /Og. 

In the sequel, we can take oy = Og = o, since the ran- 
dom process is stationary. 


Sometimes, the integral in (8) can be written in the 


form 
4 





Ry (t) = inVi-p x 
fs (12) 
g + & — 2pEiEs . _ 
g | \ _ |— Te | 1 (6%) £ (2%) dé,dé,, 


By joining (11) and (12), we obtain the following ex- 
pression for the autocorrelation function: 


{ co 
Ry (t) = Devi >) p” x 


oat gee 
x) | ok) f (ek) e AG) Hy Gs) di dp. (219) 


The double integral in (13) decomposes into the pro- 
duct of two simple ones. Then (13) may be given in the 
form: 


Ry (t) = >} ao’, (14) 
v=0 
where 
Te | -= 
2 


The orthogonal Hermite polynomials are defined by 
the relationship 


(q)Je *=(—1tHe * aa 


and equal, for various values ofv, respectively: . 


Ho) =1, Hi() =§, 426) =F —1, 
H,() = ®—3§, H,() = &—66+3, (1 
H, (§) = & — 108 — 158. 

Thus, the autocorrelation function at the nonlinear 
element's output can be determined by computing the 
series in(14). This series is absolutely convergent, since 
p(t) is a bounded function, and the coefficients a,, de- 
crease at the rate 1/v'. 

If the nonlinear device's characteristic is an anti- 


symmetric function, then (14) contains only odd terms 
of the expansion 


Ry (t) = aip + agp* + afp> + --- (18) 


If this function is odd, but not symmetric, then all 
the terms of the series are retained. With this, the co- 
efficient a equals the square of the mean valueat the 
output which results from the lack of symmetry. 

In this case, the autocorrelation function for the 
centered random process will be defined by 

co 
R(t) = D av. (19) 
v=1 

From a known autocorrelation function, we obtain 
the spectral density by a Fourier transformation. 
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We may now turn to the definition of the equivalent 
frequency transfer factor. For the determination of this 
' coefficient, it turns out that it suffices to know the ay 
and the spectral densities of the powers of the autocorrela- 
tion function at the nonlinear element's input. 

Then, by the definition of (4), K( w, ox, m,) will 
have the form: 


K (@, Sx Mz) 


Dab, (20 
v= 1 


where 
S” (w) 
Sin (o) ’ 
For example, for correlation functions and spectral 
densities of the form 


by = SY (w) +p" (2). 


2B 


Px (t) = e— Fit, Sx (@) = +p 


(21) 


the coefficients b, equal: 
by 


by 


1, 

w* + B* 
= 2a ‘oe + 4B? 

w* + BF 
=3 sy 
oad w* + B? 
b=455 168?’ 
itt w? + 8? 


bs 


(22) 


o o..4, 612 ae © & 


In this case, the transfer coefficient equals 
_o? + BF 


=5V a+ 20h + 2a3 2a ‘o? + 4B 


In accordance with (14), the coefficients a’ are 
proportional to the dispersions, and may be computed 
for the standard essential nonlinearities. 

- In the given case, the transfer factor by mean Value 
is defined by the expression 


K (@, ox, Mx) = 





+ 305 SR +--+ (23) 


km, — 4% 
m, a 


ko (mz, ox) = (24) 
where k is the slope of the nonlinear characteristic at the 
working point. 

We now consider (20) in more detail. If we assume 
that the spectral densities at the nonlinear &lement's in- 
put and output have the same form, then the coefficients 





by = 1. Asa result, we get the transfer coefficient ° 
K (ex, mz) = — Sata kM, (25) 
Sx 
v=1 


which corresponds to equivalent transfer factor (1), com- 
puted by the first approximation method, since 
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If we take only the first term of the series, we obtain 
the transfer factor 


K (sx, mz) = — =k), (26) 
which corresponds to the equivalent coefficient in (2), 
computed by the second approximation method. 

If we neglect all the coefficients a,, except ay, this 
means physically that we are not taking into account the 
components of the random process due to the nonlinear 
transformation. 

Thus, (25) and (26) serve as the mathematical bases 
for the two methods of statistical linearization of non- 
linear functions. 

We now turn to the consideration of the equivalent 
"transfer" function for a nonlinear element. Taking into 
consideration that the equivalent frequency-~transfer factor 
is taken as the ratio of spectral densities, which are even 
functions of frequency, we present it in the form of a pro- 
duct of complex conjugate factors: 


K*(jw, ox, mz)K(jo, ox, mx)=K*(w, ox, mz). (2 


By choosing the zeroes and poles in the upper half of 
the complex plane, we obtain the following expression for 
the equivalent “transfer” function: 


(T,jo + 1) (T37@ + 1)... 








K (jo, ex, ms) =ke “Piette tt). * 
where the value of 
1 Tees cee 
ke = V4 hS2+ sd+- (29) 


is found from (23) for w = 0; Ty, Tz, Ts,... are time con- 
stants, and are functions of the coefficients a,,, and of the 
form of the spectral density at the input, for which T,> 
> Te. Ts > Tes etc. 

The coefficient in (29) has a value lying between the 
values of the coefficients in (25) and (26), and can give as 
accurate a value as desired for the computation of the 
dispersion in a nonlinear system. 

However, this fact holds in the case when, at the in- 
put of the nonlinear element, there acts a random process 
with a spectral density in the form of “white” noise, or 
with a frequency band which exceeds the pass band of the 
linear filter which follows the nonlinearity. In all other 
cases, it is desirable, in analyzing dynamic accuracy, to 
use the equivalent “transfer” function of (28). 

. Since (20), (28), and (29) have been obtained in the 
form of series, it suffices for practical computations to 
limit oneself to five, and sometimes even three, terms of 
the expansion. 

It is obvious that, both for the computation of equiva- 
lent transfer factors and for the determination of the equivé 
lent “transfer” function, it is necessary to know the coef- 
ficients a,,, defined by the integral in(15). These coef- 
ficients are functions of the dispersion of the input random 
process, and depend on the form of the nonlinear function. 

The computational formulas for the a, may be given 
for any nonlinear characteristic. 
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In (15) the nonlinear function f(o €) enters, complicating 
the computation of the integral, so that we approximate 
this function by rectilinear segments. Then, the integral 
of (15) may be easily computed by using tables of the in- 
complete gamma function and of the x? probability func - 
tion [3]. 
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We illustrate this methodology by an example. We 
take a nonlinearity of the limitation type, as shown in 
Fig. 1, and write the general expression for the a,» Which 


has the form: 
co 


ke \ -= 
% = Vinvi 2, tH, (€) e d= + 





k we 
+ Vinvi «(2 — %)H, (Be * dE + 
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_ko_ me tat 
+ Vinvi J cH, (§) e 


ajc 


— = 
y (—b—of) Hy (—)e * dé. 


—oo 


k 
+ V onvi (30) 


In the integrands of (30) there appear powers of the 
variable which depend on the order of the polynomials 
in (16). 


These integrals may be reduced to the tabulated ones: 





as 4 a 
Jee d= V2" (*=*)111— Pas, n+ 1) 
ry (31) 
or 
1 r -< 
sa | te = Ply? rn). (22) 


The function P(x?, n) is a function of two arguments 
for which tables have been computed (Cf. [2]). 
In the given case, the quantity x is defined as 


aan (33) 


6 , 
where 


d=h+m,. 





Let us consider the methodology of computing the a,, - 





For the nonlinear characteristic shown in Fig. 1, 


h—m, 


u1= 





(34) 
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Using (34), set up for the nonlinear characteristics 
shown in Fig. 2, we have calculated and constructed the 
graphs for the normalized coefficients 4ny which are 
given in the Appendix in Figs. 3-28. The parameter m 
on the graphs has the scale of the quantity x’. 

In computing the coefficients a,, from these graphs, 
it is necessary to take the following relationships into 


account: 


1) for nonlinearities of the limitation type (Fig. 2,a 


and b), 


ay = any L 


2) for dead zones with limitations (Fig. 2, c), 


L 
~=4y Te" 


h 


2=—, 


r 


3) for dead zones, cubic parabolas and antisym- 
metric quadratic parabolas (Fig. 2,d,e and f), 


v 


= an, N 


Example. In conclusion, we give the example of 
the determination of the equivalent "transfer" function 


for a nonlinearity of the cubic-parabola type of Fig. 2,f. 
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10, 3B? (3a?+ a?) ays 
‘eh 
35 Oe he. 9 
(35) Ts= 35 = 30° 
Then, 
0,0524s + 1 
K (6, og, m,)=75 DEE T. (36) 


The investigation of the accuracy and the forced 
motions of closed nonlinear automatic-control systems 
leads to the replacement of the nonlinear elements by 
equivalent transfer functions, using the method of suc- 
cessive approximations for graphicoanalytic methods, 
The transfer function considered here can be used for 
determining the conditions for the investigation of sta- 
bility of nonlinear systems. 
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Fig. 5. 
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PROCESSES OF FINITE DURATION IN PULSE-WIDTH SYSTEMS 
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Equations are derived for pulse-width systems which determine the form of the transient responses which 
terminate in a finite time. For second-order systems, methods of realizing such processes are given. It is sug- 
gested that the complicated nonlinearities necessary for the realization of such processes be replaced by linear 


Today, much attention is being given to those opti- 
mal processes which terminate in a finite and, moreover, 
a minimal interval of time. In particular, such types of 
processes may be realized in relay and in pulse systems. 

The duration of such a transient response is deter- 
mined by the order of the differential equation of the 
system's continuous portion [1, 2]. Finite-duration pro- 
cesses in pulse-amplitude systems with amplitude limita - 
tions were considered in [3, 4]. 

The present work is devoted to the investigation of 
finite-duration processes in pulse-width systems, where 
the pulse height is given, and the duration is limited by 
the repetition interval. 

In order that the transient response terminate at 
some point of discrete time, t = / T, it is necessary that, 
at this moment, the error quantity and its first (n-1) 
derivatives, inclusive, equal zero, where n is the order 
of the differential equation describing the system's con- 
tinuous portion. These conditions may be satisfied by 
the proper choice of the n durations of the pulses which 
act on the system up to the moment of time t= 1 T{z =n). 
Thus, with a limited external stimulus, the transient re- 
sponse can terminate after n cycles of the pulsed (sampled - 
data) element, just as in sampled-data systems with am- 
plitude modulation. We shall call such a process an n- 
cycle one. The problem of investigating a transient re- 
sponse of finite duration reduces to the determination of 
the form of this process and of its realization. In this 
work we derive the equations which define the indicators 
of such a process and, for a system described by a second- 
order equation, we present methods of realizing it. Fx- 
periments on a type MN-7 analog computer confirmed the 
results obtained, and showed that the complicated non- 
linearities necessary for realizing finite-duration process- 
es may be replaced by simple amplifiers. 

Conditions for Finite Duration 
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elements, Results of experiments on a model are given, confirming the possibility of such replacements, 


We shall consider a control system which is com- 
prised of a continuous portion and a pulse (sampled-data) 
element which carried out pulse-width modulation (Fig.1) 
Let the continuous portion, which contains an integrating 
link, be described by a differential equation of the form 


ayz’” (t) + adn_yz”” (t) + + - + + yz’ (t) = y(t), (D) 


where 2t) is the system's output quantity and y(t) is the 
pulse element's output quantity. The operation of this 
latter element is described by the expression 
Kpsign o (kT) 
y (t) | for KAT “<Ct<c kT 4+- TQ{([o(kT))}, (2) 
0 for AT + TQ[o(kT)) <t<kT +T, 
where o(t) is the controlling stimulus and 


1 
O (op = f= |o| <x, (3) 
{ for jc|>-x. 

It is clear from (2) and (3) that the pulse element is 
characterized by three parameters: the repetition inter- 
val T, the gain k, and the sensitivity « . 

At the output of the comparator we shall have the 


quantity 
a(t) = F(t) —2(0), (4) 


where F(t) is the external stimulus. 

In the ordinary system, the controlling signal o (t) 
equals this difference, i.e., o(t) =x (t). We assume the 
external stimulus to be a step function. 

In order that an n-cycle process occur in the system 
for a step-function stimulus of height 6, it is necessary 
that the input quantity to the pulse element o(t) at the 
discrete moments of time, 0, T, 2T, ..,(n-1)T assume 
definite values o(0), ofT), ...., o(mT-T) such that, at 
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time t = nT, the system's output quantity at) take the 
value 8 and its derivatives, up to and including the (n-1)- 
th, vanish. Thus, the input quantity y(t) of the continuous 
portion is a sequence of n rectangular pulses (Fig. 2) whose 
durations are determined from the conditions for a finite 
duration of the process: 


z(nT)=8, 2i(nT)=0 .,n—1), (5) 


“aan express the £) (mT) (j = 0,1,2,...,n-1) in 
terms z**’ (kT) (j = 0,1,2,..,n-1), where m and k are ar- 
bitrary numbers, with n = m>k 20. It is obvious that 
the quantities mo) mT) can be found as the reaction of 
the continuous portion to the sequence of rectangular 
pulses shown hatched on Fig. 2, with initial conditions 
2) (KT) (j = 0,1,2,....n-1). 

If we apply the Lapace transform [5] to (1) with the 
initial conditions given above, and solve the expression 
thus obtained with respect to the transform of the output 
quantity, i.e., with respect to Z(p), we find that 


(7 


at oY 








n— 
Z(p) = Si Alp) 2" (KT) + K (p) Yun (9), 
where ~ 
Ai (p) = -t# + 9 49P + +++ +0np" : 
ayp+---++a,p" 
K (p) = 7 _ www t 
and Y)._(p) is the transform of the hatched portion of the 


sequence of rectangular pulses (Fig. 2), 
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Yim (Pp) = hy a sign [¢ (i7)] e—? - 


Using the well-known rules [5] for finding the trans- 
forms of the derivatives, we get 
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(3) = sie 43) 
L {263 (t)} 2 (p) 23) (kT) + m 
= ; __ p—pTQLoliT)] 
+ ky piK(p) I sign[o(i7)Je-¥7—0 =P . 
ik 
where 
aes + 4;,.p +-+++a,p"™ 3 | 
{jit cy... er for J} Si, 
Bii(p) = . sf (7 
—p— a for j>i. 
ap t+ eo « +a, p™ 


If we go from (6) to the time-domain function for t = 
=(m-k)T, we obtain the relationship connecting the quan- 
tity 2t) and its derivatives at time t = mT with the quan- 
tity &t) and its derivatives at time t = kT. 

n—] 


29 (mT) = NL {By (p)}tan—wy12 (AT) 


1 
n—]) 


kp NS sign [2 (iT)] Ue? (mT — iT) — 
i=k 


-— (mT —iT —TQ[s(iT)))}, 


(8) 


where ne 


function of the system's continuous portion. 
This last expression is necessary for us in our investi- 
gation of second-order systems. Here, we shall make use 
of it in order to express analytically the conditions for 
finite duration of transient response (5) for a system of any 
order. By setting m = n and k = 0 in(8), and by taking in- 
to account that 23) (0) = 0(j = 0,1,2,..,n-1), we find that 


n—l 


kp by sign o (iT) {kh (nT — iT )- 


1i=0 





) is the jth derivative of the transient-esponse 


om kat —iT —TO [o (iT)]})} = B, 


m1 


S' sign ¢ (iT) (h® (nT — iT) — 


1—0 
— h (nT 
G = 4, 


9 


hy eee 


n—1). 


(9) 


— iT —TQ{s(iT)})} = 0 


In order for an n-cycle process to be realized, the con- 
trolling stimuli o(t) must assume, at times 0,T,2T,..., 
(n-1)T, the values defined by (9). These equations are 
transcendental so that, just as in the case of relay systems, 
it is impossible, in the general case, to solve them for the 
o (iT) (i = 0,1,2,..,n-1) or to state a general method for 
realizing such a process. We therefore limit ourselves in 
the sequel to the consideration of special cases, specifi- 
Cally, to second-order systems. 


Second-Order Systems 





Let the continuous portion be described by the dif- 


ferential equation 


1,2" (t) + 2 (t) = ky (2). 


By comparing this expression with (1), we note that 


T; 


n= 2, oe hae 


from whence 


h(t) = L" {— K (p)} a L} 


ay 


1 


k 


t 


= k[:—7, (1 —e =) 


h(t) = k(4 — Fh I 


t 


p* (Tip + 1) 


k 


k? K (P) = Srp FD’ 


| = 
(10) 


(11) 


If we take the inverse Laplace transform [5] of (7), 


we find that 
+ +p 
| {Boo (p)} =" 4 T 1 [‘], 
1 2 
ra a. 
“ - k 
L™* {Bo (p)} =L tor T, i =T, (1 
%?P a — 
L-* {By (p)} =L™ {0} = 0, 
T; 
7, 


L* {By (p)} =L*4 p 














(9) 










lif- 


















(11) 


(7), 


(12) 


:) (13) 


(14) 


(15) 


If we substitute (10) and (11) in (9) for n = 2, and 
carry out some elementary transformations, we can write 
the conditions for finite duration in the following form: 


F Q(T) 
sign o (T) (e" —1) + (16) 
: Q[a(0)] —F 
+ signs (OV(e™ —i}e " a@ 
and, by virtue of (16), 
z(0)=p= kk oT sign o (7) Q[o(7)] + 
(17) 


+ kk,T sign o (0) Q2[¢ (0)}. 


In accordance with (3), the parenthesized terms in 
(16) will be greater than zero and, in order that the left 
member of (16) equal zero, it is necessary that the quanti- 
ties o(T) and o(0) have opposite signs. With this taken 
into account, we solve (16) for Q[{o(T)] and substitute 
the expression thus obtained in(17). The result is 


z(0) = No{s(0)] = khyT sign « (0) {QI (0)] — 
a a Qf o(T)] i (18) 
—~in(e hii —e m4) 
Where the function Q[{o (0)] is defined by (3). The inverse 
relationship, i.e., 


¢ (0) = No" [x(0)], (19) 
is the relationship sought between the error x(t) and the 
input function o(t) of the pulse element at time t = 0. 

At time t = T, the relationship between x(T) and 
o(T) has a different form. To find it, we use (8) for k=1 
and m = n = 2, in which we substitute (10)-(15)for t = - 
(m-k)T = T. If we substitute the values thus obtained 
for z(2T) and z’(2T), expressed in terms of z(T) and 
z'(T), in conditions (5) for finite duration, and solve the 


equations obtained with respect to z(T) and z'(T), we 
get 


2(T) = 8 — kk, T signa (T) X 
(20) 
—1)\=8 —Nyfo(7), 


Fy 
é 


x {afo(ry) — 2 


r Q(T )) 


2’ (T) = — kkpsign ¢ (7) (c ; 
= — NJo(7)}. 


—1) . @amp 


By noting that, for 8 = const, x(T) = 6 -z(T) and 
x' (T) = -z'(T), we find from these expressions that x(T)= 
Ny[o (T)], x*(T) = Ng[o(T)], or 
(7) = Ny" [x(T7)I, (22) 
3 (T) = Nz" {2' (T)). (23) 
Relationships (22) and (23) define the relationship 
between the values of the controlling function and the 
error, or the velocity of the error, at time t= T. It fol- 
lows from (19), (22), and (23) that the output quantity 


of the controlling device can be given by the following 
two methods: 


- a 
sie Lets (0)} for t= * (24) 
Ny[x(T)] for t=7 
and 
"- wee [x (0)} fort = 0, 
~ (Wat ta’ (1) = — Netz (Ty fort = 7, 


In the intervals between the cycle times (the discrete 
time moments), o (t) can assume arbitrary values. 
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Fig. 3. 


In accordance with (24) or (25), the controlling de- 
vice is comprised of two nonlinear elements and a relay, 
which switches the error signal from one nonlinear ele- 
ment to the other (Fig. 3, and b). Realization (24) (Fig. 
3, a) is worse than realization (25) (Fig. 3, b) since, if 
the external stimulus varies, not by jumps but arbitrarily, 
then at time t = T the comparator’s output quantity in the 
first case will not correspond to a two-cycle process, and 
the system will not process in two cycles that stimulus 
which was measured fort = 0. In the second case (25) 
(Fig. 3,b), at cycle time t = T, control is effected by the 
velocity of the system's output quantity, and on this the 
variations of the external stimulus have no influence. 

Example. As a numerical example, we consider a 
system in which the continuous portion is described by the 
differential equation 

Z''+2Z' = y(t), 
i.e., Ty = 1 sec, k = 1 and the pulse element has the fol- 
lowing parameters: k. = 20v, x = 20v and T = 1 sec. 

To this system, we apply the second variant of the 
realization of a finite-duration process. By substituting 
the parameter values given above for the continuous 
portion and the pulse element in Expressions (18) and (21), 
and by taking (3) into account, we find that 


z (0) = No{c (0)] = 
| (0) | 





o (0) — .Osigns (0) In (0,368 ¢ so + 0,632) 
_ for |¢(0)| < 20, 
10,3 sign o (0) for |o(0)| > 20, 


— Z'(T) = N2[s(T)) = 

}o(7)| 

= 20 sign 6 (7) (e 20 _1) for |3(7)|< 20, 
34,4 sign o (7) for |o(7)|> 20. 
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Fig. 4. 


The graphs of the relationships x(0) = N, [0 (CO) and 
-Z'(T) = Ny [0(T)] are given on Fig. 4, a and b. For the 
realization of a two-cycle process with nonlinearities in- 
verse to these, it is necessary to use control by the scheme 
shown in Fig, 3,b in the system. 


It turned out, with empirical investigation on a model, 


that these nonlinearities can be replaced by amplifiers 
which roughly approximate them. It is clear with this 
that the saturation zone can be ignored, since the pulse 
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element itself, according to (3), possesses saturation for 
values of the input signal o(t) which correspond to the op- 
set of saturation in the nonlinearities Ng and Ny. Thus, in 
the given example, the nonlinearities No {x (t) ] and 

Ny? -Z' (1) ], shown in Fig. 3,b can be replaced by ampii- 
fiers with gains approximately equal to 2.5 and 0.6 respec. 
tively. 

Figure 5,a shows the oscillogram of the transient re- 
sponse of a finite-duration process in the system for a stim- 
ulus F(t) = 8 -1[(t] = 7v, from which it is clear that the 
transient response terminates after two cycles of the pulse 
element, while the ordinary process continues for approxi- 
mately six cycles. By using the standard linear-correction 
methods, for example, the introduction of derivative con- 
trol, one can significantly improve the transient response, 
but this improved response will last longer than two cycles 
(Fig. 5,b). 

In conclusion, we make several remarks about the 
behavior of such systems for large amplitudes of the ex- 
ternal stimuli. From the oscillograms shown in Fig. 6,a 
and b, it is clear that, for 8 = 15v, the duration of the 
finite-duration transient response (Fig. 6,a) is comparable 
with the duration of the response in the system for which 
the relationship between stimulus magnitude and error ve- 
locity was selected in the very best way (Fig. 6,b). Thus, 
the method presented above for the realization of a two- 
cycle process becomes of little effect when the stimuli 
are large. 

Equation (18) gives the relationship between the quan: 
tity which the system can process within two cycles, and 
the relative duration of the first corresponding pulse 9/40) 
which cannot be greater than unity. Thus, the maximum 
amplitude of the external stimulus which the system can 
process in two cycles is determined by (18) for Q[e(0)=1,l« 


T 
s 


= kk r[i—Ftin(2a—e ™)]. (26) 


v 
® max p 7 


In our numerical example, X) » yx = 10.3. For Kt)= 
B> Xgmax» the system processes Xx during the first two 
cycles and, during the nexttwocycles, the quantity B=% my 
If this latter quantity is less than Xpmax, the process 
terminates after four cycles, otherwise it continues for 
even longer. Figure 7,a shows the approximate form of 
such a transient response. As is clear from this figure, the 
response consists of a series of two-cycle processes. In 
this case, even the introduction of a simple combination 
of the quantity and its derivatives can result it an improv 
ed transient response. This,however, occurs only for 
B> Xomax: 

it follows from (26) that the magnitude of Xp max #8 
proportional to the repetition period T.If B> xq max, the 
to obtain an even shorter finite-duration process, it is 
necessary that the first pulse not terminate at time t =7, 
but continue even longer, i.e., the system must operate 
as if the repetition period were 2T, not T. The approxi- 
mate form of such a process is shown on Fig. 7,b. The 
dashed lines on Fig. 7,a and b show the variations of the 








for 





he on- 
lus, in 
d 
ampli- 
respec- 


nt re- 

' a stim- 
the 

e pulse 

pproxi- 
Tection 
je con- 
sponse, 

> cycles 


the 

e ex- 
3. 6,a 
the 
parable 
which 
ITOr ve- 
. Thus, 
a two- 
muli 


the quan- 
es, and 
se Qfa'o), 
iximum 
>m can 
5(0)=1, Le 


(26) 
or F(t) = 


first two 
B“% mar 


cess 
2s for 


>rm of 
gure, the 
ss. In 
yination 
1 improv- 
for 


max Is 
naxe thet, 
it is 

ne t=T, 
yperate 
approxi: 
. The 
s of the 






My (rit) 








\eéey_ —_ 





a 




















re ru ft 





azit}+z) 
ree Dee eat 





ae pa 











Fig. 6a. 















if, in them, one substitutes 2T instead of T. For even 
greater values of 8, the system must operate as if the re- 














petition period were 3T, 4T, etc., as 8 increases. 

For the realization of such still shorter transient re- 
sponses, it is necessary that the controlling device incor- 
porate several controlling devices, similar to those shown 





























on Fig. 3, a and b, but designed for the repetition periods 
T, 2T, 3T, 4T, etc. as a function of the maximum possible 























ee value of the external stimulus. The switching of one or 
— ona call yam " another of the controlling devices into the system must be 
implemented as a function of the magnitude of the initial 
Fig. 6b. error x(0), which equals the external stimulating step 6, 
2(t) and may be implemented by means of some supplementary 
































device. 
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quantity y(t). The total duration of the first and second, 
and also the third and fourth, pulses can be determined 
from the same expressions as in the case when 8 < x9, 
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TRANSIENT AND STEADY-STATE PROCESSES IN AN 
AUTOMATIC RANGE FINDER 


Il. RANGE FINDER WITH AN OPERATIONAL AMPLIFIER 


F. M. Kilin 

(Leningrad) 

(Translated from: Avtomatika { Telemekhanika, Vol.21, No, 2, pp. 209-219, February, 1960) 
Original article submitted June 16, 1959 


The dynamic properties of a range finder with an operational amplifier, including an integrating and 
an inertial link, are considered. In the general case, the investigation of the processes in a range finder with 
an operational amplifier, if with this one takes into account the step-function variations of the time-separator 
parameters, leads to cumbersom mathematical computations, The general investigation given in part I of this 
work [1] and the equations therein obtained for estimating the range finder’s properties permit these difficul- 
ties to be avoided, 


1. Properties of the Range Finder with an The equations for the controlled delay circuit (CDC) 
and the tracking error are 

AS (t) = cU,(t), Aa(t) + ADd(t) = Ad, (4). 

If we take the Lapace transforms of (1) and (3) and 
System solve them for L{A a (t)}, we get 

r 

An estimate of the effect of noise on a range finder L {(Aa(t)} = Aa (p) 5 {A9, (t)} —o >> A,» 7 U;, (0) 
of the simplest construction, not containing any operation- X (p) bos (P) 
al amplifier, was given in [2]. In the present work, we (4) 
start from the general-range-finder theory presented in [1]. 





Operational Amplifier When ItIs Investiga- 





ted as a Continuous Automatic—Control 














where 
Application of this theory to practically important concrete 20K, E 
designs permits their dynamic properties to be estimated X (p) = 4a (p) + =<" Ar (Pp). (5) 
comparatively simply. We shall not treat the passage of 
random signals through the range finder in this paper. In expanded form, 4a(p) and A,,{p) are given by 
In [1] we found the conditions under which the in- 
vestigation of the processes in the range finder.can be re- Aa (p) = II (p + axx), 
duced to the investigation of processes in the continuous k=1 
system described by (1.38)*. We shall assume here that A — 8 pproa te ZO a trl (6) 
these conditions are mett. We take a definite structure rk (P) = SrkP ~ . hes wt k + 
for the operational amplifier, for example, one for which on = | ; 
(1.38) take the following form. 1 for r=k, 
The equations for the time separator and the opera- niteaiiidis z/ rk) (J = 1,2,..,r-1) are constant real coef- 
tional amplifier are prey j — 
(D 4+- a,,)U,()) = 2Ko Es Aa (t), The form of the operator function X(p), defined by 
ToT (5), depends on the structure of the element Ay). In its 
Ay,U, (t) +- (D + age) U2 (t) = 0, (1) expanded form, the expression for A,(p) is given as 
Gg,U y(t) + age 2 (t) + (D + 453) U3 (t) = 9, An, P + Gee, 0,... ( 
431, 32, Pp + Q33,-+-+9 
a Ars (P) = (1 ee ee 
> OrmU m (t) + (D + Grr) U, (t) == 0. ary, are, Orgy. - + yr, r—1 
m=1 
The parameter a,, appertains to the time separator, = Z,pr-2+Z,p’-3+ ---+Zy—4, (1) 
and is expressed in terms of the time constants % and 
Ty by 
a T—a@ 





*The Roman numeral I in parentheses means that the 
ayy = Sy : (2) formula numeration refers to part I of the given work. 
t Cf. [1], section 9 and (1.44). 
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p) 
(4) 
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xy (8) 


coef- 


ed by 
In its 
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where the 2, (j = 1,2,...,r 
We introduce the notation 
Aa = —_ Aa, + _ 


L (Aa, (t)} = — > Tomy Ue, & 
A 
L {Aa (t)} = 80) cas, (t)}. 
We assume that 
AS, (t) = nsl*. (9) 
where the 1), are constant coefficients. Consequently, 
,4a(p) 

L {Aas (t)} = sl Xp) ° (10) 


If the range finder is a stable system.then, for the 
steady-state mode, we obtain 
Aa, = 0, 


Aa, (t) = s! Ne Z Ge k I 


est Ea rs rel. 
~ Rl Lap* X (p) 
For large values of t, the greatest t influence on the 
accuracy of range-finder operation will be manifested 
by the term 





te—* 
(11) 





II a4 
Adieg (t) = HGal*, Go = = “IKE * (12) 
TT yy + ay A,, (0 
k=] 10 


On the basis of (7), we have 


91, Qo2, Q, oer 0 


Ay, (0) = (— 1)" G31, G32, 33, .--, O 
, = (— 
ary, are, Apg, - 225 Oy, _ 
(13) 


It follows from (12) and (13) that the steady-state 
accuracy of the range finder will be increased if the co- 
efficients a,, are decreased and the quantity A (0) is in- 
creased. According to (13), the increase of A,(0) can be 
achieved by means of the operational amplifier para- 
meters ask for which ja k. 

To see more graphically how the coefficients a, 

(k = 1,2,..,r) affect tracking accuracy, we consider the 
case when the following equalities hold 


ay =0 (i=r—p+i, r—p+2 eee, r). f (14) 


For this case, in the steady state, we get 


v 


_ Gy . 
Aa, (t) = s! ns 2 C— hl tv-*, v= s—p. (15) 
The coefficients Gy, are found from the formulas 


a [ at Aa,(p) 4 
Guy= {4 To. Aa, (p)= Il (p + axx). 16) 


A=) 


~1) are constant real coefficients. 


Thus, if the conditions of (14) hold, the degree of the 
polynominal in t (15), determining the magnitude of the 
steady-state tracking error, is lower by p units than that 
in (11), which provides for a decrease in the rate of error 
growth for large values of t. Formula (12) graphically 
shows how the range-finder parameters and the parameters 
of the sequence of pulses affect the accuracy of the de- 
vice's operation. 


2. Range Finder with a Simple-Integration 





Operational Amplifier 





By setting r = 2 in(1) and (3), we obtain the equa- 
tions for the range finder with the simplest operational 
amplifier: 

2K oF 
(D + ay) U, (t) = 7" 4a (t), 
10 


ay,U; (t) + (D + aos) U2 (t) = 0 


AD (1) = oU,(t),  Aa(t)+ Ad(t) = Ad, (2). 


We shall be interested only in the component A og. 
We assume that the translation of the return pulses from 
the object is described by the function 


0 

AS, (t) -{ for t< 0, (18) 
| Not mt fort >0, 

where 19 andn, are constants. 


By solving (17) with (18) taken into account, we ob- 
tain the expression for A a(t): 


2 
Aas (t) = 9[Go+ >) Seem ]+ OO 


sn 


+h [Gt + Gs + 5S sent. 


kemy © 


(17) 


The coefficients Gj, and S; are found, respectively, 
from the formulas 


20K oF, an (a4, + a9) 





G, =“ G, 





oe ae 
oe eh alkt 
Se ae PyX” (By) ( ‘ ), 
where 
(ms) = {5 XP} 
Here, 
X(p)= p?+Xipt+ Xe, Xy = ay, + aoe, 
2 = 413099 — yg Qo. 


The parameters 1; and jig are the roots of the charac- 
teristic equation 


X(u) = 0 


If the range finder is a stable syetem, then, for the 
steady state, (19) takes the simpie torm: 


Aas (t) = Go [% + mt] + mG. (23) 


147 





The function A 9;(t) is related to the function Xt)= Xz, > at, /4, the roots of the characteristic equation assume 
Ry = Ryt which describes the motion of the object with complex values, and the tracking system goes into a mode 
respect to the radar set. This relationship is expressed by of damped oscillation. With this, the lower the value of 
yy, the lower the values of X, for which this transition 
occurs. From this point of view, it is well to have a larg 
value of a3;, which may be obtained by decreasing the 
time constants Ty and 7, (and not increasing them, as 
was required above). 
Consequently, (23) can be rewritten in the form: To decrease the steady-state error, it is required to 
give a large value to the quantity X,. However, for large 


AS, (t) = 20 | 


where c is the velocity of propagation of electromagnetic 
waves. 


9 
Aa, (t) = IR +- Rt] + —. (24) values of Xq, the oscillatory process will be characterize; 
‘ by an inadmissibly large overshoot. This is shown by the 
The value of Aa, corresponds to an error in the radius example, given in the Appendix, of the given range-finde 


scheme. 


AR (t) = 9% GR + Ril + RG. (25) 





Structure of an Operational Amplifier 





If there are very limited possibilities for decreasing Providing Changes in the Range-~-Finder's 
the rate of error growth in the range-finder with a sim- 
ple amplifier, such limitations do not exist in the given 
scheme. It follows from (19), (20), and (24) that, if ags It was previously shown that an increase of the para- 
tends to zero, one can attain an arbitarily small value for meter X, leads to an oscillatory transient response with a 
the quantity 2GyR,;/c which characterizes the rate of growtt large overshoot. This imposes a limitation on the choice 
of the error. For agg = 0, when the operational amplifier of values for X,. It also limits the possibilities of decrea- 
becomes a straight integrator, the following relationships ing the tracking error in the transient and steady states. 
hold: We now show that such limitations, and to a large degree, 

—_ _ 20K obs an G,=0, G= a11 also appear when an operational amplifier with multiple 
tioT integration is connected in the range finder; due to these 
limitations, the effectiveness of employing such an ampli- 
fier in similar devices is significantly decreased. We now 





Dynamic Properties 





i.e., for the steady state, (24) in the given case takes the 


_— R r write the characteristic equation of (5) in expanded form. 
ae oe = 0 (i = 2,3,..,r) in(5), we get 
Aa, > a oe (26) By setting aj; (i = 2,3,..,r) in (5) g 





To guarantee stability of operation, one must give a 
positive value to the parameter Xa, i.e., X_ > 0 or ag; >0. 
According to (26), when ag = 0 the tracking error in ov " 
= + any 1 ee---* 
the steady state is proportional to a,,, the parameter of X (p) = Pt Xap - 
the time separator. The value of this parameter is chosen 
in accordance with the conditions which guarantee reliable 


ask Aa (p) = p’—*(p +a), (28) 


X,-,p+X,=0. @ 


The coefficients X;, are found from the formulas 





operation of the time separator. As shown by (2), to de- X, = dy, 

crease parameter a,, it is necessary to increase both time 20K ok, 

constants Ty and 1,;, Which is not always possible. Thus, Xk = ToT" Ze-1 (k= 2,3,...,7), (MM 
even if the translations of the return pulses follow (18), the 
range finder now under consideration has essential disad- In their turn, the quantities Z;, Z,,..,2, are determit 
vantages. In it, the tracking accuracy is limited by the ed from (6) in terms of the operationalamplifier's para- 


possibilities of decreasing parameter ay, of the time sepa- eters bk (j = 2,3,...,15 k = 1,2,...,j). 

rator. In the case under consideration, the operational-am- 
Parameter ay, also imposes limitations on the character plifier implements an s-fold integration. With this, the 

of the range finder’s transient response. In fact, fors=2, Wantitys is related tor by the equation s = r-1. 


a = 0 and a given value of a4;, there remains at our dis- The character of the transient response depends on 
to X, by the equation to vary, within wide limits, all the coefficients X;, of (29), 


this will thereby provide the necessary conditions for chalj 
20K, E, A : : 
— (27) ing the character of the tracking system's transient re- 
trol’ sponse. These possibilities are easily established, on the 
In view of the fact that the device is described basis of (30). 
by a second-order differential equation, it is neces- Since the coefficients X_, Xs,..,X, depend on the op- 
sary to have at least two variable parameters. The rational amplfier’s parameters, they admit variations 
presence of only one of them entails this, that for within quite broad limits. This can not be said, howevet, 





X,=— 
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of coefficient X; = ayy, the magnitude of which is deter- 
mined by the parameter ay, of the time separator. As has 
already been stated, the possibilities of varying this para- 
meter are limited. 

On the other hand, the coefficients X; may be ex- 
pressed as follows: 


r 


X, = > M,, 
k=] 
X,=M, >) Me +My >) Me +---+M,_.M,, 
k=2 k= 
; (31) 
X, = M,M,...M,. 
Here, 
M, = —Pe (k=1,2,...,r) (32) 
Consequently, 
ay; = > M, => >) ReM,. (33) 


k=1 

Thus, if it is necessary to improve the transient re- 
sponse by increasing the real parts of the roots Ly, then, 
according to (33), it is necessary to increase the value 
of parameter ay;. In the case of multiple roots, 


M=M,=M,=.---=M, 


Formula (33) takes the simplier form: ay; = rM, or aq; = - 
(s+ 1)M, Le., it isrequired to provide a greater value for 
the parameter ay;, the greater the multiplicity s of in- 
tegration in the operational amplifier and the greater 

the necessary magnitude of the quantity M. 

The basic disadvantage of the scheme under con- 
sideration is that it is not possible in practice to provide 
the variation in the time separator's parameter a,,; which 
is required for improving the transient response in the 
automatic range finder. We now show that this disad - 
vantage is eliminated if we use a different operational- 
amplifier scheme, for example, one for which the pro- 
cesses are described by the equations: 


2K oE, 
(D + ay) Us () = 5 da (t), 


Ao,U, (t) + (D + aq) U2 (t) = 0, 
GU; (t) + Agql2 (t) + (D + a3) U3 (t), 





tr 


>) Orel x (t) + (D + arr) U, (t) = 0, (34) 


k=1 


>) orgs, Ux () + ras (t) = 0. 


k=1 
Equations (34) can be obtained from (1) by adding 
r 
Ura = — QD Ort eUn, 


k=1 


which expresses the simple amplification of the signals UL 
with gains proportional to the coefficients a,.4 ,., with 
summation of the results. For (34) we find 


r 
Aa(D) = |] (D + ax), Argi..(D) = 
k=1 
= 2D * 4+ 20 *+.--4+0_D)4+ 2... 
The equations tor the controlled delay circuit and the 
tracking error when the operational amplifier described by 
(34) is connected in the device are written in the form 


(35) 


AS (t) = oV,4; (t), Aa (t) + Ad (t) = Ad, (0). (34) 


Taking theLaplace transforms of (34) and (36), and 
solving them for L {Aa(t)}, we obtain 


Pin - A, +1, x (P) 
L (deh = —2 PFO, ay 





L (ha, (t)) = tL {0% (0), 


Aa, + Aas = Aa. 
For aii 0 (i = 2,3,...,r), we have 


Aa (p) = po (p + ay), 








X(p) =p +X pr + Xgpr te + Xr og 
2cK pL, _. 
xX, = 4+ —" Zo 
20K , 
X,= et e. @enks,....0, 


Here, all the coefficients X,,Xz,..,X, depend on the 
operational amplifier’s parameters. Consequently, it is 
not required to vary the parameter a,, in changing from 
one mode of operation to another, since this transition may 
be effected by a corresponding change of the operational 
amplifiers parameters. When an operational amplifier 
of the type described by (34) is used, the parameter ay, is 
not profitably increased, as was required in the examples 
given earlier,but, conversely, should be decreased, since 
the steady-state error is the lower, the smaller the val- 
ue of aq. 


4. Automatic Range Finder with Double In- 





tegration Operational Amplifier 





In using an operational amplifier with a double inte- 
gration, the processes of which are described by (34), the 
avtomatic range finder’s equations are written in the form 


(D + ay) U; () == a(t), 
AU, (t) + (D + ag.) U2 (t) = 0, 
3,U, (t) + agg 2 (t) + (D + aay) Us (t) = 0, (39) 
4,0; (t) + gaa (0) + aes 5 (t) + Us () = 9, 
Ad (t) = o, (1), 


Aa (t) + AS (t) = A9, (2). 
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We shall be interested below only in the component 
4 o. Taking the Laplace transform of (39), we get 





L (Aas ()) = IPL (09 (4)). (40) 


Here, Aa (Pp) = (Pp + 441) (Pp + ae) (Pp + ass), 
X (p) = p> + Xp? + Xap+ Xz, 
X, = ay, + Go2 + 33 + eZo, 


Xa = G4 (G22 + G33) + Go2033 + eZ, (41) 


X 3 = 11Ag2033 + €Zo, 
200K, EF, 
Ti9T ‘ 


The coefficients Z’,, Z',; and Z’, are found from the 
equations 


6 = 





Z=— Ga, 
Zy = 21042 + G31043 — G4; (22 + Gs), (42) 
Ze = AgiAg2M43 + Go (d53449 — 39043) — 041490035. 


We set 
AS, (t) = W,t + > W,t?, R(t) = Ryt+ +Ry" (43) 


or 
Ad, (t) = AO _ 2h, 4 Ss (44) 


where c is the velocity of propagation of electromagnetic 
waves. The function Kt) describes the translation of the 
object with respect to the radar set. 

Going to the time-domain function corresponding to 


(40) while taking (44) into account, we find, for the steady 


state, 
ars (t) = 8 7, (1) + 7, (n, 
Ns (t) _ Giot + Gi, 


h(t) = mi ott 4 Se + Goo, (t > to), 


where t, is the te defining the duration of the tran- 
sient response. 

The coefficients Gyand G3, are found from the 
formulas 


, =a da (p) iain 
Gur = Tilak FW le (i=1,2). (46) 


If we set agg = agy = 0, we obtain 


(45) 


fi (t)=0, fe= a (47) 


The most favorable conditions for increasing tracking 


accuracy occur when ayy =4 33 = 0, when the operational 
amplifier acts simply as a double integrator. 


Appendix 


Example 1, To determine the transient response and 
the magnitude of the steady-state tracking error for the 


150 








following automatic range-finder data: pulse-repetition 
frequency F= 500 sec™!, duration of the return pulses from 
the object | a = 240 Ssec, length of the tracking pulses 
O% = 3- 10sec, KE, = 60V, Ty = 2- 10 “sec, Ty, = 00, the 
operational amplifier is an integrator, i.e., agg = 0. The 
object moves relative to the radar set with a constant ye- 
locity of R, = 300m/sec. 
In correspondence with (2), we find that 


a 1 
— —— — 9 sec -.  ¥ =— = 0,002 
a1 T 710 Jv F sec 





The computational results are given in the first 
table: 

The table fives the steady-state values of A oy» and 
AR. In addition, we have computed the values of the rooy 
Hy and yg, from which one may derive a representation of 
the character of the transient response in the tracking sys- 
tem. The table also gives the values of the quantities 


pil =| Rep,T I, pol’ =|Re peT |, 
(48) 
eT =|Imp,T | =|Imp.T |, 


from which one may judge if (1.44) hold, as they do in 
the example under consideration. 


As we can see, for a}, = 5 sec’, there ensues a dampe 
oscillatory mode when X, > 6.25 sec *. Figure 1 shows: 
graph of ry transient response for the case when X;, = 
16.12 sec and X, = 45.73 sec *. The first case correspon 
to the roots #4, =2.5 4 im, the second to the roots /iy)* 
2.54 i207. 


For X, = 16. 12sec %, the overshoot reaches app.oximate- 
ly 14%; as the result of this, the maximum error ro 
105m and, in the steady-state, 93m. For X, = 45. 13sec, 
the overshoot reaches a significant size - about 69%. 
Transient responses with such overshoots are inadmissible. 

We note in conclusion that the tracking error must not 
exceed |Aa|< a. If | Aa |=ca», the return pulses from 
the object fall beyond the limits of the tracking pulses, 
and normal operation of the range finder is disrupted. 
From this point of view, the mode of operation with X,* 
2.25sec™ can not actually exist. 

Example 2. To determine the transient response and 
the steady-state tracking error of an automatic range 
finder with one integrator for the following data: F = 50 
sec”, a = 2°10~*sec, a = 3°10 sec, Ty = 10 “sec, 

Ty = © ; Ry = 300m/sec. 

This example differs from the previous one in that the 
time constant Ty has half the value it previously did. 
Determining the value of the parameter aj, we find that 


a 
ann SS 10sec oe 
Tio 


the increased value of ay, creating more favorable condi- 
tions for improving the transient response. Figure 1 shows 
the graph of the transient response for X, = 50sec™. For 
this value of X,, the overshoot is about 7%, i.e., is signifi 
eantly less than in the case with X, = 45.73sec™, consider 
ed in Example 1. 
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X,sec —* t4,sec —* Bz, SEC -1 eT eT oT Aa,, sec AR, 
2,25 —0,5 —4,5 0,001 0,009 0 4 -40-6 600 
5,25 fl 5 —3,5| 0,003 | 0,007 0 1,9 -10-¢ 285 
6,25 —2,5 —2,5 0,005 0,005 0 1,6 -10-¢ 240 
10,25 —2,5+i2 —2,5—i2 0,005 0,005 0,004 0,98-10-6 447 
15,25 | —2,54i3 | —2,5~—i3 0,005 | 0,005 | 0,006 | 0,66-10-¢ 4100 
16,12 —2,5+in —2,5—in 0,005 0,005 0,006 0,62-10-® 93 
45,73 —2,5+i2" | —2,5—i2z 0,005 0,005 0,013 0,22-10-¢ 33 
ba(t) 
v, |—e_ 17, TaTi-(a1 a 0 
i 63 
Ay =4573; 0,25 S—{r7t)D+! D D L. 
4, =50, Qy=10 a 
A, = 1612, a, =5 Q 
| ay; 
Fig, 2. 
. L . Tx, SEC a Qa: ay 
0 as wv 15 tsec 
Fig. L. 2-40-* | —0,03 |} 0,08 | —0,05 
Example 3. To determine the parameters of the opera- 2-10-° | —0,20 0,18 —0,0 


tional amplifier with a double integrator which provides 
equal values for all the roots of the characteristic equation, 
and to find the steady-state tracking error in the region 

4a 0.1 a as the object moves relative to the radar set 
with a constant acceleration of Py = 3g = 30m/sec*. The 
automatic range-finder data are: F= 400sec~', « = 107*- 
sec, KE, = 60v, o = 0.5-10~° sec/v. The computations are 
carried out for two cases: a)T » = 2°10 “sec, Ty = Wi 

b) Ty = 2°10 *sec, 1, = ow. 


The steady-state error is determined from the formula 





= al 2R2an 
da = Aa, =". (49) 
The inequality |A a|= 0.1 a will hold if 
X%.>Iu= 2Reay, 2RF (50) 





O,fac 0,fety 


In finding the numerical value of Xgq for cases " a” 
and “b", we get 


4 on for ty =2-10-*sec 
sec *, 


lea 
0,4—1_ for t= 2- 10-9 sec. 
sec *, 


. We now determine the operational amplifiers para- 
meters for the value Xs = 1.5Xg,. For this, we compute 
the value of € in(41). By substituting numerical values, 
we get 


120 . = for T1909 =2-10-* sec 
sec. 


c= 


12 = for t= 2-10-* sec. 




















For the case of multiple roots, we have 1; = ly = s= 
= -M. For this, the coefficients X,, Xg, and Xs are related 
to M by the equations 

X,=3M, X,=3M*, X,=M?. 


In the given example, the magnitude of M is found 
from 


(51) 


3 ; 
ews 1,82 sec -! for Xg—=4 sec ~’, 
M=Vi5X ~{ , 
V ” 0,85 sec-1 for Xgo = 0,4 sec ~. 


On the basis of (41) we have 


211 + G29 + G33 + eZ, = 3M, 
211 (@g2 + Gg) + Gaadgy + 8Z, == 3M?, 
411429433 + eZ, = M3. 


By substituting agg = agg = 0 in Eqs. (42) and (52), and 
then substituting the values of 2, Z';, and Z from (42) 
and (52), we get 


(52) 


— 3M 3M? 
a1 = a= » 4g1@g3 + One = =" 


(53) 


1,5X 
421232043 = -——. 


We have at our disposal six variable parameters of the 
operational amplifier, these being related to M and Xg by 
the three equations of (53). Consequently, three of the six 
parameters can be given arbitrary values. For example, 
by setting ag; = agg = 1, ag, = 0, we get 


a + @,—3M e _ 3M? 
wits re hee (54) 
1,5 Xo 1.5 Ry 
43 = OO ES ST 


” U, lock o£, 


The structure of the time separator and the operational 
amplifier for the given case reduce to the scheme shown 
in Fig. 2. This scheme deserves attention, since its tech- 
nological implementation is comparatively simple. 

The results of the computations using (54) are given in 
the second table. 
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DETERMINATION OF THE SECOND-ORDER MOMENTS OF 
VARIOUS COORDINATES OF AUTOMATIC-CONTROL SYSTEMS 
BY THE USE OF ELECTRONIC ANALOG COMPUTERS* 


L. A. Tel'ksnis 
(Vil’nyus) 


(Translated from: Avtomatika i Telemekhanika, Vol, 21, No. 2, pp. 220-223, February, 1960) 


Original article submitted September 29, 1959 


An engineering method is presented for the determination of second-order moments of various coordinates 
of automatic-control systems by the method of canonical representations of random functions, using electronic 


analog computers, 


In the investigation and design of contemporary auto- 
matic-control systems (ACS), an important role is played 
by the statistical analysis of their dynamic accuracy. The 
taking into account of random factors is necessary in de- 
signing aircraft, automatic-navigation devices, ship-roll 
dampers, for determining radar-system accuracy and artil- 
lery-gear accuracy, for solving information-processing pro- 
blems, for designing modern drilling equipment, coal-min- 
ing machines, automobile suspension systems, etc. 

Analytical methods of determing second-order mo- 
ments of various ACS variables [1-3], due to the large vol- 
ume of complicated analytical computations and trans- 
formations they entail, can be used in engineering calcula- 
tions only for solving the most elementary problems. 

The method of physical simulation of random process- 
es, using electronic analog computers for ACS statistical 
analysis [3] can be used only for the solution of a narrow 
Class of stationary pzoblems. 

The methods of mathematical simulation of random 
processes, using electronic computers [2, 3], although they 
theoretically permit one to solve stationary and nonstation- 
ary problems of ACS statistical analysis, are very cumber- 
some to use in practice because of the large amount of 
complicated analytical computation and transformation 
which must be carried out in the preparation of the solu- 
tioa as well as in the solution process itself. 

For the creation of a practical method of computing 
second-order moments of various ACS coordinates, it is 
advantageous to use the canonical representation of ran- 
dom functions [1, 6] 


X(t) = D Vie, (), 


v=1 


(1) 


where the V,, are the random, mutually uncorrelated co- 
efficients of the canonical expansion, with dispersions D,,, 
and the x (t) are the canonical expansion’s nonrandom co- 
ordinate functions. 
Thanks to the introduction of the concept of normal- 
ized coordinate functions 
7 


x, (t) = D? x, (t) (2) 


one can obtain expressions for the second-order moments 
of various ACS coordinates which are readily amenable to 
computation by uncomplicated special-purpose electronic 
computers. We now write these expressions: 

a) for the output dispersion 


Dy (t) = Diy? 0; 


v=] 


b) for the output correlation (autocorrelation) function 


Ky (7). => Oy: 


v=] 


c) for the cross-correlation function of the input and 
output variables at a given moment of time 


Dey () = 225 ¥5 0; 


v=] 


(3) 


(4) 


(5) 


d) for the cross-correlation function of the input and 
output variables 


Kxy (t,2) = >) 2° (0) y* (2). 


v=] 


Here, ye (t) is the solution of the differential equation 
Lip, dy) = M(p, 22) =1,2,..., 


and the expressions 


(6) 


m 
>> bsp, 
j=0 


L(p, t) = >) apt, M(p,) = 
i=0 


d 
P=7r’ m<n, 


are the operators of the differential equation which de- 
scribes the automatic-control system being investigated. 
By using Relationships (3) - (6), one can construct a 
computer for the statistical analysis of an ACS which 
makes it possible to determine the second-order moments 
*This paper was given at the All-Union Conference on 
Computational Mathematics and Computational Tech- 
nology, at Moscow in 1959 
Tt Cn the computation of normalized coordinate functions, 
Cf. [6]. 


153 























































































































































































































































































































wy a 
z(t) yy (t) Titi) Dry(t) 
/ 2 5 4 |} 
z(t) 
Z(t) yp(t) p(t) ypltlyyit) ty(tt) 
1 ~ 2 6 5 4 p= 
y, (Y 
a(t) yy (t) _¥te) Di (yy(t) Kry (tt) 
1 sa 2 6 5 ena 
z(t) 
Fig. 1 
of an ACS duringthe transient response, when random dis- 
- turbances act on it. The equations of the system, consist- 
~2 g(t) ing of object plus controller, are written in the form: 
Pru (t) > a - 
er omé. 
@ = a, — a8 — agd + X (8) ' 
8 = w + agB + asb, 











Fig. 2 Principal diagram of special adder. - 


of various coordinates of the ACS as it is acted on by ran- 
dom disturbances, basing the calculations solely on regu- 
lar signals, namely, normalized coordinate functions. ¢ 
Such a device is a special-purpose automatic electronic 
analog computer. Figure 1 shows the block diagram for 

a computer designed to determine second-order moments. 
It contains the following basic blocks: 1 is the block for 
the memory and for the program control of the computa- 
tion, 2 is the electronic model of the ACS being investi- 
gated, 3 is a squaring device, 4 is a specialized adder, 5 
is a multiplier, 6 is a device for introducing delays. The 
devices used in the computer for squaring, multiplying 
and introducing delays are standard functional blocks used 
in electronic simulation [4]. The memory and control 
block allows input of information (normalized coordinate 
functions x § (t) and the program for controlling the com- 
putation) which is written on magnetic tape. For solving 
(7), the electronic model of the investigated ACS is used. 
Summation of the series is implemented by means of a 
special adder (Fig. 2), which is a dc amplifier, operating 
as a scale amplifier, with additional feedback in the form 
of a controlled delay line (CDL). 

It takes several minutes, on the average, to compute 
the output dispersion or the cross-correlation function of 
the ACS’s input and output variables on such a machine. 
It takes some tens of minutes to determine the output 
autocorrelation function and the ACS's output variable. 

Example. As an example, we consider the problem 
of determining the dispersion of the output coordinate 
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. 1 
b= > [i(p+po)]—2 




















+- -—->——- - —— 


















































Q 1 F d 4 Ss «& 7 4 9 tsec 


Fig. 3. Curve of the dispersion of the output co- 
ordinate ¥ of the ACS described by differential 
equations (8). 


Calculation of the dispersion Dy( t) of the ACS's out- 
put coordinate ¥ for 0 =t=T = 10 sec, when the random 
disturbance X(t) with correlation function K,{t,r) = exp 
(-| t-r] ) acts on it, was implemented in accordance with 
the block schematic (the first line thereof) shown on Fig. 
1. With this, ten terms of the canonical representation 
of the random disturbance X(t) [6] were taken into accoull 
The time for the computation was 110 sec. Figure 3 gives 
the curve for the dispersion Dy(t). 

What has been presented here shows that the use of 
canonical representations of random disturbances, elec- 
tronic models of the system being investigated, and addi- 
tional uncomplicated electronic computing hardware make 















it possible to find, quite simply from the point of view of 
engineering practice, the second-order moments of various 
coordinates of a broad class of automatic-control systems. 


LITERATURE CITED 


[1] V.S. Pugachev, Theory of Random Functions and 
Its Application to Automatic Control Problems [in Russian] 
(Gostekhizdat, 1957). 

[2] J. H. Laning and R. H. Battin, Random Processes 
in Automatic Control [Russian translation] (IL, 1958). 

[3] A. M. Batkov, "The solution of one class of non- 
stationary problems of the statistical dynamics of auto- 
matic control systems using electronic analog computers,” 
Dissertation [in Russian] (1958). 


[4] B. Ya. Kogan, Electronic Analog Computers and 
Their Use in the Investigation of Automatic Control Sys- 
tems [in Russian] (Gostekhizdat, 1958). 

[5] I. E. Kazakov, “Approximate probabilistic analy- 
sis of the operational accuracy of essentially nonlinear 
automatic systems," Avtomatika i Telemekhanika 17, 

5 (1956). ¢ 

[6] L. A. Tel’ksnis, "Expansion of random functions 
in canonical series using electronic analog computers,” 
Trudy Akademii Nauk LitSSR, Ser B 4 (20) (1959). 


¢ See English translation. 








FINITE AUTOMATA. I 


M. A. Aizerman, L. A. Gusev, L. I. Rozonoer, 


I, M. Smirnova, and A. A, Tal’ 


(Moscow) 


(Translated from: Avtomatika 1 Telemekhanika, Vol, 21, No. 2, pp. 224-236, February, 1960) 


Originai article submitted June 10, 1959 


The views of the authors are presented on the basic content and problem areas of the theory of finite 
automata and on the interrelationships of this theory with the theory of relay-contact circuits; also considered 


are the possible ways of realizing finite automata, 


INTRODUCTION 


As soon as it became clear that the propositional 
calculus was adequate for the technical problems re- 
lated to the investigations of many relay (switching) de- 
vices [1-3], the abstract concept of the “nonsequential 
switching circuit” was formulated. For many switching 
circuits, such an abstraction turned out to be natural, 
and the propositional calculus became the basic mathe- 
matical apparatus for analyzing and synthesizing within 
the limits of this abstraction. However, this abstraction 
turned out to be inadequate, even for many important 
switching circuits, since it gave, in principle, no possi- 
bility of considering circuit operation in time, and only 
by artificial methods can be used for solving certain pro- 
blems in the theory of such circuits. 

The desire to introduce time (“sequentiality”) into 
consideration explicitly, and to find an adequate mathe- | 
matical apparatus for the description of circuits operating 
through time went beyond the limits of application of 
solely algebraic, logical operations (the propositional cal- 
culus). 

"Motion" in such devices was given by special 
formulas and tables (Cf., for example, [3-7]), by geo- 
metric forms (motion of a representative point along the 
vertices of multidimensional cubes, or along the nodes 
of some graph or another) (Cf., for example, 4,6-9]) 
and, finally, by integer recursive functions (Cf., for ex- 
ample, [9-14]). In these works, account was taken, to a 
greater or lesser degree, of the peculiar features of some 
concrete method or other of technically realizing logical 
operations (in particular, the specific features of relay 
contact technology), i.e., essentially, the soil was prepared 
for the introduction of new, broader abstractions. It gradu- 
ally became understood that what was at issue was the iso- 
lating of a specific class of dynamic systems % and that 
systems constructed of elements having, not two, but an 
arbitrary, though finite, number of possible states had to 
be encompassed by these new abstractions. 

While, in ordinary dynamic systems, the space and 
time coordinates have continua for their domains, the 
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basic idiosyncrasy of this special class of dynamic sys- 
tems is that their time has a countable set for its domain, 
while their coordinates have finite sets for theirs. 


As these abstractions were developed, various names 
were given to them: logical nets [10], generalized nerve 
nets [15, 16], finite automata [16], nonprimitive circuits 
[14], sequential machines [18, 19], robots [17] etc. 


The authors prefer the name finite automata for 
these abstractions. 





Certainly, not all switching devices (i.e., those ex- 
hibiting relay functioning) are encompassed by these ab- 
stractions. For example, the abstraction "finite automata’ 
is not adequate for systems containing devices in which 
the coordinates are chosen from finite sets, while the 
time has a continuum for its domain, although these are 
relay devices. On the other hand, the abstraction "finite 
automata” corresponds organically to all relay contact 
circuits which are devoid of such devices. But with cer- 
tain natural idealizations, many other important techni- 
cal devices also correspond to it, for example, machine- 
stand automata, mass-production lines, various digital 
machines, automata for piecework production, certain 
abstract models of nervous systems (for example, the 
nerve nets of McCulloch and Pitts) ,etc. 


2. Basic Definitions and Terminology 





As was said earlier, a finite automaton is a special 
type of dynamic system. The first characteristic of this 
class of dynamic system is that, in describing its actions, 
only discrete moments of time are of importance, the 
intervals between these moments (which, in general, are 
of unequal duration) being called cycles of the finite 
automaton's operation. Let p be one of these moments 
("the present moment"), and let p-1, p-2,...,2,1,0 be 
the series of preceding moments of time. In this notation, 
0 corresponds to the initial moment of automaton opera- 
tion. We assume that some conditions or other are given 


* In this sense, the older abstraction - the nonsequential 
relay contact circuit - corresponded to a static system. 
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which unambiguously define the decomposition of the 
time axis into cyclest. 

In the general case, the dynamic systems under con- 
sideration are subject to external stimuli. In the sequel, 
these stimuli will be called inputs. 

Definition. A finite automaton is adynamic system 
which, at the discrete moments of time under considera - 
tion, satisfies the following conditions. 

1. At each of these moments of time, the system 
can be in just one of a possible finite number k of states. 
2. At these moments of time, the system's inputs 
are also restricted to be chosen from a finite number r 
of possible states. > 

3. At any of the moments of time considered, the 
state of the system is uniquely defined by the state of the 
system and the state of the input at the previous moment 
of time. 

To distinguish one automaton state from another, we 
assign a different symbol to each of the k possible states. 
We denote these symbols by x4, x ,..., KR, We call the 
set of these symbols, denoted by {x} » the alphabet of the 
states. Analogously, we write the symbols p4, pg,.--»Pr 
for the input states. They form the automaton’s input 
alphabet{ p }+. 

We denote byx(t) and p(t) the symbols for the auto- 
maton state and, input respectively, at time t. Then, by 
virtue of point 3 in the definition of the finite 
its operation is described by the recursion formula 

x (p) = F [x(p—1), p(p— Mls @) 
where F is a single-valued function which puts a definite 
symbol from the state alphabet into correspondence with 
each pair of symbols, one of which is chosen from the 
state alphabet, the other from the input alphabet. 

Thus, the abstraction introduced by the concept of 
"finite automata” singles out a class of systems whose pro- 
cesses are described, notby differential equations, nor yet 
by difference equations of general form, but by specific 
equations of the type given in(1). The theory of finite 
automata Concerns the general properties of such a type of 
system. Therefore, from the point of view of this theory, 


which we call the output alphabet, Then, to describe the 
operation of an automaton with outputs, (1) must be supple- 
mented by the “instantaneous” relationship 


h(p) = O[x(p)], (2) 


which puts a symbol from the output alphabet into corre- 
spondence with each symbol from the state alphabet. 

Equation (1) can be given by Table 1, which we call 
the fundamental table of the finite automaton. 
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This table is set up in the following way: a cell of 
the table is defined by giving one symbol each from the 
{p} and {x} alphabets; by taking these symbols to re- 
present matters at time p-l, i.e., by considering them to 
represent p(p-1) and (p-1), we determine, (p) from (1); 
we assign the value of, (p) thus determined to this cell. 

If the automaton has outputs, then (2) defines the 
table of correspondences between, and A (Table 2). 








TABLE 2 
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t It may be assumed, for example, that the continuous 
time axis is previously divided in some manner into cy- 
cles which depend neither on the automaton’s operation 
nor on its external stimuli. In other cases, the moment 
when a new cycle begins can be taken as the moment 
when some change in the automaton’s input occurs. Fin- 


all systems are identical whose descriptions, with the appli- ally, the cycles may be defined in such manner that the 


cable idealizations, reduce to an equation of type (1) with 


beginning of a new cycle will also depend on the state 


the same function F in the right member. Moreover, auto- of the automaton, or on its output (Cf. part II of the pre- 
mata described by equations of type (1) with different func- sent paper, to be published in the next number of this 
tions F in their right members can still turn out to be equiv- journal). 


alent, in a definite sense, to one another (Cf. section 4 be- 


low). 

In the actual systems subsumed under the heading of 
"finite automata,” a change of state of an automaton can, 
in some cases, be taken directly as being the result of the 
system’s operation. In other cases, the automaton has an 
“output ,” whose state is uniquely determined by the auto- 
maton's state, but is not identical with it**. 


Of course, in a finite automaton, the number of possi- 


ble output states may, but in general does not, coincide 


+ In particular, one can number the automata and input 
states. Then, the state and input alphabets would be com- 
prised of sets of numbers, 1,2,...,k and 1,2,...,r. Such 
alphabets might be convenient in investigations of in- 
dividual automata. They are inconvenient in those cases 
when automata are joined, one to another, since the same 
symbol here, say the number 3, would apply to states and 
inputs of different automata. 

** For example, in an automaton used for sending alarm 
signals, the important thing is not its state at each mo- 


with k. Let the output have J (2 <k) states. We write each ment, but only whether this state does or does not corres- 


of the possible output states by a symbolA{Aa Age.» AZ}, 


pond to an alarm signal. 
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Together with the finite automaton’s fundamental 
table (Table 1), we shall consider an infinite tape (Table 
3), characterizing, for a given initial state and a given 
sequence of input symbols, the sequence of appearance 
of state symbols as determined by (1) (or, what amounts 
to the same thing, as determined by Table 1). 
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On this tape, each moment of time t is put into 
correspnndence with a p andax«- We call it the state 
tape. 

A triple of symbols, x(p-1), p(p-1) and (p), such 
as shown on Table 3 by the heavy lines, and chosen for 
any p, is related internally by(1). We call such a triple 
a triad. 























TABLE 4 
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In cases when the system is defined by (1) and (2) (or 
by Tables 1 and 2), we shall also consider the automaton's 
output tape (Table 4), in which each moment of time t is 
put into correspondence with a pair of symbols p and A. 

We take particular note of the case when the input 
state does not change with time - automata of such type 
are Called autonomous. For an autonomous automaton, 
(1) has the form: 


x (p) = F[x(p— 1), p (0)], (3) 


where p(0) is a time-invariant symbol which can be con- 
sidered as a parameter. 

The operation of an autonomous automaton ean be 
illustrated by a graph. 

We place k points on a plane, assigning to each point 
one symbol corresponding to one of the automaton's possi- 
ble states. Equation (3) then defines the sequence of 
changes of state, i.e., the directed transitions from point 
to point. These directions can be indicated by arrows. It 
is easily seen that, in general, several arrows can lead to 
the same point. From a given point, however, at most 
one, and possible no arrows may lead. 

The set of k points, together with the arrows as just 
described, is called the graph of the autonomous finite 
automaton. 

In a nonautonomous automaton, there can be r input 
states. Therefore, by fixing and by leaving thenceforth 
invariant the input state, one can organize r different 
autonomous automata from a given nonautonomous auto- 
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maton, each of the autonomous automata having its own 
graph. The set of these r graphs completely characterizes 
the nonautonomous automaton. Figure 1 shows an example 
of such a set, comprised of three grapits; a state table 
(Table 5) corresponds to this automaton. 
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3. The General Problems Arising in Con- 





nection with the Introduction of the Ab- 





straction “Finite Automata” 





Although the abstraction of “finite automata,” as 
this term was introduced in section 2, is a far-reaching 
one, there arose with its introduction [i.e., with the iso- 
lating of the class of systems described by (1), or by (1)+ 
(2)] a number of problems, equally important for all ac- 
tual systems for which this concept is adequate. 

One of these problems, that of the representation of 
events in finite automata, was solved (inefficiently, it is 
true) by Kleene's theorems [16] (Cf. also [19-21}). 

This is not the sole problem, and a number of im- 
portant problems in finite automata theory are today not 
only unsolved, but not even rigorously formulated. These 
problems arise principally in connection with the fact that 
both the state tape and the automaton's output tape are in- 
finite. It is therefore meaningless to say that these tapes 
are “given.” All that may be given are algorithms which 
permit, in a finite number of steps, choosing of the sym- 
bol from given alphabets which must be written in an ar- 
bitrary cell of the tape. 

We denote by Ap and Ax the algorithms correspond- 
ing to the upper (p) and lower (x) rows of the state tape, 
and by A, and A) the analogous algorithms for the out- 
put tape. From this point of view, a finite automaton, de- 
fined by (1) or by (1}+(2) (i.e., by Table 1 or Tables 1 
and 2), is essentially an operator which establishes the 
relationships between algorithms A p and Ax, or between 
Ap and A). 

In connection with this, no scientific problems arise 
in the task of analysis, i.e., in the task of determining, 
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for a given automaton, the algorithm Ax (or A)) from 

Ap and conversely, and here one may only speak of de- 
veloping more or less convenient ways of implementing 
the operations defined by the “finite automaton” operators. 

Matters are otherwise when the tape algorithms are 
given (i.e., Ap and Ax, or A. and A)) and it is required 
to synthesize the automaton (i.e., to construct (1)). 

We first consider the case when it is required to syn- 
thesize the automaton from the given algorithms A. and 
Ax It is clear that this problem cannot always be solved, 
since the tape defined by algorithms A) and A, may con- 
tain triads which are mutually contradictory. 

We call two triads contradictory if, in these triads, 
the symbols p(p-1) and, (p-1) are identical, but the sym- 
bols for ,(p) are different. For example, the pairs of 
triads 








ly permit one to fill in all or some of the cells of Table 1, 
and to define the automaton with an exactitude to within 
the cells which remained unfilled and could be filled in 
arbitrarily. 

In the case when a finite number of state tapes is 
given, and it is required to find an automaton which re- 
produces all these tapes, no new problem arises, since 
the algorithm for observing contradictions and the al- 
gorithm for determining all the different triads may be 
applied to each tape in turn. New problems arise only 
when the number of state tapes is infinite. 

The matter is significantly more complicated when 
it is required to synthesize an automaton from its output 
tape. Here the additional problem arises: knowing the 
algorithms Ay and A), to determine k, i.e., the num- 
ber of states which must be realized in the automaton, 
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are noncontradictory. 

Each cell of Table I contains one symbol, There- 
fore, the existence of contradictory triads is incompatible 
with Table 1, i.e., it contradicts the single-valuedness 
of the function F in (1). 

Thus arises the first problem: knowing the algorithms 
A p and Ax, to find an algorithm for recognizing their 
“noncontradictoriness,” i.e., the absence of contradictory 
triads on the state tape. 


This problem may be formulated precisely, as follows. 


We are given the recursive functions x(p) and p(p), each 
of which can assume a finite number of values. It is re- 
quired to construct a recursive functional ¥[{ «(p),p (p)] 
(Cf. [23] Chap. 9) which equals unity on those functions 
x (p) and p(p) for which there exists the function F(u,v), 
defined on the finite set constituting the domain of u and 
v such that, for any p, the relationship x(p) = F[x(p-1), 
p(p-1)] holds, and which equals zero otherwise. If, as 
may be assumed, the problem as formulated is algorith- 
mically unsolvable, the necessity then arises of posing the 
more restricted tasks of determining a recursive function- 
al ¥ on the functions x (p) and p(p) which lie in a more 
restricted class than the class of all the recursive func- 
tions which assume the given number of values. 

Although the automaton’s state tape is infinite, the 
numbers r and k are finite, so that the total number of 
different triads on the tape is also finite. After the ques- 
tion of the “noncontradictoriness” of the algorithms is 
resolved, there arises a second problem: for given “non- 
contradictory” algorithms A, and A,, to find an algorithm 
for seeking all the different triads contained on an infinite 


state tape. The solution of this problem would immediate- viously used the “instantaneous” relationship of (2). 








and to find the function # in (2), i.e., from the given 
output tape to construct a noncontradictory state tape. 

It is precisely problems of this type which are most 
interesting in practice, since the requirements on the 
automaton to be synthesized are ordinarily given direct- 
ly in “input-output” terms. 

It is necessary to mention that the general problem 
of synthesizing finite automata from given operating con- 
ditions has thus far not been precisely formulated. The 
difficulty of a rigorous posing of this problem stems from 
the following. In the majority of questions arising in 
practice, the automaton must transform a given infinite 
class of input sequences into a required class of output 
sequences. Up to now, however, there has been develop- 
ed no sufficiently general and convenient mathematical 
apparatus in terms of which an adequate description of 
these infinite classes of sequences could be given. In 
the present case, a description in terms of sets of tables 
[15] is unsatisfactory, since the original (verbal) given 
conditions of operation are posed in other ("logical") 
terms (closer to the language developed in [22]), while 
no algorithm has been established for making the transi- 
tion from the language of sets of tables to the “logical” 
language. 

All algorithmic problems of this nature vanish, and 
the synthesis problem becomes much easier, in the case 
when the automaton's tapes are finite. 


4. Mappings of Finite Automata and 





Their Equivalence 





In considering an automaton with outputs, we pre- 





This relationship can be considered as a rule for trans- 
forming state symbols to output symbols. An abstract 
device which realizes (2) is called a symbol transformer. 
There is, of course, an analogous symbol transformer at 
the automaton's input; it will transform the symbols of 
some external stimulus alphabet to symbols of the auto- 
maton's input alphabet. 

We now assume that two automata are given: auto- 
maton A, operating with input alphabet { a} and with 
state alphabet {a}, and automaton B, with input alpha - 
bet {8} and state alphabet {b}. We further assume 
that there are two symbol transformers: one of them 
(the input transformer) transforms the symbols { a} to 
{8}, and the second (output) transforms symbols { b} 


to {a}. 
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Fig. 2 

We now consider the system consisting of the input 
symbol transformer, automaton B and the output symbol 
transformer (Fig. 2). This combined system has inputs 
from the alphabet { a} and outputs from the alphabet 
{a}, the same as automaton ATT. a 

If we can choose input and output symbol trans- 
formers such that, for any sequence of input symbols 
from alphabet { « } and for any initial state of auto- 
maton, A, there would exist an initial state of automa - 
ton B such that there would appear, at the output of the 
combined system (automaton B and the two transformers), 
the same sequence of symbols from the alphabet {a} as 
the sequence of state symbols in automaton A, then we 
shall say that automaton B is a mapping of automaton A 
(or that automaton A is mapped onto automaton B). 

We symbolize this fact by: 


ACB or BDA. 
If simultaneously 


ACB and ADB, 


then the automata A and B will be said to be equivalent, 
and we symbolize this circumstance by A = B. 

We now give a precise definition of the concept of 
mapping. 

The automaton A ({ a}, {a}) is mapped onto auto- 
maton B({b}, {8}) (and automaton B is a mapping of 
automaton A) if, to each pair (al, oy) of symbols from 
the alphabets { a} = {aj,..,a,} and {a} ={ oy,..,0,}, 
one can set up a correspondence of such a pair({b;}, 
{8,}) of sets of symbols from the alphabets {b} = { b,,.. 
b,} and {6 } = {8,,...8,} such that, if 











F 4 (4, 4;)=4,, (4) 
then; for any symbols 6, € {ob} and Bp € {Bj} the fol- 
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lowing condition holds 


(5) 


Here, Fa and Fp are the functions in the left 
bers of equations of the type of (1) for automata A and B 
respectively. 

Two equivalent automata do not differ practically 
from one another — their differences are external, con- 
sisting only in the methods chosen for coding the automata 
states and inputs. In spite of this, two equivalent automata 
can be described by completely different equations of type 
(1). Even the variables in these equations may be taken 
from different alphabets. Thus, for example, two auto- 
mata with identical numbers of states and inputs, differ- 
ing only by the notation of the states and the inputs, i.e., 
the alphabets from which the symbols are taken and the 
manner in which these are associated with states and in- 
puts, are known to be equivalent, although the equations 
of type (1) for them are different. 

If A C & but AB, this means that there exists a 
symbol transformer by means of which automaton B can 
be “forced” to operate like automaton A, but not con- 
versely. 

The concept of the mapping of finite automata, and 
of their equivalence, is closely related to the topological 
peculiarities of their graphs. 
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As an example, we consider the graphs of the four 
autonomous finite automata, given in Fig. 3. The auto- 
mata of Fig. 3, b, and c are equivalent, since they differ 
only by the nomenclature of the states. The automaton 
of Fig. 3a,is mapped onto each of the automata of Fig. 3, 
b, and c. This becomes clear by supposing that a trans- 





tt To be sure, such a combined system might not be a 
finite automaton, since its output symbol at time p might 
not be uniquely determined by its output and input sym- 
bols at time p-1. We shall assume that the combined 
system under consideration is a finite automaton. 
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former is placed at the output of the automaton of Fig. 3,b 
which places each symbol b; (i = 1,...,8) into correspon- 
dence with any of the a; (i = 1,...,8). 

All three of the automata thus far considered are 
mapped onto the automaton of Fig. 3,d. The output trans- 
former which it is necessary to append to the automaton 
of Fig. 3,d for it to operate the same as the automaton of 
Fig. 3,b must place each of the symbols c; into correspon- 
dence with b; for i = 1,2,6,7,8,9,10,11, and 12. Symbols 
Cs and cg are made to correspond with the same symbol 
bs, symbols cg and c’, with by and symbols cg and c", with 
bs. In addition, symbols C4g,Cyg,...,Cy are placed into 
correspondence with any of the symbols by,...,bg. 

The problem now arises: knowing the fundamental 
tables of two automata, or their equations of type (1), to 
determine whether these automata are equivalent, or 
whether one of them is mapped onto the other¢t. In 
other words, it is required to find criteria for the mapp- 
ings and equivalences of automata. 

The concept of equivalence of automata, and of their 
mappings, allows one to introduce the concept of the ab- 
stract aggregating (assembly) of finite automata, i.e., of 
assembling a finite automaton from standard elementary 
automata. 


5. Abstract Assembly of a Finite Auto- 





maton 


We now consider a device with s input wires u,..., 
u, and n generalized coordinates X,,.....X,. At the times 
0,1,2,3,...,p, each input and each coordinate may have 
only one of a finite number of values, which we shall 
symbolize by symbols from some alphabets. At any mo- 
ment of time, each coordinate is determined uniquely 
by the values of all or some of the inputs and coordinates 
at the preceding moment of time. 

The operation of such a device is described by the 
-ecursion relationships 


xi (p) = filzi(p— 1),...,2n(p — 1); 


u;(p—1),...us(p—A))(i=1,.....n) (6) 


or, in vector form, 
x (p) = /(x(p—1), u(p— 1)). (7) 


If an xj can have x; different values (i = 1,...,n), 


n 
then the vector x(p) can have k = I] yi different val- 
i=1 
ues. We call the vector x the state of the device being 
considered, and we choose its k values from the symbols 
K1,...%, Of the alphabet {x }. Analogously, the input vec- 


s 
toru, with r= I] F values (where & 4 is the number 
j=1 
of values of uy) is called the input state, and we relate 
its values to the symbols p,,...,pr of the alphabet {p}. 


With this notation, (7) reduces to (1), which means that 
the device just described is a finite automaton. 

Equations (6) can also be understood as another form 
of equation of the type (1). Different equations of type 
(6) can correspond to one and the same equation of type 
(1), but they all describe the processes in equivalent 
automata. 

In the definition of a finite automaton, given in sec- 
tion 2, and in the equation of type (1) which correspond - 
ed to it, we had to do only with automata states and in- 
puts. Now, in making the transition to equations of type 
(6), we introduced coordinates and input wires (input 
lines). We shall call the transition from (1) to the equiva- 
lent (6) a choice of structure for the finite automaton, and 
(6) will be called its abstract structure (AS). A technical 
device which, with the suitable idealizations, is described 
by equations of type (6), will be called a technical reali- 
zation of the automaton’s structure. 

Let there now be given several systems of equations 
of type (6), each of which is the abstract structure of a 
corresponding automaton. Each of these systems of equa- 
tions has its own coordinate and input-line alphabets. 

We now introduce into the discussion symbol trans- 
formers, each with several input lines and one output. 
The symbols from the transformer’s output alphabet ap- 
pear at its output at the same instant that the symbols 
from their respective alphabets are applied to the input 
lines. 
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We use transformers as means for the interconnection 
of AS (as, for example, on Fig. 4), choosing them in such 
manner that the corresponding AS alphabets and transformers, 
are properly linked. We call such a set of AS, intercon- 
nected via symbol transformers, a net. The coordinates of 
the net are comprised of all the generalized coordinates x; 
of all the AS, including those not applied to any transform- 
er. The net's input lines are both the AS input lines on 
which no transformer acts, and the transformer input lines 
on which no AS coordinate acts. 

In order to obtain the equations for a net, it is neces- 
sary to use the given equations for the transformers for 
eliminating from the systems of equations of type (6) the 
corresponding variables for the AS input lines, replacing 
them by the AS coordinates and the variables of the trans- 
formers’ input lines. As the result, the net as a whole is 
also described by a system of recursion relationships of 


tt Certainly, what we have in mind here is not the tri- 
vialsolution, obtained by the exhaustive trial of all mean- 
ingful correspondences of symbols. 
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type (6), i.e., the net itself is the abstract structure of 
some finite automaton. 


From what has been said, it follows that a finite auto- 


maton can be constructed from other finite automata by 
joining them in a net. The question naturally arises as 
to the possibility of creating a set of automata, or even 
all meaningful automata, from a small collection of ori- 
ginal automata - “elements,” Such construction of auto- 
mata we Call abstract assembly. The technological value 
of this concept is obvious. 

The introduction of the concept of abstract assembly 


leads immediately to the posing of the following problems: 
to construct algorithms which answer the question as to the 


possibility of constructing nets from a given set of AS and 
transformers, these nets being mappings (in the sense of 
section 4) of given automata, and to specify an algorithm 
for the efficient implementation of this construction. 


A set of AS and transformers from which one can con- 
struct nets mapping any given automaton is called a com- 


plete set. The problem then immediately arises of deter- 
mining criteria for the completeness of a set. 

A complete set which is very important for techno- 
logy is the set consisting of the following elements: a) 
a binary delay element, described by the following equa- 


tion x(p) = 2(p-1) 


where x and z are given in an alphabet which contains 
only two symbols (for example, 0 and 1); b) some set of 
logical elements which permits the realization of any 
logical function (this set may consist for example, of con- 
junction, disjunction and negation elements). 

From this set, a net which is the mapping of any 
finite automaton can be constructed. Indeed, let a finite 
automaton be given by Eq. (1), i.e., let the function F 
and the alphabets {4} and {p} be given. 

We choose we numbers n and s such that the condi- 
tions 2 =k and 2 =r are satisfied. We introduce the 
variables Xty--eXp and uj,...,U,, Where each of these vari- 
ables assumes only the values 0 and 1. 

We consider Table 6. On its left side are enumerated 
all possible combinations of all values of Xpy--eeXrs Ugy- eee 


us, i.e., Table 6 has 2™8 rows. Each different combination 


of variables x;,..,X, is placed in correspondence with its 
own symbol from the alphabet {x}. Since the number of 
these symbols is k = 2", there may remain combinations 
X4,---sX, for which the alphabet {~} cannot supply sym- 
bols different from other symbols already used. We as- 


162 





sign arbitrary.previously used symbols from { x} to these 
combinations Kjy-+-yXp- 

We proceed analogously with the combinations of 
variables u,,..,U,, designating them by symbols from {p}. 
We write the selected symbols from {x} and {p } in the 
corresponding row of columns« (p-1) and p(p-1) of 
Table 6. 

Thus, to each set of variables Xgye-yX pi Uy. -U there 
corresponds a pair of symbols from the alphabets {x } and 
{p }. By taking these pairs of symbols as standing for 
K(p-1) and p(p-1), we find x(p) from (1), and write these 
values in the corresponding rows of column, ( p) of Table 
6. We can then immediately fill up the right side of 
Table 6, since we have already related a definite combi- 
nation of values X,,...,.X, with each symbol from { }ees 

The table thus constructed determined the value of 
all x{ at time p from the given x; and u, at time p-l, 
i.e., it gives n logical recursion relationships of the type 
of (6): 


vi (p) = Li [xz (p — 1),...,2n(p— 1); 


u,(p—1),...Us(p — 1)] 
(i =4, .. . si). 


(8) 


For any combination of values of the variables, the 
values of all the functions L; are read off from Table 6. 

It is easily seen that the obtained equations (8) are 
the abstract structure of an automaton which maps the 
original automaton which was given by an equation of 
type (1). 

Each equation in (8) can be realized in the form of 
an elementary one-cycle lag, with equation 


ri(p)=2(p—1) i= 4,--- 0) 

and a logical,instantaneously acting transformer 
a = Ly (24, ..-5 Zn; aS (ia i,... np 
(10) 


which, in the given case, is also a transformer of the co- 
ordinate x; and input uj symbols to input z, symbols of 
the lag of (9). In the given special case, all the symbols 
Xin Uy and Z, are chosen from the same alphabet, name- 
ly, { 0,1}. The entire set of equations (8) is a net com- 
posed of n lags of the form (9) and n transformers of the 


*** If k < 2" and if several combinations, x,,...,X,, Were 

related to one and the same symbol of {x }, then we choos 
any one of them in any case when the corresponding sym- 
ool arises. 
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form (10). Since any logical transformer is realized from 
the set of logical elements given by us in the conditions 
of the problem, this shows that this set of logical ele- 
ments, supplemented by a delay element, is a complete 
set of the elements necessary for constructing any auto- 
maton. 

It can be shown analogously that a set of elements 
allowing any function of q-valued logic (q > 2) to be 
reproduced, and supplemented by a delay element of 
the form of (9), wherein the input and output assume 
values from a q-symbol alphabet, also comprises a com- 
plete set. 

In this section we first used the concept of a delay 
element. This concept requires the following additional 
elucidation. 

When any symbol from the input alphabet is applied 
to the input of delay element (9), the corresponding sym- 
bol from the output alphabet appears at its output after 
a one-cycle delay. It was stated, at the very beginning 
of this paper, that the cycles could be nonuniform, and 
that the division of the time axis into cycles was uniquely 
determined by some additional conditions. So long as we 
had to do with abstract structure and abstract assembly, 
i.e., essentially with the replacement of (1) by equations 
of the form (6), this circumstance was of no significance 
since in (1) and (6) the sequentiality (i.e., the character 
of the variation of the argument p) was considered as 
given. In using the complete sets, described above, for 
the technological realization of automata, it is assumed 
that the automata contain special devices - sequencing 
transducers which, in correspondence with given condi- 
tions, process the signals for cycle beginning and apply 
them to delay elements. This stimuli of the delay ele- 
ments are not external in the sense in which the term 
"external stimuli” was used above; they are not encoded 
by symbols from {p } in setting up (1), and are not ger- 
mane to the stimuli uj in(6). Therefore, even an auto- 
nomous automaton may contain these stimuli, which are 
essential only for the technological realization, and are 
not taken into consideration when one deals with abstract 
structure or abstract assembly. 
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Methods are presented for quadratic approximations for functions given either analytically or graphically, 
Simple formulas and relatfonships are derived for finding the rule for argument partitioning. For obtaining a quad- 


ratic relationship, a variant of a diode element is considered. 


The paper describes the circuit of a universal function generator using quadratic approximations and capa- 
ble of being implemented in the form of an attachment to a type EMU-8 analog computer. 


Universal diode-function generators [1-5] based on 
the principle of piecewise-linear approximations are wide- 
ly used in modern electronic analog computers for repro- 
ducing various nonlinear functions. 

When it is necessary to reproduce a function with a com- 
paratively large range of variation of the absolute magni- 
tude of the first derivative, function generators of this type 
require the use of a large number of diode elements, which 
complicates the circuit and increases the error. In addition, 
discontinuities in the value of the first derivative of the re- 
produced function are introduced which, in many cases, is 
highly undesirable. 

The disadvantages just cited can, to a large extent, be 
removed if nonlinear functions are approximated by seg- 
ments of curves, in particular, by segments of parabolas 
[6-8]. 

However, the use of quadratic approximations entails a 
significant increase in the amount of computational work 
necessary to determine the partitioning along the argument 
axis, which has resulted in work being carried out to de- 
velop simplified analytic and graphical methods of calcu- 
lating functions, and a variant of a universal function gen- 
erator, based on the principle of quadratic approximations, 
has been developed. 


Method of Piecewise-Quadratic Approxi- | 





mation of Nonlinear Functions 





The problem of the piecewise-quadratic approximation 
of a nonlinear function {x) amounts to the finding of a num- 
ber of fixed argument values xq, Xg, X3, | which provides 
a second-order approximation polynomial with vertices 


f(X1), {X%q),--. 


P, (xz) = a, + ayx + >} bi (x — 23)’, 


i=l 
where 
ae (; for <u, 
i= 
B;= const for «>, 
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Fig. 1 


to the curve f(x) to be reproduced such that, on the entire 
interval of argument variation, the absolute error of the 

approximation function does not exceed some previously 
given value € (Fig. 1). 

With respect to the function f(x), it is assumed that iti 
single-valued and continuous on the given interval of arg 
ment variation, and that it has a continuous second deriv 
tive everywhere on that interval with the possible excep- 
tion of a finite number of first-order discontinuities. 

For quadratic interpolation, the remainder term in New 
ton’s interpolation formula has the form [9, 10): 

R(2) = £9 @—2,) (2 —2,) (x — 25), 
where f**"(€) is the value of the reproduced function's 
third derivative at some point € inside the interval (xj, 
Xs). The modulus of the error Xx) has its maximum val- 
ue when f*"* (E) and the expression | (x-x,) (x-X_) (x-%)| 
simultaneously attain their maxima: 


| R (x) | max == 
” (E)| 
—_ fee (x — x4) (x — 22) (x — Z3) |max. 
We introduce the new variable u = (x-x,)/(xg-x,) and 
let k = (X_-X) /(Xg-Xy) and h = Xg~Xy. 
Then, 











rl 


ntire 
the 
usly 


at itis 


f ange: 
Jerive 


cep 


n New 


's 
(Xy, 
1 val- 


x-1)| 


ax. 


) and 





17” (€) | max 


|R (u, k)| max = ——,—— | 


| 
3 


u(u — 1) (u—hk) |max hf’, 


The function Wu-1) (u-k), in the region of variable 
variation 0<u=1, 0=k=1, attains its greatest value, equal 
to 4/27, at the points (u=1/3, k=1) and (u=2/3, k=0). 

Thus, for this worst case, we have 


4h*|f” (E)|max 7A*17” (€)1 max 

R(u kk) max = ——ay = 81 . 
By assuming that max | Ku,k)| = €, we obtain the ex- 

pression for determining the length of the partition seg- 


ment h = Xg-xX, in the form 





3 


h~3,43)/ © _. 
\ Tmax 

The worst case was taken in estimating the approxi- 
mation error. In the general case, a more accurate con- 
sideration is very difficult. Experiment shows that, for a 
wide class of frequently arising functions, we may consider 
that 

3 


ir ®l 


max 
Example. It is required to reproduce the function 


h= 4y 


" f(x) = -10~ in the domain 0 =x <100. We find the law 


of argument partitioning when the approximation error is 
taken as € = + 0.35. 

The third derivative of the given function equals the 
constant quantity 6°10~4. Therefore, all the segments will 
be the same, and equal to 


0,35-104 
3 ’ 29 a 
h,. = 4 \/ errkk ag — 33,6. 


The total number of partition segments is three. 

If the function f(x) is given graphically, then the three- 
fold graphical differentiation necessary for the computation 
of the value of f*""(€) is very difficult, and entails large 
errors. In this case, a graphical method should be used to 
determine the law of argument partitioning. This method 
is based on the replacement of a quadratic approximation 
to the given function f(x) with error € by a piecewise-lin- 
ear approximation to the first derivative f'(x) with error €'. 
The problem reduces to the finding of the relationship €' = 
¢( €,x). Figure 2 shows one segment of the quadratic sec- 
ond-order approximation polynomial P(x). For x = x,, 

X = X, and X = Xg, the function f{ x) to be reproduced and 
the polynomial P,(x) have the common points: 


J (1) = P2 (x), 
-P, (xs), f (a2) = P (as). 

The greatest error of approximation occurs when x =x" 
and x = x'’. It is obvious that for this 


I (2') = Py (x'), ' (2) = Py’ (2’). 
In Fig. 3 there has been constructed the derivative 
f (x) of the given function, as well as the derivative P’s(x) 


of the approximating polynomial. The errors €'y, €', 
and €'s occur respectively for x = Xy, X = Xg, and X = Xz. 


} (%2) 


fur) #,(z) 











Fig. 2 


From the approximation conditions, it follows more- 

over that 
Sanc, = Scope = Sper = Spex. = | 2}. 

At the point x = x’, the angle between f'(x) and P"s(x) 
on the left, a , is not greater than the angle on the right, 
o%: a, =a + Sa, where Aazd. 

By replacing the curve f¢x) by straight lines on the 
segments AC and CE we find, from the equality of the areas 
S AABC™ S ACDE™ | €| , that €'; S's. 

Analogously, we can convince ourselves that €"s 2 ¢'s. 

We now consider F'(x) on the segment x ,x"’. Here 
there is a one-sided linear approximation. Therefore, [11], 
A, (z’x")* 
ee 
where A, is the maximum value of the third derivative 
f***(x) in the interval x’x"’. 

By assuming that point x, divides segment x'x"’ in 
the ratio! , we find that 


Eo 


» A, (2’2*)* 


gx Per, 
° A ‘x"\y? ’ 
0 ES <a (1 
Thus, 


1 1 
e= > Scre=> (Sacer + Sacer + Sazro) = 








= +5 [e+ Mey +4 — 1)? «5] e: 
= == 6; (3/* — 31 + 2). 


From whence 


ar 4e 
2 = ye" (38-31 +2)" 





The expression in the denominator assumes its min- 
imum value for = . With this, €', = 16 €/5a, where 
a=x'x''. This last formula gives the connection between 
the given error € and the tolerance margin €' and size of 
a. 
alll reproducing various functions, the lengths of neigh- 
boring partition intervals differ little from one another. 
Therefore, by measuring the quantity a on a constructed 
segment, we find from the last formula what the tolerance 
will be for the following segment. Inside the tolerance 
band €, constructed with respect tof (x), we carry out a 
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two-sided linear approximation, as a result of which we 
find the size of the following partition interval. 

It fs easily shown that, with such a method 
of partitioning, the actual error will not be larger than the 
admissible approximation error €. Indeed, let the admis- 
sible error €" be computed from the previous segment. To 
determine the size of the partitioning segment, we con- 
struct the tolerance band with respect to f"(x) (Fig. 4). We 
then carry out a two-sided linear approximation. 

Since ACB=s DCE, DFEs_ GPK, thensS BC 
SSAcpDE = € and Sarg = Sapp = & i-€-, with this 
we obtain a certain underestimate of the length of the 
partition interval, which leads only to a decrease in the 
error of approximation of the nonlinear function f(x) to 
be reproduced. 

As the second derivative f"'(x) increases, the curva- 
ture of the function f(x) also increases and, consequently, 
the lengths of the segments of the argument partition de- 
crease. So that the magnitude of the error, defined by 
the previous segment, not exceed the tolerance on the 
given interval, the partitioning should be carried out from 
large segments to small ones, i.e., on the side of an in- 
creasing second derivative f'"(x). 

For the finding of the first segment, we may recom- 
mend a simplified method of successive approximations, 
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shown on Fig. 5 for the function (x) = x*-10~, where the 
graphical construction of the first segment is given. Let- 
ting the magnitude of a, = 20v, we construct the tolerance 
margin €*; = 16 €/5a; = 0.16 €, and we then carry out a 
piecewise linear approximation. We obtain the quantity 
ag > ay = 20. We construct the new tolerance margin 

e's = 16 €/5ag, and carry out a piecewise linear approxi- 
mation which determines ag. It is obvious that the opti- 
mal length of the segment will be ay = ag =a Say. As 
the interval length varies within these limits, the maxi- 
mum actual error €'; decreases from €'; to €'g, varying 
along a certain curve AB which is close to being a straight 
line. At the same time, the admissible error varies along 
a hyperbola which passes through points C and D, which 
correspond to the maximum values of €', for a; = 20 and 
ag = 24. If we now replace the hyperbola by the straight 
line CD we find, on Fig. 5., the point 0 of intersection 
which, with an accuracy sufficient for a graphical meth- 
od, gives the magnitude of the error €' for the first and 
second intervals of the argument partitioning. 


Circuit of the Quadratic Diode Element 





For obtaining segments with quadratic dependencies 
in the present work, we used the diode element represent- 
ed on Fig. 6, a[8]. To the diode are applied input vol- 
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Fig.6 


tage E,, reference voltage e and high-frequency sawtooth 
voltage U,. The diode element's transconductance char- 
acteristic is determined by the conductance Y. 

If the high-frequency voltage U, is not applied to the 
diode, the dependence of the current i on the input quan- 
tity E, has the form shown on Fig. 6,b by the dotted line. 
For E. = F,.° eY./ Y,» diode D is opened, and thereafter 
the current i increases proportionately to input voltage E,. 

If a sawtooth voltage of amplitude A is applied to 
the diode via Y,, the diode will open for lower values 
of the input voltage: 

Y 
E, = E,. =egt— Ay 

By using the principle of superposition, we find the 

average value of the diode element's current: 


0 for E: < Ex, 
i = | a , A 
av rave (Y.Kx— Yee + Ys A)* for Ex > Eng. 


Here, 


st 2 Mie ial 
lng = ey — AZ 


a 
Y,=¥.+Y,.+K +Y. 

For getting a sawtooth voltage, a generator circuit, 
shown on Fig. 7, was developed. A blocking generator 
was executed based on triode V,, while pentode Vs was 
used as a discharge tube. Charging condenser C, is con- 
nected in the anode circuit of Vs. During the blocking 
process, the condenser is charged up to the potential U,.,, 
close to the anode voltage. Tube V, is cut off and con- 
denser C, discharges through tube Vg. Because of the cur- 
rent-stabilizing action of the pentode, the discharge of C, 
occurs linearly. As soon as the voltage becomes equal to 
Uc, = Vunb.v,’ tube Vy opens and the blocking process 















Fig, 7 


repeats. The generated voltage is applied to the grid of 
cathode follower Vg. Part of the output voltage is tapped 
off by the voltage divider formed by resistors Ry, Ry, and 
Ry, and is applied to the cathode of Vg. The magnitude 
of the positive feedback can be controlled by resistor Rg, 
thus varying, within very wide limits, the linearity of the 
output voltage. 

To stabilize the sawtooth amplitude, we use a re- 
ference voltage of +100v, limiting the growth of poten- 
tial on capacitor Cy. The frequency of the generated vol- 
tage is determined by the capacitance of condenser C, and 
by the magnitude of the anode current of pentode Vs. 

With the nominal values shown on Fig. 7, the circuit 
provides a generated sawtooth voltage of frequency 25 kcps 
and output voltage amplitude of 90v. With line-voltage 
variations of + 10%, the sawtooth amplitude varies by + 5%. 


Complete Function-Generator Circuit and 





Its Technical Characteristics 





Vacuum tube diodes type 6Kh2P were used as the func- 
tion generator’s diode elements. The transition to vacuum 
tubes resulted from the great difficulties which arose when 
germanium diodes were used, due to their low back imped- 
ance and their instability under various temperature condi- 
tions. 

The functional schematic of the universal diode-func- 
tion generator is given in Fig. 8. In it are included two 
operational amplifiers Ay and Ay. Applied to the diode ele- 
ments DE are a high-frequency sawtooth voltage and a ref- 
erence voltage from a low-ohmic divider. The low-ohmic 
divider is executed in the form of individual sections, each 
of which consists of six wire resistances: ry = 60 ohm, & = 
120 ohm, tg = 240 ohm, m = 420 ohm, rt = 840 ohm and 
tg = 1680 ohm. For setting up the law of argument par- 
titioning, part of these resistances are short-circuited by 
means of a twelve-position switch with inserts in its panel. 
The total impedance of the low-ohmic dividers is always 
so chosen that a current of 16.7 ma flows through them. 
With this, the increment of the argument between the’ 
switching points of any neighboring diode elements is con- 
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trolled in the limits of 0 to 56v, varying by discrete steps 
of 1 volt. 

In the function generator, there are eight quadratic 
diode elements, for which the switching voltage can be 
controlled, and two linear elements LE, made up of pas- 
sive-resistor circuits, with switching voltages by zero. By 
applying the voltage after the diodes, either directly to A, 
or through sign-inverter Ay, we can obtain increments of 
the function's slopes in all four quadrants. The magnitude 
of these increments can be controlled discretely and smooth- 
ly by, respectively, switch S and potentiometer P (Fig. 8). 

The function generator’s complete circuit also con- 
tains an assembly for establishing the function's initial 
value. 


SUMMARY 


1. The function generator constructed in accordance 
with the principle of quadratic approximations has the fol- 
lowing basic advantages in comparison with function genera - 
tors, in which piecewise-linear approximations are used. 

a) A significant reduction in the number of segments 
of the argument partitioning which are necessary to repro- 
duce the function with the same error. This permits one 
to simplify the universal function-generator’s circuit, to 
decrease markedly the number of vacuum tubes, to ease 

_ significantly the set-up of the functions and the adjust- 
ments of the generator. 

b) The function generator provides smooth repro- 
duction of functions, and eliminates discontinuities in the 
first derivative. 

2. With this, the universal function generator con- 
structed on the basis of the principle of quadratic approxi- 
mation possesses the following disadvantages. 

a) The device's pass band is determined by the saw- 
tooth ‘generator’s frequency which, due to technical diffi- 
culties, can not be made higher than 30 to 50 kcps. There- 
fore, the function generator cannot have a larger band than 

1 keps. 
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b) The sawtooth voltage has a high-frequency com- 
ponent, so that the mounting block requires careful screep- 
ing. 






Despite these disadvantages, the circuit developed fa 
the universal function generator can find practical employ 
ment in electronic models and in various computing de- 
vices. 
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RELIABILITY CRITERIA FOR AUTOMATIC RELAY DEVICES 
WITH RADIOACTIVE EMITTERS 


A. G. Vasil'ev, I. S. Zhitomirskii and K. S. Klempner 


(Khar*kov) 


(Translated from; Avtomatika i Telemekhanika, Vol. 21, No, 2, pp. 245-253, February, 1960) 


Original article submitted July 17, 1959 


Criteria are given for the reliability of automatic relay devices with sources of radioactive emission, The 
relationship is given connecting the parameters of the relay devices with the probability that a relay remain in 


Calculated relationships were given in[1] which al- 
low one to choose the minimum necessary activity of 
emitters in relay devices from the given limits on the ad- 
missible number of standard deviations of the operating 
utilization of the detector from the threshold level. In 
the work just cited, one limited oneself, in dealing with 
the statistical characteristics of the process, to an approxi- 
mate representation of the potential distribution at the in- 
tegrator relay’s output as a Poisson distribution, or as a 
Gaussian distribution close to a Poisson one. The recom- 
mendations based on these considerations turn out to be 
of practical utility in determining the parameters of or- 
dinary industrial relay devices. However, in mawyy relay 
devices in which high reliability or speed of action are 
required, the accuracy of the computations based on the 
methodology of [1] may turn out to be inadequate. 

Moreover, in the methodology given, the reliability 
of a device's operation was given only by the probability 
of a relay's remaining in the proper states. For many de- 
vices, in addition to this criterion, it is of great practical 
value to use, as a reliability criterion, the mean number 
of "false" transfers to be expected in unit time. 

In the present work we render more exact the limits 
of applicability of the first approximation of [1] for var- 
ious values of relay parameters. Methods are given for 
the numerical calculation of reliability to any desired de- 
gree of accuracy. Reliability is defined in terms of the 
probability of “false” relay transfers. Methods are given 
for the computation of the mean number of "false" trans- 
fers in unit time to any desired degree of accuracy. Grapl 
are given for the determination of reliability from the 
probability of remaining in a “false” state, and the num- 
ber of “false” transfers for various combinations of the 
parameters of relay devices. 

As was shown in [2] and [3] the characteristic distri- 
bution function of the random variable 
V (1) 

ae 
for a Poisson distribution of the frequency of the applied 


signal nt) and reaction function (t) of the device has 
the form: 


T 
9 (h) = exp{ \ [erorr—o 4] n(t) dt \ . (2) 


a given state and, in addition, with the mean number of *false® transfers of the relay in unit time. 


From (2) we immediately obtain the expression for 
the cumulants of the distribution 


T 
k, = | + (7 —t)n(t) dt. (3) 
—0oo 
We shall consider devices in which an RC network 
is used an the integrator. Then, the reaction function 
will have the form: 


t 
OD (i) = e RC, (4) 

In (1), V is the potential at the RC network's output 
and Vy is the increment of potential at the RC network's 
output at the moment when a pulse is applied. In the 
steady state, V» is a constant. During transient responses, 
Vp remains practically constant if the pulse amplitude is 
very much greater than the size of V» over its entire range 
of variation during the transient response. This condition 
holds with sufficient accuracy in all relay devices which 
operate as velocity measurers with linear scales. 

A direct use of (2) for calculating a probability den- 
sity and a distribution function would lead to extremely 
cumbersome computations. We shall consider two series 
expansions for the desired probability density p(x) and 
the distribution function F(x), one of which is useful for 
large values of 

v=nRC, (5) 
and the other for small values of v. 

It was shown, in[2, 3] that, as v--oo, the desired dis- 
tribution tends to a normal one. For computing the cor- 
rection to the normal distribution, Rice suggested the use 
of a series first introduced by Edgeworth (an expansion in 
terms of the derivatives of the normalized density of a 
normal distribution). For the probability density, this 
series has the form: 


po) = defeat) —ar ate (Cah) + 
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and, for the distribution function, 
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In using the Edgeworth series, the greatest difficulty 
is in computing the distribution cumulants in the case 
when n is not constant. If the velocity of counting does 
not depend on time, then 


T »(t—T) RC 
- \ 2 (4) rm rn 
00 


and the solution is essentially simplified. 

We shall consider such relay devices for which the 
relative fraction of the time spent in a stationary state 
is quite large. Then, any consideration of the reliability 
of the relay’s operation can be reduced to a consideration 
of its reliability in the stationary states.* For this case, 
(6) and (7) are written in the form 


p(z) = y= (3) -“ o, (n+ ! 


+ (95) + 
as.) 9,87:10-t@ ©. (4) 


v2 V2 
__ 3,58. (4. 
Ay (7a) * 
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a manok ®, (; 


V2) 
(9) 
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(10) 





where a 


® (z) = —a\e *av (11) 


at 


is a normalized Laplace function, $,(zA/2) is a derivative 
of the error probability integral and 


7.—_— ky (12) 
z= an im 
ky? 
is the multiplicity of the standard deviation. 
It is easily shown [1,4] that 
r=n,RC, (13) 
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ky =nRC, (14) 
puff (15) 
After substituting (13)-(15) in (12), we get 
5 on See (16) 
V2RCn ’ 
where ny is the relay’s operation threshold, i.e., the count- 


ing speed for which both relay states are equally probably. It 
is obvious that for the states in the absence of a controlled 
object on the control line (state 1) the magnitude of n will 
exceed ng, and these states will be characterized by nega- 
tive values of z. For the other state (state 2), the quantity 
z will be positive. 

‘The essential disadvantage of (9) amd (10) is their 
nonuniform convergence. These series obviously have 
poor convergence for large values of z and small values 
of v. Forv =0, the distribution function degenerates 

to the Heaviside unit function n(x): 


0 for <9, 
1 for x>0. 


q(x) = ss 


Therefore, for large z and small v, it is meaning- 
ful to use expansions based on the function n(x). 

One is led to this type of expansion by a method 
analogous to that used by Maslov and Povzner [5]. 

As one may easily convince oneself, the desired 
probability density p(x) is the solution of the equation 





1 
ou u(xr—z’, y) — u(z, y) , hee ’ (18) 
=~ = — v(In x + Ro) a 


0 
with initia] conditions 
u(x, 0) = 4(z) 
and for the parameter value y = 1. 
For y = 1, the distribution function F(x) also coincides 
with the solution of (18), but with the initial conditions 
u(z, 0) = y(z). (20) 


With initial conditions (20), the solution of (18) can 
be obtained either by the methods of finite differences, or 
in the form of a series 


F " = eA [4 (z)] = ge + A[n(z)] + 
+ FAq (al t-.-+ SA" na), 


where the letter A denotes the sisi operation 


(19) 


(21) 





(In 2’ + nz) dz’. (22) 


zx 


1 
Au(z)) = \=S—=)—*©) 
0 


The convergence of (21) is also nonuniform. For small 
x, the series converges rapidly, and it is meaningful to use 
it. For large x, it is more convenient to use finite-differen- 
ces methods. 


* A criterion for reliability of relay operation during the 
transient response could be the probability of just one trait 
fer during the response time. 
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Computing the distribution function by some method o1 
other, we can find the average time that the relay remains 
in a given state. 

In addition, as has been stated, it is necessary to com- 
pute the mean number of "false" transfers per unit time. 
For an inertialess relay, this problem leads to the deter- 
mination of the number of intersections of the current po- 
tential value with the potential value V_ for which the 
relay changes state. For state 1,| V| > | V,|, and the 
problem reduces to finding the number of intersections 
coming down from above (S} ). For state 2, | V|<|V,|, 
it is necessary to compute the number of intersections go- 
ing up from below (S ¢ ). 

The general formula, well known in the literature 
[6], for finding the stated quantity ("mean number of 
blips”) requires, for its practical use, a tremendous amount 
of computation. 

In our case, finding the expression for the mean num- 
ber of intersections poses no difficulty. 

We consider a small interval of dimensionless time 
t =t/RC,(t,7 + Ar). The average number of inter- 
sections of some level x), both coming down from above 
and going up from below, in this interval is an infinitesi- 
mal quantity of the order of Ar. As is known [6], this 
quantity can be replaced, with an accuracy of second-or- 
der infinitesimals in Art, by the probability Pe Ar of one 
intersection going up from below and the probability 
P. Ar Of one intersection coming down from above. Our 
aim is to find the following limits: 
St=lim + Plas, 


At 


' 1 wi 
At-0 At * 


Thus, the quantities Ploa. and Phas can be 
found with an accuracy of second-order infinitesimals in 
Ar. 

Each of the realizations of the process under consi- 
deration increases by a unit jump when a pulse is applied, 
and decays exponentially in the intervals between pulses. 

A realization of the process can intersect the level 
X_ from above to below if, and only if, the quantity x(rT ) 
at time T lies between the limits X, and xe and, in the 
time interval T, T+ln (X/Xg), no pulse is applied. Thus, 


nas . 
Pias= \ p(x)(1—v(s) In 2 \ dz =~ xop (ap) At, 
Zo/ 
y (24) 


where the last equality holds with the accuracy required. 
It immediately follows from (24) that 


S+ = Zop (Zo). (25) 


If we now make the transition from dimensionless 
time to ordinary time, and take (1), (9), and (12)-{16) in- 
to account, we obtain the definitive formula 


S+ = nop (z). (26) 









Intersection of level x» from below upwards is pos- 
sible if, and only if, xfr) lies between xX,~-1 and X» and a 
pulse is applied during an interval less than ln[x /(x-1)]. 
This means that 


x, 


Plae= v(t) At ( p(x) dx + v(t) x 


(x,—1)e4* 
——_ x (27) 
x \ p(x) InZo—j dae = 
x—1 
= v(t) At p(x) dz. 
Whence, sy 
St=v(r) | p(x)dx =v(2)[F (x) — F (ao — 4). 


alli (28) 
For values of Xp > 1 (this corresponds, as is easily 
seen, to ngRC > 1), we can set 


St=y (t) P (Zo). (29) 
After simple transformations, we finally obtaint 
St = np(z). (30) 


Starting from what has been presented, we can write 
the expression for the average blip duration. For state 1 
we get 


F (x) 
6 ae cuneate 
r nop (z)’ (31) 
and for state 2 
, 1 — F(z) 
Y ame ‘Geseeeee 
T apis) ° (32) 


It is of interest to find the expressions for the distri- 
bution of durations in “false* states and for the distribution 
of the number of “false” intersections. This is a compli- 
cated problem, and will not be treated in the present work. 


At the same time, for example, it is easily seen that, with 
the condition 


i i 
gotta ae t (33) 


the distribution of the number of “false” intersections is 
close to a Poisson distribution. 

We now consider concrete results of computations 
based on the method presented here. 

There is great interest in determining the quantities 
1— F(z), F(z), S+, St, T+, Tt fora wide inter- 
val of values of v and z. 

Of practical interest are the solutions for large val- 
ues of z, corresponding to automation and control sys- 
tems in which “false” transfers lead to very undesirable 
consequences, and for small v. Such solutions are ob- 
tained effectively by means of computers, using the 
method presented above [(17}-{22) ]. 





t For reliability characteristics based on the number of 
“false” intersections, one can also use functions of S which 
take on values from 0 to 1, for example, P(S) = 1/ (1+S). 
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500 


1000 =1500 §=2000 
Fig. 1. Values of the probability of 
the potential remaining in a “false” 
state. 


For widely used relay devices working with low val- 
ues of z and large values of v, one can use the Edgeworth 
series (9) and (10) for the direct computation of the quan- 
tities listed above. The values of ©, (z/ V2) are com- 
puted by means of Hermite polynomials: 

2 z* 


Zz 

®,, { 75) Vane? Hnrlz). (34) 

Figure 1 shows the values of the quantities F(x) and 
1-F(x) for several values of z and n RC « 2500. These 
quantities equal the mean relative time the potential re- 
mains in a"false” state for, respectively, state 1 and state 
2. The comoutations were made with all the terms actual- 
ly written above in (10). For convenience in comparing the 
results with those previously obtained [1], the product n RC 
has been laid out along the axis of abscissas, being defined, 
as follows from (5) and (16), by the expression 


noRC =v 4+ 2/2. (35) 


The first two terms in(10) do not depend on v, and 
give a solution in a Gaussian approximation (v = a). As 
is clear from Fig. 1, in determining the relative time 
spent in a “false” state, the deviation from the Gaussian 


n—l 


F142 © 
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approximation is the larger, the larger is z and the small- 
er is the product ngRC. 

It should be mentioned that, for state 1, the compu- 
tations carried out in the Gaussian approximation lead 
to an increase in the actualreliability of state 1 in com- 
parison with the given value. With this, the value of the 
minimum necessary activity of the radiation source of 1 
is somewhat overstated. In connection with this, as is 
clear from Fig. 1, the reliability of state 2 will be lower 
than that which follows from the computations in the 
Gaussian approximation , from whence the value of the 
limiting admissible activity of the radiation source turns 
out to be somewhat overstated. However, in devices in 
which the controlled object provides a significant atten- 
uation of the emission, the magnitude of the limiting ad- 
missible activity significantly exceeds the value of the 
minimum necessary activity, as is clear from the calcu- 
lations made in [1]. For these devices, one need not take 
into account, in calculating source activity, the devia- 
tions from the Gaussian approximation for state 2. 

The deviations changing, by an order of magnitude, 
0.5- @(z) occur at Z =+2 for ngRC <10, atz =:3 for 
ngRC <30 and at z = +4 for mRC <50. Therefore, the 
Gaussian approximation, as was already noted [1], en- 
compasses a wide circle of problems, but is inadequate 
when one deals with devices designed for high relia- 
bility (large z), high speed of response (low RC) in con- 
junction with high sensitivity of the relay to radiation 
(small ng). 

At z = -6, the magnitude of F(-6) is (0.2 to 0.5)107 
for values of ngRC from 500 to 2500. In calculating the 
value of F(-6) for nsRC< 500, the number of terms written 
in (10) turned out to be inadequate. 

The table gives the upper limits on the error of com} 
puting the quantities 1-F(x) and f(x), using all the terms 
actually written in(10). The errors are stated to within 14. 
As is clear from the table, the error in computing the rela- 
tive time in a “false” state falls with decreasing v and 
with increasing z. 







































The accuracy of the determination of the number of 
"false" transfers is determined by the accuracy of deter- 
mining the quantity p(x), as is clear from (26) and (30). 
The accuracy of computing the quantities Stand S4 varie 
in approximately the same manner as that for the relative 
duration of the relay in a “false” state (Cf. the table). 







For various values of z and RC, Fig. 2 shows the num 
ber of “false” intersections of potential per unit time as 
a function of ng. At z = 6, the number of terms of (9) take! 
into account permit the computations to be made only for 
relatively large values of ngRC. For example, the compt- 
tations give S#® const © 2°10 for z = 6, RC = 50 and 
Ng < 200. 










As is clear from Fig. 2, the number of intersections 
S | and S # do not equal each other for the same multipli- 
city of the standard deviation. This is particularly appat- 
ent for large values of z. 
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TABLE Error in Computing the Relative Time the Potential 


Error less than 4. % 
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Fig. 3. Average time of "false" intersections. 


Figure 3 shows the average time the potential remains 
in a “false” state for various values of z, ng and RC, calcu- 
lated by (31) and (32). 


It should be mentioned that the inertia of the electri- 
cal (for example, lamps or bells in the anode circuit of 
power amplifier relays) or electromechanical relay output 
(electromechanical relay) can lead to a certain lack of 
coincidence of the relative time the relay contacts (lamps, 
bells) remain in a “false” state as compared with the val- 
ues Computed by the method presented above. 


On the one hand, a sufficiently inertial output at relay 


Operation will not react on intersections less than some val- 


ue Tp, determined by this inertia. On the other hand, re- 
lays with long release times will lead to a drawing out of 
the time that the relay device remains in a “false” state 
after the “false” intersection of the potential has terminated. 
With respect to “false” relay transfers, it may be re- 
marked that their number can be somewhat less than the 
number of “false” intersections of potential, due both to 
intersections lasting less than Ty and to the appearance of 
more than one intersection while the contacts are to be 
found in a “false” state. These circumstances act to in- 
crease the relay's reliability, defined in terms of number 
of intersections, but limit its speed of response to a well- 
known extent. 





Based on the results given above, one can calculate 
the relative duration of relay contacts in “false” states 
and the number of “false” contact transfers for concrete 
relays with known characteristics. Since the number of 
known relay types is very large, this question requires 
special consideration. 

In conclusion, the authors wish to thank L. K. Tato- 
chenko for the great interest shown in the present work, 
and for a number of very useful comments. 
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ESTIMATES, BASED ON USAGE, OF THE SERVICE LIFE OF 


CONTACT APPARATUS 


D. A. Abdullaev 
(Tashkent) 


(Translated from: Avtomatika i Telemekhanika, Vol. 21, No, 2, pp. 254-259, February, 1960) 


Original article submitted April 15, 1959 


The region of useful life of contact apparatus is determined as a function of its modes of operation and of 
the given admissible number of transfers, Starting from the service life, we estimate the increment of use ex- 


penditure per element, 


Hitherto, relay contact apparatus has been widely used 
in systems for automation, remote control, communicatioa, 
etc. The breakdowns (out-of-order conditions) of individ- 
ual elements and blocks¢ which are unavoidable under con- 
ditions of use require significant time for the detection and 
replacement of the faulty element. As a rule, the system 
must stand idle during this. The system's idle time can be 
reduced to a minimum by the opportune determination and 
replacement of the element or block at fault. On the other 
hand, this measure eliminates the possibility of faulty sys- 
tem operation. 

It is well known that the probability of failure of the 
electromechanical (receiving) portion of a relay element 
is very much less than the probability of failure of its 
contacts [1]. This allows one to estimate, with a definite 
degree of accuracy, the useful state of an element from 
the probability of failure of its contacts. 

A general methodology for estimating the reliability 
of operation of contacts in a given circuit as a function of 
such factors as contact pressure, possibilities of sparking 
or arcing, changes in the form of the contacts’ surfaces and 
other physical reasons for their breakdown are given in[1, 
2}. The essential obstacle to using this methodology in 
usage computations is that it is very difficult, in practice, 
to determine the data necessary for the computations and, 
consequently, for the estimation of the state of the con- 
tracts during their usage. 

The method, described by N. Taylor [3] ( the marginal 
checking method), for inspecting allows one to detect ele- 
ments which are on the point of breaking down and to re- 
place them. However, this method has been developed as 
applied to electronic apparatus. Its practical use for the 
estimation of contact apparatus is not worthwhile, since it 
entails superfluous operations. 

The common feature of these methods is that testing 
and verifying of element states must be carried out direct- 
ly at the location where the elements are mounted or by 
removing them which, in many cases, can turn out to be 
impossible or inadmissible in practice. 

It follows from this that it would be of interest to de- 
velop a method for the indirect estimate of the state of 
contact apparatus. Such a method, which would permit 
one to determine, if only approximately, the length of 


serviceable life a given element, would make it possible 
to replace elements in good time, without detriment to 
the system's operation. 

In actual systems, the modes of operation of different 
elements and blocks are diverse. Therefore, their break- 
downs do not occur simultaneously. If the conditions im- 
posed are such that a given system must operate properly 
for a period of H years, then the replacement of elements 
or blocks before the expiration of this period entails an in- 
crease in the material expenditure, a"use"expenditure in 
the given case. It is of interest to estimate the variation 
of these latter expenditures. 


Admissible Duration of Continuous Opera- 





tion of Contact Apparatus 





In order to replace an element at the breakdown 
point, it is necessary to know beforehand what the admis- 
sible duration of its operation is in the given mode of work. 
It is well known that, under given conditions and modes 
of operation, an element's useful service life is determin- 
ed by the number of transfers of its contacts[1]. Since 
failure of these latter is a random phenomenon, one of the 
principal problems amounts to determining the distribution 
law for the probability of failure of the contacts of an ele- 
ment under given conditions. We note that certain authom, 
for example, V. A. Moroz [4], in their works start with the 
assumption of a normal distribution law. 

In this work we shall adopt the convention that, on 
the basis of tests and statistical data on the operation of 
elements whose mode of operation is given by the para- 
meters and structure of the scheme under consideration, 
we shall have determined the distribution law of the prob- 
ability of contact failure, and that from this we shall have 
found the admissible number of transfers of the latter, 
Nadm- 

If we start from this condition, then we can deter- 
mine the admissible duration of continuous operation 


+ By the term “block” in the present paper, we under- 
stand a set of elements implementing a definite func- 
tion and which are found under identical operating con- 
ditions. 
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of an element for a given speed of transfer of its contacts 
during a ca}l cycle. For electromagnetic elements, this 
is 


N 


where v is the speed of contact transfers, in units of sec™! 
Definite types of reed frequency elements in the pro- 
cess of signal reception, transfer the contacts in each per- 
iod of current (voltage) oscillation, coinciding with their 
natural frequency (Cf., for example, [5]). If oscillations 
of frequency f, to which the frequency element is tuned, 


exist during time t, then the total number of contact trans- 


fers during this time will be ft. Consequently, 


N a3 
7 , —. _ adm 
ladm= 3600 fx (2) 
where r is the period of a call cycle, cycle, in seconds ¢. 
Thus, the admissible duration of continuous operation 
of a reed element depends both on its natural tuned fre- 
quency, and on the time during which the signal exists. 


Domain of Proper Operation of Contact 


Apparatus 


If a given circuit element operates N,, times during 
one call cycle with transfer velocity v, then the duration 
of its continuous operation with transferring contacts will 
equal 





. Ne (3) 
Te = 36000 
We assume that n,, call cycles are made during one 
hour. Then, the time of continuous element operation 
during one hour is defined by the expression 


R= thc m, (4) 


: ~ 3600 » 

Multiplying (4) by the number of hours in a year, we 
obtain the time of proper element operation during one 
year: 

T =nNe. 8760 hr (5) 
Y =.3600 v ° 

From the design conditions, the useful life of the sys- 
tem must be not less than H years. During this period, the 
element under consideration must be found in the mode 
of continuous operation during 


HnyNc 
3600 v 
where H is a dimensionless quantity expressing a quantity 
of years. 

With the condition that 


Tad 

Hr > 1, (7) 
the given element operates properly until the end of the 
proposed time of proper functioning of the entire system. 


It is therefore not necessary to replace it. Otherwise, re- 
placement would be necessary. 





T’ =Hhy = 8760 hr, (©) 
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By substituting (1) and (6) in(7), we get 


N am 





8760 Hm, > 1. 7) 
Whence, 
N adm - 
m™ S g760 Ht (8) 






The left member of (8) is the number of transfers of 
the element in one hour. We call it the hourly intensity 
of the element's operation, and we use the notation 







Q= np N.- 





Then, (8) can be given in the form 





Nadm 8 
Q S 8760 H ° (8) 






Thus, we have established the possibility of proper 
element operation with a given N.4_, as a function of the 
intensity of its operation. The appended figure gives the 
curve for Q = ¢(H), constructed from the condition of 
equality of the left and right members of (8"). It is there 
fore a boundary line: the region below it characterizes 
proper operation, and the upper region characterizes im- 
proper functioning of the element during the given per- 
iod H. It should be mentioned that the character of vari- 
ation of the curve Q = 9(H) coincides with the character 
described by Feugel for the experimental reliability char 
acteristic as a function of time [6]. This confirms the 
correctness of the initial premises of our analysis. 

In a similar way, we can obtain the expression for 
estimating the useful life of reed frequency elements 
(Cf. the Appendix) 

















Q@< aan (+). (9) 


It follows from this that, to increase the useful life 
of a reed element, tuned to frequency f, it is necessary 
to shorten the signal's transmission time t. This is pos- 
sible only if the speed of response of the element itself 
is increased. 

Estimating the usage state of contact apparatus by 
the method presented here is most effective in the case 
of mass-produced devices, since the value of N » a 
tablished just once on an experimental basis for individud 
elements and blocks, can be used as a basic parameter 
for all the apparatus produced in the same production rut 









Influence of Useful Life of Apparatus on 











Its "Use" Expenditure 











The economic efficiency of employing a system is 
usually estimated by the time taken to recoup its costs 


[7] 
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T = (10) 
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tHere we have in mind parallel transmissions of fre- 
quency -indication pulses. 
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Boundary of the proper functioning characteristic 
of a relay contact element. 
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where £ K is the total expenditure and Ey is the annual 
saving. 

The total expenditure, together with the one-shot in- 
itial capital outlay and the wage fund for the personnel, 
includes deductions for amortization and maintenance, 
i.e., operating costs ("‘use” expenditures), a variable quan- 
tity which depends on maintenance and replacement per- 
iods of the system's individual elements and blocks. 

Operating costs per element, in monetary units, we 
denote by K, /A, where A is the total number of contact 
elements in the system. As a percentage of the sum of 

costs,** this quantity.amounts to 
K ,xK (11) 


Ka = A100: 


By substituting the value of 2K from (10), we obtain 


K, 
a aa 
Ka=ama Ty. = 
whence, 
K’,A- 100 
Katy 


The condition for proper apparatus functioning dur- 
ing the time taken to recoup the system's cost is 


l <H,, (12) 


where H,, for a given mode of element operation, is de- 
termined from the proper curve of Q = ¢(H). We assume 
that 

T = H. + AH, 


Obviously, element replacement is necessary if it is 
the plus sign which stands before AH. We now analyze 
how the quantity K a varies with this. 

From the figure, we have 


AQ 1 
ew mee: (18) 
or 
dQ 
Wn =tana, 


Bearing in mind the condition for the construction 
of the curves Q = 9(H), we find from (8") that 


dQ Nadm 


en (14) 
"iw H*:8760 * 
By substituting the value of tana in(13), we get 
Alf w 2, (15) 
Nadm 
If we adopt the convention that the values of K , *de- 


termined from (11') correspond to H, = T then, for f <T, 
the magnitude of the operating costs is expressed by the 
formula 


Kin = — 


AA EH 
soa Fy4 (14 ore 

The relative increment of operating costs per ele- 
ment, as a function of the mode and intensity of its 
operation, is defined by the expression 

K* 

Ar 4 4 22% 8760, 
K, adm 


x — 8760). (16) 








8 = (17) 


We can, by a similar method, obtain the expression 
for the case when reed elements are used 


@=1+4 (+ ) 8760. (18) 


Nadm 





SUMMARY 


i, The curves Q = 9(H) for a system's individual ele- 
ments and blocks allow timely replacement of elements on 
the verge of breakdown, and thereby reduce to a minimum 
the system's down time. 

2. The expressions obtained permit the approximate 
estimation of and, if necessary, the correction of, initially 
assumed values for, the operating costs as applied to appara- 
tus operating in identical conditions. 

The development of a method for estimating the sys- 
tem as a whole is the subject of a special investigation. 

The author takes this opportunity to express his pro- 
found gratitude to B. §. Sotskov for posing the problem and 
for giving his attention to its solution. 


Appendix 
Derivation of Fxpression (9) 





The numbers of transfers of a frequency element dur- 
ing one call cycle is 
c = ft. 
If the given object is called times an hour, then the 
total number of transfers of the element's contacts which 
perticipate in the selection of the given object, is 


Nn = "No =,"pft. ay 
It is obvious that the quantity nt = T), is the dura- 
tion of continuous transfer of the element's contacts during 
one hour. From (19), 
**As in Russian. The correct percentage should be 


Ka‘100 _ publisher. 
AIK 





Ne *h hr. 
f-3600 

The time of continuous transfer of the element's con- 
tacts during one year is 


_ % 
Ty F-58700 br 


th = (20) 


(21) 


Furthermore, 


(22) 


HN, 
r= my = 


8760 hr 


By denoting Q = Nomh = ftm,, and dividing (22) by 
(2), we obtain (9) 


0<7-im (+). 
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A NUMERICAL-TO-ELECTRICAL CONVERTER 


Vol, 21, No. 2, pp. 260-265, February, 1960) 


The block diagram is described of a transducer which converts numerical quantities to proportional values 
of voltage or current, the conversion being based on the principle of an intermediate transformation to time inter- 


vals, The converters operation is analyzed, and some general considerations relating to the design of such de- 
vices are given, 


‘ 


For the purpose of transforming numerical quantities 
to voltages or currents proportional to them, one can, in 
many cases, make very effective use of the principle of 
intermediate transformation of the numbers to time inter- 
vals. Below we describe a device based on this principle. 
Its operation is characterized by this, that in each trans- 
formation cycle one counter is used to develop two suc- 
cessive signals whose durations are such that the time for 
the entire cycle remains invariant while the duration of 
one of the signals is proportional to the value of the in- 
put quantity. This signal is used for modulating a stabil- 
ized dc voltage or current. After being smoothed, the 
duration-modulated voltage (current) is converted to the 
electrical equivalent of the numerical quantity. 

To elucidate the principle of operation, we consider 
the block diagram shown in Fig. 1. It consists of the out- 
put register of a computer (triggers Tj -T},), gates for un- 
complemented code (g) -g,), gates for complemented code 
(g')-g'n). 4 counter (triggers with counter inputs T,-Tp), 
controlling trigger T, key K and smoothing device SD. 

The rectangles with the sketches of the bipolar pulses in 
them in Fig. 1 denote differentiating networks. One of the 
pulses in each rectangle has a stroke through it. This 
means that the pulses of this polarity will not henceforth 
be used. 

The converter's operation is comprised of two cycles. 
In the first cycle, via gates g)-g,, the code for the num- 
ber N to be converted is written in parallel in the counter. 
Pulses continuously applied to the counter’s input with fre- 
quency f are added to the number written there, so that 
the counter is returned to its zero state at some moment of 
time. At this moment, an overflow pulse is emitted by its 
last stage (trigger T,)- The emission of the overflow pulse 
terminates the first cycle of operation. 

At the beginning of the second cycle, the overflow 
pulse reverses the state of trigger (flip-flop) T, opening key 
K. Simultaneously with the key's opening, a signal from 
the differentiating network is applied to gates g’)-g",,, by 
means of which the number N is again written in the coun- 
ter, but this time in complemented code. Input pulses a- 
gain are added into the counter until it is reset to zero. At 
the moment when the counter returns to its zero state, a 
second overflow pulse is emitted, transferring T to its orig- 






inal state. This leads to the closing of the key and to the 
writing of the number N in uncomplemented form in the 
counter. After this, the aforementioned two-cycle process 
is repeated. 

It is easily seen that the duration T of the first cycle 
is linearly related to the number N, the complement of 
the number N, namely 





te (1) 
a 

The equality just written is the mathematical expres- 
sion of a well-known property of an n-place counter, sta- 
ticizing input pulses in the number system, with radix b, 
which has been adopted. This property is formulated as 
follows. If some number is staticized in a counter, then 
the counter will be returned to its zero state after one 
adds to it a series of input pulses, the number of which 
would bring the counter total to b" (i.e., the number of 
pulses is the 10's complement in the decimal system, the 
2's complement in the binary system). This added num- 
ber equals the complement of the given number (i.e., 
its 9's complement in the decimal system, its I's comple- 
ment in the binary system) plus one. Therefore, rt can 
be defined as the product of (N+1) by the size of the in- 
terval between two neighboring input pulses. 

As an example, we consider two special cases, when 
the numerical quantity to be converted, expressed in the 
binary system, has its maximum and minimum values. 

If N=0(N=2" -1), then the number written in the 
counter in the first cycle will equal zero. Consequently, 
the overflow pulse appears at the counter’s output after 

2 input pulses, i.e., in the given case T= 2 /f. If the 
number to be transformed has its maximum value, that is, 
N = 2"-1(N = 0), then all the counter triggers will be in 
their "1° state at the beginning of the first cycle, and the 
application of the first input pulse will lead to the emis- 
sion of the overflow pulse, i.e., for this case, tT = 1/f. We 
may show analogously that for any intermediate value, 
the magnitude of T is determined in accordance with (1). 

The number written in the counter at the beginning of 
the second cycle equals N. Therefore, the duration r of 
the second cycle can be obtained from the expression 


ca tt, (2) 


‘= 
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On the basis of (1) and (2), the time T of a complete 
cycle, equal to the sum rf + Tr, is defined as U 
Uy = => (N (7) 
T= NtNt? (3) aia 2 
On the other hand, it is well known that the addition 
of an n-place binary number N to its complement N re- 


sults in a constant which is independent of the absolute 
magnitude of N, namely: 


Or, taking (2) and (5) into account, 






It is clear from (7) that the value of Uoy, is strictly 
proportional to the converted numerical quantity N. The 
constant term connected with the "1° in the parentheti- 
cal expression (N+1) can always be compensated, either 
by writing the number N-1, instead of N, in the computer's 


























N+ W =2"—1 (4) a shay or by biasing U5y, by the quantity 
Us (2 +2). 
Consaqnentiy, om 4 Depending on the concrete requirements, the smooth- 
T = 7 " (5) ing device can be executed in various ways. In particular, 


it may turn out to be advantageous to use the correspond- 
On the basis of (5) and (2), we can draw the very im- ing filters, special electronic devices or even simply point- 
portant conclusion: in a converter executed in accordance ert instruments indicating the average value of voltage (cur- 
with the block diagram of Fig. 1, the duration of a full rent). A consideration of all the possible practical variants 
cycle (composed of the two part cycles) does not depend 
on the value of the numerical quantity N to be converted. 
With this, each cycle consists of two part cycles; the dur- 
ation of one of them, corresponding to the open state of 
key K, is a linear function of the number N. 


The key K is connected in series in a current genera - ‘ le 
? 
ZI I | [) {1 


UY 


tor circuit, or in parallel to a de voltage generator. In 4 
the first case, there will be at the smoothing device's in- 
put, duration-modulated current pulses and, in the sec- 


ond case, voltage pulses. Depending on which quantity . 
is smoothed, the numerical quantity will be converted 7] Cl] 
either to current or to voltage. For the sequel, this dif- 


on 

ference is immaterial so that, for the sake of definiteness, may z oy 

we shall assume that the number is converted to a voltage. 

Figure 2,a shows the curve of voltage variation at the {? Uy Up 

smoothing device's input as a function of time u = {{t). —— _ $ 
The output voltage U,,,,, of the smoothing device is jit Ue | a8 U, . am oie 

proportional to the mean value of the input voltage. If Y 


we denote the modulated de voltage by U,/2, we can Us 
write 
Uo = ay ; (6) Fig. 2. The curve for u = f(t) and the fre- 


quency spectrum for @ = 0.5. 























of the smoothing device is not part of our present task. It 
is relevant here to cite only several general considerations 
which should be kept in mind in designing it. For this, it 
is first necessary to look into the magnitude of the individ- 
ual harmonic components of the curve u = f(t). 

In accordance with the converter’s principle of opera- 
tion, u = f(t) for a given N is a periodic function. At any 
of its points it satisfies the Dirichlet condition, and may be 
expanded in a trigonometric series. We choose the origin 
of coordinates as is shown in Fig. 2,b. In this case, the 
function is an even one and, when expanded in a Fourier 
series, will have, in addition to the constant term defined 
by (6), only cosinusoidal terms. With this, the amplitude 
of the vy th harmonic U, is computed as the coefficient in 
the Fourier series 

+T/2 
U. == | f (t) cos = vt dt. (8) 
—TrR 


If we bear in mind that the function f(t) has the con- 
stant value of U,/2 in the interval (-r/2+0) =t=(r /2-0) 
and equals zero in the intervals -T/2 <t=(-1r/2-0) and 
(r /2+0) =t= T/2, we can write Expression (8) in the form 

+t/2 
U. = \ cost wae. (9) 


T T 
—t/2 


By carrying out the integration in (9), we get 
. 1 
U, ox 72 sin wy -<-, (10) 
Ty 7 
By denoting the duty factor, equal to the ratio r/T, 
by @ we can write 





U,= To sin xv0. (11) 
gv 


If we differentiate (11), we arrive at the expression 
dU 


7 = U, cos rv). (12) 
Fquating the right member of (12) to zero: 
U, cos rv8 = 0, (13) 


we can find the values of @ which correspond to maxi- 
mum values of U,, from the condition 


k 
w=, (14) 


where k is any odd integer. 

It must be mentioned that not all solutions which 
satisfy (14) should be taken into consideration. In fact, 
it is clear from (2) and (5) that 

N +1 
244 
Since N is included in the interval 





(15) ~ 


0< N <2*—1, 
then 








serene GOR anor, (16) 
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Fig. 3 





























Consequently, only those solutions of (14) have phy- 
sical meaning for which the value of 6, for a given v, 
satisfy (16). In particular, forv = 1, only the one value 
of@ equal to 0.5 is meaningful; for v = 3, the three val- 
ues of @, equal to 1/6, 1/2, and 5/6; forv = 5, the 
five values of @ equal to 0.1, 0.3, 0.5, 0.7, and 
0.9, etc. It is characteristic of all this that there exists 
a unique value of @, equal to 0.5, for which the peak 
voltage of all the harmonics, including the first, is a 
maximum. This can be confirmed by setting 6 = 0.5 in 
(14). After this substitution, we obtain the equation 

v =k, (17) 
which is satisfied for any odd v. Therefore, the compu- 
tations of the admissible pulsations at the output of the 
smoothing device should be carried out for the case @ = 
0.5.which is the most disadvantageous from the point of 
view of the size of the harmonic component. The fre- 
quency spectrum for this case, constructed on the basis 
of (6) and (11), is given in Fig. 2,c. 

Figure 3 shows the envelope of the frequency spec- 
trum for any physically meaningful values of v and 6. 
Here, the quantity U, /o has been laid out along the 
axis of ordinates, instead of U,, and the product v@, in- 
stead of v, has been taken as the independent variable. 
As a result of these replacements, we can use one graph 
to show, in general form, the relationships between or- 
dinal number of harmonic, magnitude of the individual 
components’ amplitudes and the duty factors, which are 
proportional to the number N to be converted. The graph 
was constructed on the basis of (11), which can be written 
in the form 






























































































































































U, sin xv@ 18) 
> =u: ( 
- The harmonic components cause pulsations at the out- 





put of the smoothing device. In its turn, the size of the 
admissible pulsations is closely connected with the accus 
racy of the converter’s operation. Indeed, if it is neces- 
sary to convert a numerical quantity into a voltage pro- 
portional to it, this then means that, over the entire range 
of variation of @, the absolute error A at the smoothing 
device's output must not exceed the quantity which corre- 
sponds to one unit. In other words, 4 must be less than 
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component. In particular, for a binary system, the re- 


quirement on the magnitude of the absolute error, taking 


(6) and (16) into account, is expressed by the inequality 
Uo 

A < gn+i 42 : 

If the attenuation of the constant component is in- 

variable in the smoothing device, then A is determined 

in practice only by the pulsations. With this, and based 

on what was presented earlier, it is not necessary to de- 


sign the device for the entire range of variation of 6 .The 
necessary condition for satisfying (19) in this range is that 


it be satisfied at the point corresponding to 6 = 0.5. 


In practical designing, it should be kept in mind that 


the speed of a concrete converter, for given f and n, is 


almost completely determined by the frequency proper- 
ties of the smoothing device. These, in turn, depend on 
the device's form, circuit and parameters. With this, as 
the maximum number of significant places of the numeri- 
cal-quantity converter decreases, the requirements on the 


permissible pulsations at the smoothing device's output, 


the fraction i/b" of the maximum value of the constant 


(19) 


in correspondence with (19), relax sharply, which permits 





a significant increase in speed of action. The speed of 
the converter can also be increased by decreasing the dy. 
ation T of a complete cycle which, as is clear from (5), 
is inversely proportional to the frequency f of the input 
pulses. But increasing f leads to the use of higher-fre- 
quency elements. Therefore, one should significantly in- 
crease the frequency of the input pulses only after one 
has exhausted the possibilities of increasing speed of ac- 
tion by improving the frequency properties of the sinooth- 
ing device. 

Among the advantages of the converter herein de- 
scribed are its simplicity, its lack of precision elements 
or assemblies, requiring fine adjustment, the ease of 
executing the device with transistor or ferro-transistor 
elements, the absence of special requirements on fre- 
quency stability of the input pulses, or on the duration 
of a complete cycle. This latter advantage is due to the 
fact that any variation of frequency automatically leads 
to an inversely proportional variation both of r and of T 
and, in accordance with (6), thus has no effect on the 
magnitude of U,,,.,. 















AUTOMATIC TELEMETRY SYSTEM FOR OIL-WELL FLOW 
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(Translated from: Avtomatika i Telemekhanika, Vol. 21, No. 2, pp. .266-270, February, 1960) 


Original article submitted October 3, 1959 


The basis for establishing the correct operating con- 
ditions for exploiting oil wells is the measurement, and 
taking account of, the quantity of oil and water extracted 
from the wells. 

An optimal technological mode of operation of a 
well can only be established where there is sufficient 
information concerning the well's output and its varia- 
tions under various operating modes. Variations in the 
output of fluid, observed by measurements, allow the 
necessary measures to be taken to guarantee normal well 
operation. 

In the southern petroleum regions of the Soviet Un- 
ion, flow measurements are made by grouped measuring 
stands. The use of such stands permits a sizeable short- 
ening of the pipelines on the fields. At these stands, 
measurement of well flow is carried out manually. This 
increases the staffs of the fields, does not allow frequent 
measurements to be made (the measurers work on one 
shift) and does not permit systematic observation of well 
flow. Moreover, manual measurement of well flow re- 
tards the transition to complex automation and remote 
control of the fields. 

In the present paper we describe a system of auto- 
matic telemetry of oil-well flow, developed at the Re- 
search and Planning Institute (NIPI) "Neftekhimavtomat,” 

The system developed for the automatic measure- 
ment of oil-well flow permits the measurement of oil 
and water flow individually from a well, with transmis- 
sion of the readings to a control point and their record- 
ing on special cards. This system is one of the basic 
devices (assemblies) of the system of complex automa- 
tion and remote (central) control of oil fields, developed 
jointly with the Institute for Automation and Remote Con- 
trol (IAT) of the AN SSSR and the All-Union Petroleum 
Research Institute. In 1957, on the fourth field of the 
"Ordzhonikidzeneft’ " oil-field management, workers 
of the fields and of the NIPI equipped 24 wells with this 
system, these wells being still successfully exploited to- 
day. 

Figure 1 shows the control desk (a) and a measuring 
stand (b). Each measuring trap (measuring tank) serves 
eight wells (the number of wells can be increased). 

Figure 2 shows the scheme of the system for the auto- 
matic telemetering of oil-well flow based on grouped 
measuring stands. Fach such stand is equipped with dis- 


tributor valves with pneumatic-membrane executive mech- 


anisms (MEM) 8, a manometer 17, and position signaller 
(PS) 21. 


To measure the flow of any well, it is necessary to 
switch it from the common collector to a measuring tank. 
This is implemented by applying compressed gas to the 
MEM of the given well. With this, the well is switched 
from common collector 10 to measuring tank 5, the 

MEM distributor valves are controlled by the pneumatic 
stepping distributor (PSD) which has 16 positions; eight 

of these correspond to a connection of a well to the meas~- 
uring tank, and eight to an empty measuring tank. De- 
tails as to the position of the PSD are obtained by means 
of the pneumoelectric transformer 21. The PSD drive Is 
a bellows mechanism. 

When compressed gas is fed to the chamber of the 
PSD’s bellows mechanism, the PSD is switched from the 
given position to the following one. The supply of com- 
pressed gas to the bellows mechanism is implemented 
by electropneumatic valve 16. When this latter re- 
ceives an electric impulse, it sends gas to the bellows 
mechanism's chamber, thereby switching the PSD. In 
addition to the eight valves which switch the wells 
from the common collector to the measuring tank, there 
is common valve 7 for controlling the measuring tank's 
drain line. The MEN. of this valve is subject to com- 
bined control: by the PSD and by electropneumatic 
valve (EPV) 20. The MEM of the valve on the drain line 
is maintained in the closed state during the period when 
a well is connected to the measuring tank by means of 
the PSD and, after the well is disconnected from the 
measuring tank, by means of EPV 20. To control the 
emptying of the measuring tank, zero-level indicator 6 
is connected in its lower portion and, to signal overflow 
of the measuring tank, trouble indicator 14 is connected 
in. Compressed gas, obtained from the same measuring 
tank, is used to supply all the system's pneumatic ele- 
ments. For the pressure stabilization and cleaning of the 
gas, the system contains pressure reducer 18 and filter 19. 
Measurement of the oil and water levels in the measur- 
ing tank is implemented by level transducer 1 with two 
floats, 3 and 4, inside of which constant magnets are 
placed. One of the floats rests on the “water-oil" divid- 
ing level, the other rests on the oil level. The floats 
freely move along the vertical brass tube (or a tube of 
some other nonmagnetic material) inside of which is to 
be found indicator 2, set in motion by the transducer's 
synchronous motor. 

Control of the automated measuring tank's opera- 
tion, and implementation of the interconnections be- 
tween the stand's individual elements is carried out by 
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relay apparatus, located in the measuring tank's elec- 
trical block 13. At the central (control) point is to be 
found the control console, including program device 25, 
program-device commutator 26, relay apparatus 28, pro- 
viding-automatic system operation, fixed frequency gen- 
erator 29, flow recorder 30, and array 27 of signal lamps, 
keys and knobs. 

The program device automatically gives the time 
for filling the measuring tanks by the output of each 
well in accordance with a previously established central 
control program within the limits of 10 minutes to 4 hrs, 
at 10 min intervals. The time for filling is established 
by the position of a contact pin on the program device 
commutator for each well individually. The program- 


ing device is built around stepping selectors. An in- 
dividual stepping selector is assigned to each measuring 
tank. 

The relay apparatus, consisting of several assem- 
blies, connects the basic blocks and elements of the cen- 
tral control point, and provides automatic cyclic opera- 
tion of the entire measuring system. The output record- 
er is a number-printing device for recording the serial 
number of the measuring tank, the well and the levels 
of oil and water. The electronic frequency generator is 
a two-channel one. There are 12 fixed frequencies in 
each channel. 

In the system, a frequency method is used for selec- 
ting the objects (measuring tanks), and for transmission of 


























ng 


-en~ 


fa~ 


eee 









































































































































































































































































































































quyjog TONUCD [PrIUseD 


yn of 














yoold s xX x yoolg 
uoy] ° a tye yd yoorg 
-dnisjq rere — — 9 — Oo s[eUusqS = - jonuoy 
-e7 ; i P| AUeL Ea 33} enuey 
ura3sAS °| | 
gn —eeroe ¢ Oh ome Oo . 3 } i 
+4 gp mnerem: omen r | 4 
—_—_—S—-s" = == ——_— a = GD G2 GD Ge ec aay ey 
| AIRE aS Gee Pee ot: eee eee “pa r 
eo¢ ! i aoe ¥oauw _— a = tT ] 
Licit™ - a 
_- oy 19 - 
yoorg yore a ord = Spee yoord Buy tea, 
a cone we Ett) -youtas uaa 
auitnd> wed aseay | TI] -picoey 112M 
— -Olg 
3 3 ' 
$ * ° } TT eee an c | | 
‘ | ° | cesta aeeitiacion 1 4 | 
¥ ° 7 occu © em 6 am 0 am 0 an © emoen Oem ee=ee, + 
Il: I]. — o>—o—-e _ call ~ Ee lin clei al 3 
x ° 4 | 
ae eee it; ! { 
ee 1 Le 
‘ 90 ~~ e  H on ot Oy 
Feeeeros 2 os Sop FIT] -e80 ones Sanene 
-39025 ope j +103 
Ln ne “7 = -omy ml 

















w 
fo °) 
cm) 

















halle yoO1g 








Suyxt4 jonuo; 

















- ————} yo01g uo 
oak eee eee -eZIoyINy 























21¥1S 
duryq 
wwe x BuyAr 

















ee AI 
A 


} 

















yoorg yoolg 





Burxt quoul 

















UORIIsOg -omseow 




















poe Kee es ee ap Smee eee eee eee ae 








5 2 


























4oOg 





Suto Msuy 














1U}JOg petjonu0D 























—— — — Zuri 01 peyorias, 
—_— Ii X= x= _JUSUIOINSeOY 
0 0 Ome om Om (DINDT] OU) 
. peddois Butta. 
ame ome om ome (10198910 PINbTT) 
-poddois Butta. * 
- Ardara you yUeL, “q 
doccccccccccccees AIduId yuelL, “8 
sjeuS]s Buyiomsuy “py 
—* — x — 1 —— ¥— x -- OIMSPAF, “OD 
—-—-—-—— -— Bunty dois, *q 
— — — ——, Bur TY 01 youtms, “te 
s[Teus{s purwui0> *¢ 
.[eusts uot eBouewy, °z 
UOTIO@UUOD IUFISUOD *T 


:puese7 














the measurement signals, a pulse-time method of selec - _1. Interrogation as to the state of all the measuring 
tion. tanks (periodically, at short intervals of from 2 to 5 min- 
A given measuring stand is a continuously acting com- utes) with the aim of observing measuring tanks which 
plex system, since each of the 12 measuring tanks compri- are emptied subsequent to a prior measurement. 
sing the system is, as it were, an independent unit, and bas 
a closed operating cycle which is virtually independent of 
the state and operation of the other measuring tanks. 
The system's functioning can be decomposed into _3, Halting the filling of a measuring tank after the 
several operations. time for filling has elapsed. 


_2. Switching an emptied measuring tank, observed 
during an interrogation period, to the filling mode. 
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4. Measuring and recording the results. 
_5. Switching to the emptying mode after completion 
of the measurement cycle. 

Figure 3 shows the block schematic of the central con- 
trol point and one of the 12 controlled points of the sys- 
tem of automatic telemetering of oil-well output. We 
shall consider the operation of the system as shown on the 
block diagram in accordance with the technological cycle 
of measuring oil-well output. 

For convenience, we begin the explanation with the 
interrogation cycle. 

At the beginning of the interrogation cycle, the in- 
terrogation block sends interrogation signals over the line, 
i.e., signals of 12 fixed frequencies are successively sent 
over the line. Each measuring tank has its own fixed se- 
lection frequency. Each frequency selects its correspond- 
ing measuring tank, which sends its answering signal back 
to the central control point. 

If the answering signal “tank not empty” is received, 
the interrogation block ceases sending signals to the given 
tank, and switches to the interrogation of the next tank's 
state. This terminates interrogation of the state of the 
given tank. If the answering signal “tank empty” is re- 
ceived, then the order is sent over the line from the cen- 
tral control point, “switch tank to filling.” At the control- 
led point, by means of the pneumatic stepping distributor, 
the next well in turn is switched from the common collec- 
tor to the measuring tank, and filling of the tank begins. 

The time that each well is connected to the measur- 
ing tank, i.e., the time alloted for filling, is established 
beforehand on the program device by the dispatcher (cen- 
tral controller). After the filling time established for the 
given well has elapsed, the program device sends a signal 
over the line to halt filling of the tank. This signal, by 
means of the pneumatic stepping selector and valve, ef- 
fects switching of the well from the measuring trap to the 
common collector. The flow of liquid from the well to 
the measuring tank ceases. The filled tank awaits its turn 
to be measured. With this, the process of preparing for a 
flow measurement terminates. 


It should be mentioned that, in those cases when the 
well does not supply fluid, although the measuring tank 
was in the fill state for the established time, after the well 








is switched from the measuring tank an answering signal, 
“filling stopped” (no fluid) is sent. This signal is fixed at 
the central control point and, when it again is time to in- 
terrogate the given tank, it is set up to be filled from the 
following well. 


A measuring tank filled with liquid remains in that 
state until the recorder is released from operation. As soon 
as the recorder is freed from the previous measurement and 
if there are tanks filled and ready to be measured, the pro- 
cess of measuring the flow in the next tank begins. The 
measurement and fixing of results of the measurement is 
implemented by one recorder, common to all the tanks en- 
tering into the given system. 


The measurement process is begun by the order, sent 
along the line, "begin measurement.” At the controlled 
point, the transducer is switched in, and the measurement 
process begins. At the end of the process, the recorder, hav- 
ing received the measurement pulses, prints on a special 
card the serial number of the tank, the well and the well’s 
output of water and oil separately. 

After completing the measurement process at the con- 
trolled point, the tank automatically switches to emptying, 
i.e., the liquid in the tank is routed to the common collec- 
tor, and the process of interrogating the state of the next 
tank begins. 

Thereafter, the operation of the system repeats cycli- 
cally. 





In addition to automatic operation, the system is e- 
quipped for manual control, manual calling of tanks and 
out-of-sequence flow measurement of any well, at the 
dispatcher’s option. Both at the central control and at 
the controlled points, there are synchronization blocks, 
controlling the synchronization of system operation. If 
synchronization is disrupted, these blocks signal the dis- 
patcher of the trouble in the system. 

There is a special block for fixing the disruptions in 
the system's operation, notifying the dispatcher of the char- 
acter of the disruption if such occurs, and preventing false 
operations. 

In the apparatus presently used in the system for the 
automatic measurement of well flow, the recorder is a 
specially developed device for the given system. 
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